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Preface 



The series of workshops on "What Comes Beyond the Standard Model?" started 
in 1998 with the idea of organizing a real workshop, in which participants would 
spend most of the time in discussions, confronting different approaches and ideas. 
The picturesque town of Bled by the lake of the same name, surrounded by beau- 
tiful mountains and offering pleasant walks, was chosen to stimulate the discus- 
sions. 

The idea was successful and has developed into an annual workshop, which 
is taking place every year since 1998. Very open-minded and fruitful discus- 
sions have become the trade-mark of our workshop, producing several published 
works. It takes place in the house of Plemelj, which belongs to the Society of Math- 
ematicians, Physicists and Astronomers of Slovenia. 

In this ninth workshop, which took place from 16 to 26 of September 2006 at Bled, 
Slovenia, there were some changes: the date for this year's workshop was moved 
due to participants' other obligations from a customary mid July to September, 
and several members of prof. Sannino's group were present thanks to bilateral 
Slovene-Danish collaboration project. They delivered talks, write-ups of some 
of which you can read in this volume and enriched our discussions. This ninth 
workshop differs from the previous ones in the fact that it is a very short period 
between the workshop and the deadline for sending the contributions for the 
proceedings. Because of this many a participant has not succeeded to send the 
contribution in time. We promise to include those, which were received too late 
to be included in this proceedings, in the next year proceedings. Also the discus- 
sion section, which usually is quite the rich one, is in this time missing several 
contributions, from the same reason - too short time. 

We have tried to answer some of the open questions which the Standard models 
leave unanswered, like: 

• Why has Nature made a choice of four (noticeable) dimensions? While all the 
others, if existing, are hidden? And what are the properties of space-time in 
the hidden dimensions? 

• How could Nature make the decision about the breaking of symmetries down 
to the noticeable ones, if coming from some higher dimension d? 

• Why is the metric of space-time Minkowskian and how is the choice of metric 
connected with the evolution of our universe(s)? 

• Why do massless fields exist at all? Where does the weak scale come from? 

• Why do only left-handed fermions carry the weak charge? Why does the 
weak charge break parity? 
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• What is the origin of Higgs fields? Where does the Higgs mass come from? 

• Where does the small hierarchy come from? (Or why are some Yukawa cou- 
plings so small and where do they come from?) 

• Where do the generations come from? 

• Can all known elementary particles be understood as different states of only 
one particle, with a unique internal space of spins and charges? 

• How can all gauge fields (including gravity) be unified and quantized? 

• How can different geometries and boundary conditions influence conserva- 
tion laws? 

• Does noncommutativity of coordinate manifest in Nature? 

• Can one make the Dirac see working for fermions and bosons? 

• What is our universe made out of (besides the baryonic matter)? 

• What is the role of symmetries in Nature? 

We have discussed these and other questions for ten days. Some results of this 
efforts appear in these Proceedings. Some of the ideas are treated in a very pre- 
liminary way. Some ideas still wait to be discussed (maybe in the next workshop) 
and understood better before appearing in the next proceedings of the Bled work- 
shops. The discussion will certainly continue next year, again at Bled, again in the 
house of Josip Plemelj. 

Physics and mathematics are to our understanding both a part of Nature. To have 
ideas how to try to understand Nature, physicists need besides the knowledge 
also the intuition, inspiration, imagination and much more. These fundamental 
questions also receive quite an attention by the general public — see for exam- 
ple articles on these topics in several science magazines. Among them we should 
perhaps mention the article 'Create your own universe' by Zeeya Merali in the New 
Scientist magazine (Issue 2559, 10 July 2006), mentioning the work of our contrib- 
utor Eduardo Guendelman (see the article by Ansoldi and Guendelman on this 
very topic in this issue). 

The organizers are grateful to all the participants for the lively discussions and 
the good working atmosphere. Support for the bilateral Slovene-Danish collabo- 
ration project by the Research Agency of Slovenia is gratefully acknowledged. 
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1 Child Universes in the Laboratory 

S. Ansoldi 1 and E.I. Guendelman 2 

1 International Center for Relativistic Astrophysics (ICRA) 
and Universita degli Studi di Udine, Udine, Italy 

email: ansoldiOtrieste . inf n . it* 

2 Eduardo I. Guendelman 

Ben Gurion Univeristy, Beer Sheva, Israel 
email: guendel@bgu .ac.il 

Abstract. Although cosmology is usually considered an observational science, where there 
is little or no space for experimentation, other approaches can (and have been) also consid- 
ered. In particular, we can change rather drastically the above, more passive, observational 
perspective and ask the following question: could it be possible, and how, to create a uni- 
verse in a laboratory? As a matter of fact, this seems to be possible, according to at least 
two different paradigms; both of them help to evade the consequences of singularity the- 
orems. In this contribution we will review some of these models and we will also discuss 
possible extensions and generalizations, by paying a critical attention to the still open is- 
sues as, for instance, the detectability of child universes and the properties of quantum 
tunnelling processes. 

1.1 The studies so f ar . . . 

The world Cosmology stems from the Greek word cosmos, which meant beauty, 
harmony, and is the name of that branch of science which studies the origin and 
evolution of the universe. Thus, considering its name and the object of its study, 
it is, perhaps, natural to take a "passive" point of view when dealing with cos- 
mological problems, where we use the word passive to emphasize that our experi- 
ence of cosmology is mainly observational in nature. This may undoubtedly be a 
condition that seems hard to change in practice: after all we are dealing with prob- 
lems, as the birth of our universe and its evolution in the present state, which do 
not appear suitable for a direct experimental approach. On the other hand, we do 
not see any reason why this should prevent us from changing our attitude toward 
the problem, switching from a contemplative to a more active one. In our opin- 
ion, a stimulation in this sense is coming already from the theory which first gave 
us the opportunity to address cosmological problems quantitatively, i.e. General 
Relativity. General Relativity raises for the first time the concept of causality as 
a central one in physics. This means that, taking a very pragmatic point of view, 
we have to admit that only a subset of what exists in our universe can be experi- 
enced/ observed by us. This is not because of our limited capabilities as humans, 

* Mailing address: Dipartimento di Matematica e Informatica, Universita degli Studi di 
Udine, via delle Scienze 206, 1-33100 Udine (UD), Italy 
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but, more fundamentally, because of the restrictions imposed by the spacetime 
structure on the causal relations among objects. At the same time causality also 
brings a challenge to cosmologists in connection with the large scale spacetime 
structure; this is because the simplest models of the universe which are built ac- 
cording to General Relativity and whose late time predictions have a reasonable 
degree of consistency with what we observe, seem doomed to have an initial sin- 
gularity in their past so that field equations break down exactly where we would 
like to set up the initial conditions. This undesirable situation looks even more 
disappointing after the observation that many parameters describing the state of 
the early universe are quite far from the domain of "very large scales" which char- 
acterizes the present observable universe. Let us, for instance, consider a Grand 
Unified Theory scale of 1 1 4 Gev: the universe could then emerge from a classical 
bubble which starts from a very small size and has a mass of of the order of about 
1 Kg (by using quantum tunnelling the mass of the bubble could be arbitrarily 
small, but the probability of production of a new universe out of it would be re- 
duced). The density of the universe would, admittedly, have been quite higher 
than what we could realize with present technologies, but the orders of magni- 
tude of the other parameters make not unreasonable to ask the question: might 
we have the possibility of building a child universe in the laboratory? 

As a matter of fact, a positive answer to this question was already envisaged 
some years ago (for a popular level discussion see [28J). In particular Farhi et al. 
suggested an interesting model able to describe universe creation starting from 
a non-singular configuration and involving semiclassical effects. This proposal, 
actually, leaves some open issues, for instance about the semiclassical part of the 
process and the global (Euclidean) structure of the solution. Although since then, 
a few more proposal have appeared, addressing in more detail qualitative issues, 
it is interesting to observe that most of the problems which emerged in the ear- 
liest formulation are, somehow, still open. It is our hope that the present review 
of the different approaches which have been developed along this interesting re- 
search line, will stimulate to study in more detail and with systematic rigor these 
problems as well as other realistic answers to the above question. In our opinion, 
this question is not a purely academic one, and might help not only to change our 
perspective (passing from an observational to an experimental one) in address- 
ing cosmological problems, but also to shed some light on the importance of the 
interplay of gravitational and quantum phenomena. 

We would also like to remark how, this complementary perspective can be 
considered much more promising nowadays than some years ago, thanks to the 
results of recent observations. These observations are helping us in focusing our 
field of view back in time, closer and closer to the earliest stages of life of our uni- 
verse and are providing us with a large amount of data and information that will, 
hopefully, help us in sharpen our theoretical models. This has already allowed 
tighter constrains on the parameters of models of the early universe, giving us 
the chance for a more decisive attack of most of the still open problems. This will 
be a great help also for a "child universe formation in the laboratory" program; it 
can make easier to identify the fundamental elements (building blocks) required 
to model the creation of a universe that will evolve in something similar to the 
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present one. At the same time, it will help us to narrow our selection of the fun- 
damental principles that forged the earliest evolution of the universe, and, as we 
said already before, strengthen our hope to enlighten a crucial one, which is the 
interplay between General Relativity and Quantum Theory. 

This said, in the rest of this section, keeping in mind the above preliminary 
discussion, we are going to give a concise review of the state of the art in the 
field and to to make a closer contact with some of the models for child universe 
formation; in particular we are going to review some of the existing works on the 
dynamics of vacuum bubbles and on topological inflation (both also considered 
in a semiclassical framework). 

Callan and Coleman initiated the study of vacuum decay more than 30 years 
ago [10 9J; after their seminal papers the interest in the subject rapidly increased. 
The possible interplay of true vacuum bubbles with gravitation was then consid- 
ered [11 25]. More or less at the same time and as opposed with the true vacuum 
bubbles of Coleman et al., false vacuum bubbles were also considered. The clas- 
sical behavior of regions of false vacuum coupled to gravity was studied by Sato 
et al. I33I23I31I27I32I22I and followed by the works of Blau et al. and Berezin 
et al. [5 6|. The analysis in |7| clarified some aspects in the study of false vacuum 
dynamics coupled to gravity; in it, for the first time, the problem was formulated 
using geodesically complete coordinate systems: this made more clear the issue 
of wormhole formation, with all its rich sequel of stimulating properties and con- 
sequences. 

The presence of wormholes makes possible a feature of false vacuum bub- 
bles that is otherwise counterintuitive, which is that these objects can undergo 
an exponential inflation without displacing the space outside of the bubble itself; 
this could seem strange at a first look and is due to the fact that they have an en- 
ergy density which is higher than that of the surrounding spacetime and which 
is responsible for keeping the required pressure difference. Because of this, child 
universe solutions appear as expanding bubbles of false vacuum which discon- 
nect from the exterior region. Apart from the already mentioned wormhole, they 
are also characterized by the presence of a white-hole like initial singularity; the 
simplest example can be obtained by modelling the region inside the bubble with 
a domain of de Sitter spacetime and the region outside the bubble with a domain 
of Schwarzschild spacetime. These two regions are then joined across the bubble 
surface, using the well known Israel junction conditions [21 4J; Einstein equa- 
tions, which hold independently in the two domains separated by the bubble, are 
also satisfied on the bubble surface if interpreted in a distributional sense; they 
determine the motion (embedding) of the bubble in the two domains of space- 
time. Although there are various simple configurations of this system (as well as 
more elaborate generalizations) that are appropriate to describe the evolution of a 
newly formed universe (i.e. they are such that the expanding bubble can become 
very large), these classical models present also some undesirable features. In par- 
ticular it turns out that only bubbles with masses above some critical value can 
expand from very small size to infinity. But then these solutions necessarily have 
a (white-hole) singularity in their past; in fact, for all of them the hypotheses of 
singularity theorems are satisfied. 
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In connection with the restriction on the values of the total mass, the situa- 
tion could be improved in theories containing an appropriate multiplet of addi- 
tional scalars[19jj then all bubbles that start evolving from zero radius can inflate 
to infinity if the scalars are in a "hedgehog" configuration, or global monopole 
of big enough strength. This effect also holds in the gauged case for magnetic 
monopoles with large enough magnetic charge: in this way the mass requirement 
is traded for requirements about the properties of magnetic monopoles. 

A possible connection of this approach with the problem of the initial singu- 
larity appears, then, from the work of Borde et al. [8]: they proposed a mechanism 
which, by means of the coalescence of two regular magnetic monopoles (with be- 
low critical magnetic charge), is able to produce a supercritical one, which then 
inflates giving rise to a child universe. This idea might help addressing the sin- 
gularity problem and in this context it is very interesting the work of Sakai et al. 
(30l : in it the interaction of a magnetic monopole with a collapsing surrounding 
membrane is considered; also in this case a new universe can be created and the 
presence of an initial singularity in the causal past of the newly formed universe can 
be avoided. 

To solve the problem of initial singularity, there are also other approaches 
which make a good use of quantum effects. Needless to say, these ideas are very 
suggestive because they require a proper interplay of quantum and gravitational 
physics, for which a consistent general framework is still missing. This is the main 
reason why most of these investigations try to obtain a simplified description of 
the system by requiring a high degree of symmetry from the very beginning. In 
particular, if we describe the bubble separating the inflating spacetime domain 
from the surrounding spacetime in terms of Israel junction conditions [21 4 [, un- 
der the additional assumption of spherical symmetry, the dynamics of the sys- 
tem is determined by the dynamics of an effective system with only one degree 
of freedom: this is called the minisuperspace approximation; in this framework the 
problem of the semiclassical quantization of the system, even in the absence of 
an underlying quantum gravity theory, can be undertaken with less (but still 
formidable) technical problems using as a direct guideline the semiclassical pro- 
cedure with which we are familiar in ordinary Quantum Mechanics. This has 
been the seminal idea of Farhi et al. Q2) and of Fishier et al. 1141131 . One addi- 
tional difficulty in these approaches was in connection with the stability of the 
classical initial state. Interestingly enough, this could be solved by the introduc- 
tion of massless scalars or gauge fields that live on the shell and produce a clas- 
sical stabilization effect of false vacuum bubbles. By quantum tunnelling, these 
bubbles can then become child universes [ 15J and, at least in a 2 + 1 -dimensional 
example [18J, it has been shown that the tunnelling can be arbitrarily small. 

1.2 ... and their future perspectives 

From the above discussion, we think it is already clear that there are many inter- 
esting aspects in the study of models for child universe creation in the laboratory. 

1 The subject of inflation assisted by topological defects was also studied later in 1 35 1 and 
(26). 
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We would also like to remember how most of these models are based on a very 
well-known and studied classical system, usually known as a general relativistic 
shell [21 4J. The classical dynamics of this system is thus "under control", many 
analytical results can be found in the literature and numerical methods have also 
been employed (see the introduction of 1 1 1 for additional references). On the other 
hand there has been little progress in the development of the quantized theory, 
which still remains a non-systematized research field. We stress how a progress in 
this direction would be decisive for a more detailed analysis of the semiclassical 
process of universe creation. 

Before coming back to the quantum side of the problem, let us first consider 
what could be done on the classical one. We will concentrate mainly on the works 
of Borde et al. |8J and of Sakai et al. 11301 , which suggest many interesting ideas for 
further developments. For instance, it is certainly important to extend the anal- 
ysis in [8], which is mainly qualitative in nature, to take fully into account the 
highly non-linear details of the collision process by means of which a supercriti- 
cal monopole is created (this is certainly instrumental for a quantitatively mean- 
ingful use of the idea of topological inflation). Also the study performed in l30l 
should be extended; to obtain some definitive conclusion about the stability of 
the initial configuration, it is, in fact, necessary to study the spacetime structure 
of the model for all possible values of the parameters; it could then be possible 
to determine if stability is a general feature of monopole models or an accident 
of some particular configurations. From the classical point of view, in both the 
above models another central point is the study of their causal structure; it can 
be obtained by well-known techniques, but, again, a full classification of all the 
possibilities that can arise is certainly required to gain support for the proposed 
mechanisms. Known subtleties which require closer scrutiny (as for example, the 
presence of singularities in the causal past of the created universe but not in the 
past of the experimenter creating the universe in the laboratory or, sometimes, 
the presence of timelike naked singularities) make a discussion of the problem of 
initial conditions not only interesting but necessary, especially in this context 

A suggestive complement to the classical aspects discussed above, is rep- 
resented, of course, by the quantum (more precisely semiclassical) ones, where 
quantum effects are advocated to realize the tunnelling between classical solu- 
tions. If (i) the classical solution used to describe the initial state can be formed 
without an initial singularity and is stable, (ii) the classical solution which rep- 
resents the final state can describe an inflating universe and (iii) we can master 
properly the tunnelling process, then we could use the quantum creation of an 
inflating universe via quantum tunnelling to evade the consequences of singular- 
ity theorems. The construction of proper initial and final states has already been 
successfully accomplished. The stability of the initial classical configuration has 
been, instead, only partly analyzed [15J and it would be certainly interesting to 
consider the tunnelling process in more general situations, where, for example, 
the stabilization can be still classical in origin. Although there is some evidence 
1 30 ] of a general way to solve this issue in the context of monopole configurations, 

2 The proper analysis of the Cauchy problem will, in fact, involve resolution or proper 
handling of these singularities. 
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as we mentioned above, the analysis should be extended to the whole of the pa- 
rameter space. At the same time a complementary possibility is that semiclassical 
effects might stabilize the initial configuration. In particular, closely related to the 
problem of instabilities present in many models, is the fact that the spacetime sur- 
rounding the vacuum bubble has itself an instability due to presence of a white 
hole region (see, for instance, |34[). Also in this context quantum effects might 
stabilize the system and help solving the issue. This approach could require the 
determination of the stationary states of the system in the WKB approximation, a 
problem for which a generalization of the procedure presented in [1J (where this 
analysis was performed for the first time in a simplified model) could be useful. 

Another equally (if not more) important point for future investigations is 
certainly related with the still open issues in the semiclassical tunnelling proce- 
dure. We will shortly discuss this by following, for definiteness, the clear, but 
non-conclusive, analysis developed by Farhi et al. [12]: it is shown in their paper 
that, when considering the tunnelling process, it is not possible to devise a clear 
procedure to build the manifold interpolating between the initial and final classi- 
cal configurations; this manifold would describe the instanton that is assumed to 
mediate the process. According to the discussion of Farhi et al. it seems possible 
to build only what they call a pseudo-manifold, i.e. a manifold in which various 
points have multiple covering. To make sense of this, they are forced to intro- 
duce a 'covering space' different from the standard spacetime manifold, in which 
they allow for a change of sign of the volume of integration required for the cal- 
culation of the tunnelling action and thus of tunnelling probabilities. It would 
be important to put on a more solid basis this interesting proposal, comparing 
it with other approaches which might help to give a more precise definition of 
this pseudo-manifold. In particular we would like to mention two possibilities. A 
first one uses the two measures theory [16]; considering four scalar fields it is pos- 
sible to define an integration measure in the action from the determinant of the 
mapping between these scalar fields and the four spacetime coordinates; there 
can, of course, be configurations where this mapping is not of maximal rank and 
if we then interpret the scalar fields as coordinates in the pseudo-manifold of [12[, 
then the non-Riemannian volume element of the two measures theory would be 
related to the non-Riemannian structure that could be associated to the pseudo- 
manifold. In this perspective, non-Riemannian volume elements could be essential 
to make sense of the quantum creation of a universe in the laboratory and it could 
be important to develop the theory of shell dynamics in the framework described 
by the two measures theory. 

A second one, likely complementary, can come from a closer study of the 
Hamiltonian dynamics of the system. Let us preliminarily remember that the 
Hamiltonian for a general relativistic shell, which we are using as a model for 
the universe creation process, is a non quadratic function of the momentum (this 
comes from the non-linearities intrinsic to General Relativity); this makes the 
quantization procedure non-standard and quite subtle too. Moreover, although 
it is possible to determine an expression for the Euclidean momentum and use 
it to reproduce [2] standard results for the decay of vacuum bubbles (as for in- 
stance the results of Coleman et al. HU) this momentum can have unusual prop- 
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erties along the tunnelling trajectory; some of these inconsistencies disappear if 
we consider the momentum as a function valued on the circle instead than on 
the real line [3 J but further investigations in this direction are required; they will 
likely help us to obtain a better understanding of the semiclassical tunnelling cre- 
ation of this general relativistic system and, perhaps, show us some interesting 
properties of the interplay between the quantum and the gravitational realms. 
In this context, it should be also explored how the Euclidean baby universes |29| 
could be matched continuously to the real time universes and in this way provide 
new ways to achieve spontaneous creation of real time baby universes 

To complement the above discussion, we would now like to provide some 
additional contact points between theoretical ideas and experimental evidence. 
We start considering if all creation efforts might end in a child universe totally 
disconnected from its creator or not. Of course, there is not a definitive answer 
also to this problem yet, since this is tightly bound to the child universe creation 
model. Nevertheless, it it is certainly stimulating to address the question if, in 
some way, the new universe might be detectable. There is an indication in this 
direction from the analysis performed in [24J: here a junction with a Vaidya radi- 
ating metric is employed, so that the child universe could be detectable because 
of modifications to the Hawking radiation. Generalizations that apply to solitonic 
inspired universe creatiorH can be important, especially from the point of view 
of a quantum-gravitational scenario in which the exact and definite character of 
classical causal relations might be waved by quantum effects. 

Other issues that could be tackled after having a more detailed model of child 
universe creation, are certainly phenomenological ones. They would also help to 
better understand the differences between purely classical and partly quantum 
processes, which is also a motivation to consider them explicitly and separately. 
Also the physical consequences of different values of the initial parameters char- 
acterizing the child universe formation process (initial conditions) should be an- 
alyzed^] and in this context we would also like to recall the idea of Zee et. al l20l , 
i.e. that a creator of a universe could pass a message to the future inhabitants of 
the created universe. From our point of view this is can be a suggestive way to 
represent the problem of both initial conditions and causal structure; this could 
be of relevance also for the problem of defining probabilities in the context of the 
multiverse theory and of eternal inflation. 

A final point of phenomenological relevance would be in connection with 
observations that suggest the universe as super-accelerating. This seems to sup- 
port the idea that some very unusual physics could be governing the universe, in 
the sense that standard energy conditions might not be satisfied. In the context of 
child universes creation in the laboratory in the absence of an initial singularity, 
it might very well be that a generalized behavior of the universe to try to raise 
its vacuum energy would manifest itself locally with the creation of bubbles of 
false vacuum (as seen by the surrounding spacetime), which would then led to 

3 It is, for instance, certainly possible to extend the metric describing the monopole, i.e. 
the Reifiner-Nordstrom spacetime, to the Reifiner-Nordstrom- Vaidya case. 

4 In particular different ways of creating a universe in the laboratory could lead to differ- 
ent coupling constants, gauge groups, etc.. 
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child universes. In [17 1 a proposal, based on the two measures theory, to avoid 
initial singularities in a homogeneous cosmology has already been put forward. 
It would then be desirable to apply it to the non-singular child universe creation 
also. 

To conclude we cannot miss to point out how all the above discussion about 
the possibility of producing child universes in the laboratory could take a com- 
pletely new and concrete perspective in connection with the possible existence 
of new physics at the TeV scale in theories with large compact extra-dimensions, 
physics that might become available to our experimental testing at the colliders 
which will shortly start to operate. 
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Abstract. We derive a relation between the three finestructure constants in the Standard 
Model from the assumptions of what we call "multiple point principle" (MPP) and "An- 
tiGUT". By the first assumption we mean that we require coupling constants and mass 
parameters to be adjusted - by our multiple point principle - to be just so as to make sev- 
eral vacua have the same cosmological constants (from our point of view, basically zero). 
By AntiGUT we refer to our assumption of a more fundamental precursor to the usual 
Standard Model Group (SMG) consisting of the Ng ert -fold cartesian product of the usual 
SMG such that each of the three families of quarks and leptons has its own set of gauge 
fields. The usual SMG comes about when SMG 3 breaks down to the diagonal subgroup 
at roughly a factor 10 below the Planck scale. Up to this scale we assume the absence of 
new physics. Relative to earlier work where the multiple point principle was used to get 
predictions for the gauge couplings independently of one another, the point here is to in- 
crease accuracy by considering a relation between all the gauge couplings (i.e., for U(l), 
and SU(N) with N=2 or 3)as a function of a N-dependent parameter dN that is a charac- 
teristic of SU(N) groups. In doing this, the parameter dN that initially only takes discrete 
values corresponding to the "N in SU(N) is promoted to being a continuous variable (cor- 
responding to fantasy groups for N g Z). By an approprite extrapolation in the variable 
dN to a fantasy group for which the |3-function for the magnetic coupling g 2 vanishes we 
avoid the problem of our ignorance of the ratio of the monopole mass scale to the fun- 
damental scale. In addition to increasing the accuracy of our predictions for the gauge 
couplings by circumventing the uncertainty in our knowledge of this ratio, we interpret 
our results as being very supportive of the multiple point principle and AntiGut. 

2.1 Introduction 

In earlier work [1] we invented our Multiple Point Principle / AntiGUT (MPP 
/ AntiGUT) gauge group model for the purpose of predicting the Planck scale 
values of the three Standard Model Group gauge couplings. These predictions 
were made independently for the three gauge couplings. In this work we test an 

* Invited talk by D.L. Bennett 
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alternative method of treatment of MPP/AntiGUT in which we seek a relation 
that would put a rather severe constraint on the values of the SMG couplings. 

An important ingredient for the calculational technique in this paper is the 
Higgs monopole model description in which magnetic monopoles are thought of 
as particles described by a scalar field c{) with an effective potential V e ff of the 
Weinberg-Coleman type [2j,|3j. The MPP is implemented by requiring that the 
two minima of V e f f are degenerate. This requirement results in a relation between 
the square of the monopole charge g 2 and the self -coupling A that defines a phase 
transition boundary between a Coulomb-like phase (with < cf) >= 0) and a phase 
with a monopole condensate (with < cf) >^ 0). For Abelian monopoles this phase 
boundary has been presented in earlier work |4],|5|. This phase boundary condi- 
tion which has a term quadratic in g 2 (as well as terms linear and quadratic in A) 
consists of tuples (A, g 2 ) of critical values of A and g 2 . A characteristic feature of 
this boundary is that it has negative curvature as a function of A and hence a max- 
imum value Sum) crit mQX of g 2 . It is important to keep in mind that this phase 
transition that we find by applying MPP to V e ff is assumed to be at the scale of 
the monopole mass. In the present work, the V e ff is generalized in such a way 
that it also embodies nonAbelian SU(N) monopoles using the technique of Chan 
& Tsou in which a SU(N ) monopole is described as a N under a dual Yang-Mills 
vector potential [6]. The corresponding phase boundary now has coefficients 
that depend on N (i.e., "N" as in SU(N)). But otherwise the phase boundary for 
a SU(N) monopole is qualitatively the same as that for the Abelian monopole the 
important difference being that the maximum value of g 2 in the Abelian theory is 
different than in the nonAbelian SU(N) theory: g 2 (1 , crit mQx + g 2 u(N) crit mQx . 

In fact a major trick in the present article is to consider a correspondance 
between Abelian and non- Abelian monopoles (the latter in the understanding of 
Chan- Tsou to be explained in section 2 below). This correspondence is defined 
using a N -dependent parameter C = C(N) defined by 

□ 2 

£^ » SU(N ) crit max (21) 
»U(1 ) crit max 

We can think of this definition of parameter C as a definition of the Abelian the- 
ory with 9 u nj CTit max that corresponds to the nonAbelian SU(N) theory with 

9sU(N) crit max" 

As shall be seen soon we need an unambiguous way to define the Abelian 
magnetic charge Sum corresp. su(N) corresponding to any nonAbelian SU(N) 
magnetic charge g|u(N)' Our definition of such a correspondence is simply: 

C^ , 2 §iu(N) , (2.2) 

9u(1) coTresp SU(N) 

where gf; U ( N ) can have any value (not necessarily critical or critical maximum). 
We can thereby think of say the fundamental scale nonAbelian dual - i.e. mono- 
pole - coupling also as an Abelian one. 

We now go to the Planck scale where the MPP / AntiGUT model was origi- 
nally invented as a simple way of relating the experimental values of the SMG 
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gauge couplings gu(i ) , 9su(2) an d 9su(3) (extrapolated to the Planck scale in the 
absence of new physics underway) to the critical values of these three coupling 
as determined using lattice gauge theory. This MPP/AntiGUT relation in terms 
of the (critical values of the) gauge couplings gt (i e {U(1 ), SU(2), SU(3)}) is in 
this work reformulated in terms of the dual (critical values of the) of the mag- 
netic charges using the Dirac relations gg = 2n and gg = 4n for respectively 
Abelian and nonAbelian monopole theories. Recall that we already have a con- 
vention for calculating the (squared) Abelian magnetic charge that corresponds 
to a (squared) nonAbelian magnetic charge cjsu(N) usm g me parameter C: 

2 _ 9su (N ) , . 

9u(1) coTTesp. SU(N) ~ q • K*" ) 

Now assuming for the moment that our MMP/AntiGUT model is in fact a law 
of Nature it would not be unreasonable to discuss whether the Abelian corre- 
spondent couplings g u(1 ) COTTesp su(2) coupling is smooth or not as a function of 
a gauge group characteristica such as N for SU(N) groups. Also the U(l) gauge 
group can be taken into consideration using our nonAbelian to Abelian corre- 
spondence relation in the special case 

n 2 = (2 4-1 

9u(1 ) corresp. U(l ) — ~Q~~] \ > 

which will be seen below to correspond to d N =0 

Actually for later convenience we shall instead of N use an N -dependent 

parameter dN as our independent variable and the quantitiy 



j~2 

■JSUJ 1 ) corresp.SUU) 



instead of g u(1) corTesp su(2)) as our on analytically dependent variable. We 

have now three points (d N , 39u(1 " : ° r ; esp - SU(N) ) belonging to our hypothesized 
in d N analytic function namely the points 



(0 3 ^ U(1) ) (d 2 ^Hlj CQTTesp. SU(2) ^ ^ 3 9 U (1) coTTesp. SU(3) ^ ^ ^ 

7T 7T 71 

Now we make the guess that the analytic function in dN that we seek is the 
parabolic fit obtained using these three points. It must be emphasized that our hy- 
pothesized function analytic in dN lives at the Planck scale while the phase tran- 
sition boundary discussed above lives at the (unknown) scale of the monopole 
mass and consists of critical A and g 2 values that both run with scale. So the prob- 
lem is how to connect hypothesized Planck scale physics in the form of our in 

dN analytic function — "'" c °^"p- su ' 2 ) with say critical values of the function 
9u(i ) coTTesp. su(N) cTit at tne unknown scale of the monopole mass. 

There is one value of g 2 for which this connection would be trivial namely 

the special point ( dN s p e c , 3gs 7 ^" !C ) lying at the intersection of our in dN analytically 
continued function with the function defined by requiring that the (3 -function for 

? 3n 2 

g L vanishes. I.e., |3g2 =0. Just this value of — is the same at the Planck scale 
and at the (unknown) scale of the monopole mass. We describe now briefly how 
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we find the "fantasy" (and nonexistent!) gauge group '"SU(N sp e C )"' for which 
the corresponding g 2 does not run with scale. But first a little digression on how 
(3g2 depends on 

Starting from the definition of the Abelian magnetic charge correspondent to 
a nonAbelian magnetic charge we use the nonAbelian Dirac relation to obtain 



9su(N) 1 16tt 2 



9u(1 ) corresp.SU(N)— p ~ p J 

Gel 



9su(N 

We now require that the Dirac relation remain intact under scaling; i.e., 

16tt 2 d , 1 , 16tt 2 1 d , 
- = lc 

JSU(N 

d r~-l 1 — 1 1 N 



(2.6) 



^SU(I) corresp.SU(N) " —^tf— ") " ~ 4^ a SU(N))- ( 2 - 7 ) 



Using that ^(a su(N) ) = \j£ (which is just the usual Yang-Mills contribu- 
tion to the |3-function for <x<T u ( N f using t=lnu 2 ) we get the Yang-Mills contribution 
to the running of g 2 (1) corresp SU(N) : 

^S 2 umcorTesp . su(N) lY.-M. contTll) -— — — -— -d N . (2.8) 

So even though there is of course no Yang-Mills contribution to (3u(i ) we see that 
Pu(i) corresp. su(N) inherits a dependence on P g 2 u(N) through the requirement 
that the Dirac relation remain intact under scale changes. Using the known (3- 
function for § u( ] ) (to 2-loops): 

g 4 g 6 

P 9 2 um ~ 48t? + (16tt 2 ) 2 (2 ' 9) 
which leads to the (3 -function for g 2 (1) corresp SU(N) : 

g 4 Q 6 

Pa 2 = Pa 2 + d N = T^Hr + 7T7-T7 + d N • (2-10) 

0U(1 ) corresp, SU{N) 9u(l) 487t 2 (167t) 2 

The intersection point (with no running of g 2 )is readily found as the intersection 
of 

Pc^ =0 (2.11) 

9 U(1 ) corresp, SU(N) 



and our in dN analytically extrapolated function ^ u(1 1 corresp. su(N) ^ ^ ^ 
intersection point we have 



(d N) — ) = (-0.57,14.54). (2.12) 

7T 

So any intersection point would have a g 2 that necessarily has a RG trajectory 
parallel to the A axis in the space spanned by (A, g 2 ) where the phase transition 
boundary lives. But only one of the horizontal (i.e., parallel to the A axis) RG 
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trajectories can be tangent to the phase transition boundary and the point of tan- 
gency must neccessarily be at the top point of the phase transition curve. The 
interesting result of this work is that our intersection point (which depends of 
course on the experimental values of the SMG finestructure constants) to high 
accuracy is the value of g 2 with the RG trajectory that is tangent to the phase 
boundary at its top point where g 2 — § u m) crtt max = 15.1 1. Had our intersec- 
tion point singled out any other RG trajectory our MPP would have been falsified 
(see Figures |2~T1 and l2~2t . 

2.2 The Chan-Tsou duality and monopole critical coupling 
calculation 

Investigating nonAbelian theories, we have used the quantum Yang-Mills the- 
ory by Chan-Tsou [6] for a system of fields with chromo-electric charge g and 
chromo-magnetic charge g (monopoles). This theory describes symmetrically the 
non-dual and the dual sectors of theory with nonAbelian vector potentials A,^ 
and A H covariantly interacting with chromo-electric ) VL and chromo-magnetic ) yl 
currents respectively. As a result, the Chan-Tsou nonAbelian theory has a dou- 
bling of symmetry from SU(N ) to 

SU(N) x SU(N) 

and reveals the generalized dual symmetry which reduces to the well-known 
electromagnetic (Hodge star) duality in the Abelian case. 

We want in principle to consider three phase transitions connected with a 
single nonAbelian monopole which in the philosophy of the Chan-Tsou-theory 
to be described below is an N-plet under the by Chan-Tsou introduced dual Yang 
Mills four vector field A H , namely 1) a confining phase, 2) a Coulomb phase, and 
3) a phase with monopole condensate. According to the Multiple Point Principle 
the coupling constants and mass parameters should then be adjusted in Nature 
to just make these phase degenerate (i.e. same cosmological constant). Earlier we 
have used two loop calculations to obtain the Abelian gauge group a phase tran- 
sition between the monopole-condensate phase and the Coulomb phase using the 
Coleman- Weinberg effective potential technique, which led to a relation between 
the self coupling A for the monopole Higgs field and the monopole charge g. To 
determine the monopole mass we should, however, in principle involve one more 
phase, i.e. the monopole confining one, but that would need a string description 
in the language used for the two other phases and doing that sufficiently accu- 
rately for the fit towards which we aim in this article is not undertaken here. 
Hence we shall assume only that the ratio of the mass scale of the monopole con- 
densate or approximately equivalently the monopole mass to fundamental scale, 
taken here to be the Planck scale, is an analytical function of some group charac- 
teristic, which we shall take to be the quantity dN that we shall return to shortly. 
The basic point is that we derive by the Coleman-Weinberg effective potential 
an a priori scale independent phase transition curve in as far as the monopole 
mass drops out of the relation describing the phase border between the Coulomb 



2 Relation between Finestructure Constants at the Planck Scale 



15 



phase and the monopole condensate one, so that the only scale dependence of 
this relation comes in via the renormalization group. The lack of a good tech- 
nology for calculating the mass scale of the monopole therefore means that we 
have troubles in calculating the renormalization group correction of the by the 
Coleman- Weinberg- technique calculated relation between A and g 2 to run it from 
the monopole mass scale to the fundamental scale. The major idea of the present 
article now is that this would be no problem if the beta-function for the monopole 
coupling g had been zero. The trick now is to effectively achieve that zero beta 
function by extrapolating in the gauge group so to speak to a "fantasy" group 
having zero beta-function. 

A priori magnetic monopole couplings for different gauge groups cannot 
be compared, and so to make the statement that the phase transition coupling 
is analytic as a function of some group characteristic, call it dN say, is a priori 
not meaningfull. This is so because in principle one could vary notation from 
group to group, and such a choice of notation would not a priori be analytical. 
We shall primarily be interested in the phase transitions from a Coulomb phase to 
the monopole condensed phase as obtained from studying the effective potential 
as a function of the norm of the vacuum expectation value of the monopole scalar 
field in the manner of Coleman- Weinberg. We decide to take the a priori arbitrary 
ratio between the ratio of a gauge coupling for an SU(N) gauge group and the 
coupling for the corresponding Abelian U(l) theory to be the same as that for 
the critical values of these couplings. We take the Lagrangian densities for a U(1 ) 
theory and an SU(N) Higgs Yang Mills theory respectively as 

and 

c = + H\4> Q I 2 - ^ 2 I4> Q I 2 - ^(I4> a l 2 ) 2 (2.14) 

where c() a is a monopole N.-plet, is the covariant derivative for dual gauge 
field and g is magnetic charge, such that the meaning of the mass \x and the 
self coupling A becomes related in as far as we decide to identify as having corre- 
sponding meaning of the length squares of the fields; i.e. we identify 

N 
a=1 

as is natural, since from the derivative part in the kinetic term respectively \ | 3 ^ cf) | 2 
and J 1 5" n't' 1 1 2 we can claim that a given size of |c()| 2 and I4> a | 2 corresponds to a 
given density of Higgs particles, a number of particles per unit volume being the 
same in both theories. Accepting ([2.151 as a physically meaningfull identification 
we can also claim that the A and the \x in ([2.131 and (|2.14) l are naturally identified 
i an N-independent way (i.e., N as in SU(N)). 

Denoting l|2.15t by just |c{)| 2 one can write - as is seen by a significant amount 
of calculation or by using Coleman-Weinberg [2J and Sher [3] - the one-loop ef- 
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fective potential for U( 1 ) and SU(N ] gauge groups as 

v eff = -Vi4>i 2 + ^r 14,14 



2^ 4 |VKI 64n 2 

[ 3Bg 4 ln^ + (-^ + 3A|4 ) | 2 ) 2 ln- ^2 + 3A|ct,|2 



where 
and 



M 2 v ^ Wl ; M 2 
+ A 
M 2 

A = B = 1 for Abelian case, (2.17) 



A(-^ + A|<t>| 2 ) 2 ln ^+ A|4) ' 2 ], (2.16) 



A = 2N-1, (2.18) 

(N-1)(N 2 +2N-2) 
B = ^ for SU ( N ) § au § e group. (2.19) 

The SU(N) formula we used here were derived using an N-plet monopole and 
with a convention for the covariant derivative 

D^ = 9^-iAj 1 AV2 (2.20) 

in the convention with absorbed coupling where the generators A'/2 were nor- 
malized to 

TT(yy) = (2-21) 

while for the Abelian theory we used the convention 

D^ = 3 H -iA^. (2.22) 

We have stable or meta-stable vacua when we have minima in the effective 
potential l!2.16l l which of course then means that the derivatives of it are zero 
there: 

avefftw 2 ) , . ._ (223) 

3|4>| 2 lmxn 1 - ( ' 

where i enumerates the various minima. 

Now our multiple point principle asserts that there should be as many degenerate 
vacua as possible - i.e., the more degenerate vacua the more intensive parameters 
that become finetuned by the requirement of being at the multiple point (in pa- 
rameter space). In the case we consider here there are just two degenerate vacua 
at say <£> = 4) m t n i and 4> = 4) m t n 2. This means that if we take the degenerate 
minima to have zero energy density (cosmological constant) that 

Veff (Mminl ) = Veff (WSvtra) = 0. (2.24) 

The joint solution of equations ([2.241 and (12.23ft for the effective potential l!2.16l l 
gives the phase transition border curve between a Coulomb phase and monopole 
condesced phase: 

3Bg£ t + (5 + A)A 2 t + 167t 2 A p . t . = 0. (2.25) 
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All of the combinations (A, g 2 ) satisfying (12.25ft are critical in the sense of separat- 
ing phases. The maximum value of g u(1 ) CTit - we have called it g^, crU max 
turns out to be interesting for us. Let us now find this top point of the phase 
boundary curve (|2.25b (see also |4] and |5]). 



d§4 | ■ -o 



which gives 
and 

From Eq. J2.27b we obtain: 
for U(l) group: 
A=l, B=l, 

for N=2: 
A=3,B = ^, 



~4 

9crit 



*crit 



(16tt 



2\2 



12(5 + A)B 

]6n 2 

~2(5 + A)' 



9l 



rit,U(1) 



8n 2 
3V2 



18.61, 



9crit,SU(2) 



8n 2 
3V2 



37.22, 



for N=3: 
A=5, B = 



11 

36 ' 



~2 
9cri 



^108 8tt 2 



rit,SU(3) y g 5 ^ 

In general we define the parameter C such that 



ll. 99. 



gcrit,SU(N) = C 9crit,U(l)> where C 



(5 + A)B 



(2.26) 
(2.27) 
(2.28) 



(2.29) 



(2.30) 



(2.31) 



(2.32) 



We shall assume that this relationship between the Abelian and nonAbelian cou- 
plings is also valid when the couplings are not critical. These results are given at 
the scale of monopole mass or VEV. 
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2.3 Compilation and correction of finestructure constants in 
AntiGut model 

Recall that MPP is to be applied to the N g en -fold replication of SMG. For Ng en = 
3 we have 

(SMG) 3 = U(1 ) 3 x SU(2) 3 x SU(3) 3 (2.33) 

that breaks down to the diagonal subgroup at roughly the Planck scale. It is the 
couplings for the diagonal subgroup that are predicted to coincide with experi- 
mental gauge group couplings at the Planck scale [1J: 

a 1 ,exp = 6 a l CTi . t , 
a 2,exp = ^ a 2,crit> 

"iT.lxp = 3a 3, 1 crif ( Z34 ) 

According to the Particle Data Group results [7J, we have: 

oT.Uptepi) « 55A > «I,U^ « 49.0 3 ; ot^ exv [ii n ) « 53.00. (2.35) 
In the Abelian theory the Dirac relation is 

gg = 27m where n € Z (2.36) 

which leads to 

4 



act = \ (2.37) 



or 

a 2 

a" 1 = 4a = — . (2.38) 
In the nonAbelian case the Dirac relation is 

gg = 47tn where n £ Z. (2.39) 

This leads to 

eta = 1 (2.40) 

or 

a" 1 = a = |- (2.41) 

47t 

From Eqs. < |234l and ll2~35l l we have at the Planck scale: 



7T 



(,) -27.7, 



3 9sU(2)/Z 2 



7t 

=2 



196, 



3q2 

y S u ( 3)/z3 w212a (2 42) 



7T 
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But the correction of the running of finesructure constants from AntiGUT group 
l|2.33b in the interval At = \/40 gives: 



«J,U ^pO « 53 - 3 '. a J,l X p ( Wn) ra 59 - 4 - ( 2 - 43 ) 



T — 2 

The corresponding Abelian values of are: 
for U(l): 



for N=2: 



for N=3: 



3 9up ) 



7T 



27.7, (2.44) 



3 9 UI 



71 



106.6, (2.45) 



3 9 U (i) _.,.„„ /65 



59.4 • W— « 184.32. (2.46) 



7t ' V 27 



2.4 The d N -parameter. 

As we suggested in the introduction we shall consider a correspondance between 
the nonAbelian and Abelian Chan Tsou monopoles and in this connection have 
a correspondance of couplings so that the ratio of corresponding monopole cou- 
plings is like that of the critical couplings for the same two theories; i.e., as in 
equation d2.32| |. Since we want to avoid having to try to use our bad knowledge of 
the ratio of the monopole mass scale at which the phase transition couplings are 
rather easily estimated and the Planck or fundamental scale, we are interested in 
beta-functions. So the most important feature of the gauge group for the purpose 
here is how the monopole coupling will run as a function of the scale. For this pur- 
pose we use |3g2. Now if we consider - as is the simplest - an Abelian monopole 
we strictly speaking would have no group dependence of the beta function (3 § 2 , 
but now we imagine there there actually is an effect of the selfcouplings of the 
Yang Mills fields in the "electric" sector and include that. We shall do that by 
the assumption that there will be a term corresponding to it in such a way as to 
make the Dirac relation (or its replacement for a non Abelian monopole) | |2.37| | 
be valid for the running couplings at all scales. In J2.32|l we have the relation be- 
tween the Chan-Tsou SU(N) coupling g| U ( N i for the critical coupling, which we 
here by definition of the relation between corresponding theories extend also to 
non-critical couplings: 



9 S U(N) = Cg u(1 ) corresp SU(N). where C = J (5 + A)B • ( 2 ' 47 ) 

In the Chan-Tsou formalism the replacement for the Dirac relation is that 



gg = 47m, neZ. 



(2.48) 
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Combining | |2.47|| and j2.48|l and taking n to be unity we get for the Abelian mag- 
netic charge corresponding to that of a nonAbelian SU(N) 

n 2 _ dfat 2 ) (249) 

9U(1) cott SU(N) — p 2 ' 
L 9su(N) 

With the postulate - but that is really true - that the running of the couplings shall 
be consistent with the Dirac relation we can take the scale dependence of this 
equation 112.491 on both sides to obtain 

|3 5 2 =-(16tt 2 /C)-R 2 /g| u(N ,. (2.50) 

Su[l] corresp SU(N) USUIN) 1 ' 

Now there is the group dependent contribution to the well known beta function 

ft I _4 UN 

Pg| u(N) lY.M.coixtTibutioix - 9sU(N)'4g^2 V" OL ) 

which to keep the Dirac relation valid at all scales must be tranfered to also exist 
in the beta function for the square of the monopole charge 

Pgily.M. contribution = (167t 2 /C) ■ = ~$q~ = ^N (2-52) 

2.5 Our monopole coupling versus d curve and successful 
agreement 

Let us make it quite clear what we mean by our intersection point and our so called 
fourth point. 

Our intersection point is the point at which two functions that live in the 
space spanned by the variables (dN, g u(1 j) intersect one another. One of these 
functions is our in dN extrapolated curve of "experimental" values of function 

9U(1] corresp SU(N ) • The ° ther function Consists Of values (d N , §u( 1 ) corresp SU(N] 

that satisfy the condition 

a 4 a 6 

Pa* =13,2 + d N = -f-j + -f-y + d N = 0. (2.53) 

B U(I | corresp, SU[N] H U(1) 4o7t (lo7t) 

The value of g u(1 ) corTesp su(N) at tne intersection point remains the same 
under scale changes of course since its (3 -function vanishes. 

Our fourth point - we could denote it as (dN 4th .gu(i) crit max) " uves m me 
space of variables (dN , g 2 )- By definition the second coordinate is the maximum 
value value of cj 2 ^ 1 crit on the phase transition bounday which lives in the space 
of the variables (A, g 2 ^ , ). We have above denoted this maximum alias top point 
by the symbol §um crit max- ^ ne nrst coordinate of the fourth point - i.e., dN 4tH 
- is the value of dN obtained when g u(1 j crit ma> , is substituted into Eqn 1 12.531 ) 
above. We get for the fourth point 



(dN,8 U (i) crit ma J = (-0.62,15.11). (2.54) 
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The success that we have in this paper is that the intersection point coincides 
with the fourth point to very high accuracy. 

Maybe it's instructive to think of choosing a bunch of §uni values satisfy- 
ing Eqn l|2.53b (corresponding of course to a bunch of different cIn values). With 
this bunch of §m-n values we know how to RG run them back and forth be- 
tween Planck scale and monopole mass scale - also in the space spanned by (A, g 2 ) 
where the phase transition curve is located with its top point 

(Acrtt, 9u(1) crit max) = (-7-13, 15.11). 

This bunch of cjy , , > values don't RG run at all by definition so they must be par- 
allel to the A axis in the space in which the phase transition boundary lives. What 
happens to these parallel to A RG tracks of the bunch of g~u (1 ) values depends 

only on the height (i.e., ) cr it max) anc ^ not on me details (including the un- 
known scale) of the phase transition curve. Our intersection point follows the 
one horizontal RG track that can become a tangent to the phase transition bound- 
ary and the point of tangency is necessarily the top point. RG tracks below the 
one corresponding to our (fortuitous for MPP) intersection point would hit the 
phase boundary below the top point and therefore corresspond to having only 
a monopole condensate phase. And being within the condensate phase and re- 
moved from the phase transition boundary so that the Coulomb phase is energet- 
ically inaccessible would violate MPP. The (horizontal) RG tracks of g 2 ^, j values 
above the top point would miss hitting the phase transition boundary and hence 
correspond to being in the Coulomb phase more or less energetically prohibited 
from being in the monopole condensate phase depending on how far above the 
top point, that the RG track is. This is also in violation of MPP. The only RG tra- 
jectory allowed by MPP is the one that goes through the fourth point. And our 
intersection point singles out just this RG trajectory. 

Actually, MPP is put to a very stringent test here. Had our intersection point 
picked out any other RG trajectory then the one that hits the fourth point our 
MPP would have been falsified. 

2.6 Conclusion and outlook 

In the figures we see the plot of the values of 3g 2 /7T - for the "corresponding 
" U(1) monopole coupling - versus our group characterising quantity dN- The 
extrapolation to the point where it crosses the curve of (dN , 3g 2 /7t) combinations 
for which the |3-function for g is zero (the crossing point is drawn both for the 
one loop and two loop calculation and lies of course exactly on the extrapolated 
curve but the intersecting curve (3g2 = is not drawn) it is remarkably close 
to being critical in the sense that its ordinate is very close to the critical value 
3g 2 /7t = 3- 15.11/7T = 14.46 (the two loop critical value for the case of zero 
beta function of g 2 is 15.11 corresponding to the top point of the phase transition 
curve). We have drawn on the plot the point with this ordinate calculated using 
the relation |3gi =0. 
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Fig. 2.1. Magnetic coupling (ordinate) as a function of dN (abcissa). The curve is de- 
termined to be the parabola that passes through the values of corresponding to at 
"experimental values of gauge couplings (at positive dN -values). The points at negative 
value of dN lying off the parabolic fit (solid line) correspond to maximum values of Sg 2 /n 
on phase transition curves calculated to one and two loops (see closeup in Fig. l2.2t . 




-0.6 -0.4 -0.2 



Fig. 2.2. Closeup of ^j— vs. dN for slightly negative values of dN • The two points (at 
negative d) lying on the parabolic fit (solid line) correspond to to the intersection of the 
parabolic fit with the curves |3 g 2 =0 calculated for one- and two loops the latter corre- 
sponding to the point at dN ~ —0.57. The two points (also at negative values of dN) not 
on the parabolic fit correspond to the maximum values of 3g 2 /7t on the phase transition 
curves calculated to one and two loop. The latter corresponds to the point at dN =~ —0.62. 
The points for highest order (i.e., at dN ~ —0.57 on curve and at dN ~ —0.62 off the curve) 
are seen to converge towards one another relative to the remaining two points calculated 
to one loop. 
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It is this coincidence which is either accidental or because our theory of Anti- 
GUT conbined with MPP is working so well that it only deviates in say the in- 
verse U(1 ) fine structure constant at the Planck scale by about 2-1.5 = 3 units. 

If the Planck scale had not been taken to be upi = 1 .2 • 1 1 9 GeV, as was the 
case for the figure but rather to the value obtained if we had used 87tG instead of 
G as the quantity from which to determine by dimensional arguments the Planck 
scale, the latter would have been a factor \/&n smaller. This would correspond 
to subtracting 1.61 from the In |j.pi. or 3.22 from the In [ip V With this lower value 
for the fundamental energy scale the U(1 ) running coupling would be shifted up 
in value from the 55.4 of equation (23) by 1.84 to 57.2 which in turn would shift 
the 55.4/2=27.7 value plotted on the curve up by 0.92. This is seen by eye to shift 
the curve even closer to going through the fourth point than was the case for 
the Planck scale calculated simply using G. Have in mind that the SU(2) point 
with almost compensating |3 -function contributions from the Yang Mills and the 
fermion contributions would only be moved very little by a slight change in the 
Planck scale so that this point would be almost unchanged by the replacement of 
G with 8nG, while the SU(3) point would go the opposite way meaning down 
in 3g 2 /7t-value, so that the curve would essentially be tilted so as to rise more 
slowly with increasing • 

Assuming that the correct fundamental scale should be obtained from 87tG 
rather than from G, the deviation would in terms of the inverse U( 1 ) fine structure 
constant be only about one out of 27.7 meaning about 1 in 55.4, which is only 
about 2%. We can not meaningfully expect better coincidence unless we were to 
calculate a three loop critical value since going from one to two loop in the critical 
monopole coupling squared gave a 20% correction so that 4% would be expected 
from three loops But 4% in the critical monopole coupling square would mean 
about 2% in the U(1 ) finestructure constant (inverted) at the Planck scale. 

We may also be concerned that we have not yet made two loop calculations 
of our correction factors C giving the correspondance between the nonAbelian 
couplings and the corresponding U(1 ) coupling. But it is our experience numeri- 
cally that the remarkable result of the intersection point having critical coupling 
is rather insensitive to the exact C-system used. We therefore believe that even 
with the C's only computed to one loop the accuracy of the calculation is already 
effectively a two loop calculation. 

The question of the 87tG versus the G in determining the Planck scale or 
fundamental scale must be considered basically an uncertainty at least until after 
a very detailed discussion of this point. 

Our result brings to mind the proverbial story of trying to find a needle in 
a large haystack. Actually we could even tell a better story. Start by standing at 
a distance from the haystack with an electric torch with a very narrow beam. 
Initially the torch is turned off. Now imagine that our task is to aim the torch in 
the dark it before we turn the torch on - we are not allowed to move the torch 
after turning it on. Now ask about the likelyhood of capturing the needle in our 
fixed narrow beam of light. We could improve the story by claiming that we stand 
with our torch at the Planck scale and must aim it at a haystack at the scale of the 
monopole mass before turning it on. There is even the added complication that 
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we a priori don't know exatly what we're looking for. But we miraculously find 
the needle in our narrow beam of light. But maybe with the guiding light of our 
model this is not such a miracle after all. Nature may be telling us that we're on 
the right track with our MPP/AntiGUT model. 

In this work we have only concerned ourselves with two of the three inter- 
sting phases for monopoles, namely the Coulomb phase and the phase with a 
monopole condensate. These phases have been treatable using the Higgs mono- 
pole description of monopoles as particles. We would also like to have the mono- 
pole confining phase at our multiple point. For every new phase we bring to the 
multiple point there is one more intensive parameter that becomes finetunned by 
MPP. This could be ratio of the monopole mass scale to that of the Planck mass 
if we had the confining phase. However taking the confining phase into account 
would require appending a string scenario to our approach here and doing this 
with the the high accuracy we otherwise have in this work is not yet possible. 
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Abstract. The approach unifying all the internal degrees of freedom - proposed by one of 
us| 4 2 3 4 5 6 7 8 9 10 1 - is offering a new way of understanding families of quarks and lep- 
tons. Spinors, namely, living in d ( = 1 + 1 3 ) —dimensional space, manifest in the observed 
d(= 1 + 3) -dimensional space (at "physical energies") all the known charges of quarks 
and leptons (with the mass protection property of the Standard model that only the left 
handed quarks and leptons carry the weak charge while the right handed ones are weak 
chargeless included), while a part of the starting Lagrange density in d {— 1 + 13) trans- 
forms the right handed quarks and leptons into the left handed ones, manifesting a mass 
term in d = 1 + 3. Since a spinor carries two kinds of spins and interacts accordingly with 
two kinds of the spin connection fields, the approach predicts families and the correspond- 
ing Yukawa couplings. In the paper|ll| the appearance of families of quarks and leptons 
within this approach was investigated and the explicit expressions for the corresponding 
Yukawa couplings, following from the approach after some approximations and simplifi- 
cations, presented. In this paper we continue investigations of this new way of presenting 
families of quarks and leptons by further analyzing properties of mass matrices, treating 
quarks and leptons in an equivalent way. Since it is a long way from the starting simple 
action for a Weyl spinor in d = 1 + 1 3 to the observable phenomena at low energies, and 
we yet want to make (at least rough) predictions of the approach, we connect free param- 
eters of the approach with the known experimental data and investigate a possibility that 
the fourth family of quarks and leptons appears at low enough energies to be observable 
with the new generation of accelerators. 

3.1 Introduction 

The Standard model of the electroweak and strong interactions (extended by as- 
suming nonzero masses of the neutrinos) fits with around 25 parameters and con- 
straints all the existing experimental data. However, it leaves unanswered many 
open questions, among which are also the questions about the origin of the fam- 
ilies, the Yukawa couplings of quarks and leptons and the corresponding Higgs 
mechanism. Understanding the mechanism for generating families, their masses 
and mixing matrices might be one of the most promising ways to physics beyond 
the Standard model. 
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The approach, unifying spins and charges[4 23456789 10 11], might by 
offering a new way of describing families, give an explanation about the origin 
of the Yukawa couplings. 

It was demonstrated in the references ll6l8l9H0l that a left handed SO(1 , 13) 
Weyl spinor multiplet includes, if the representation is analyzed in terms of the 
subgroups SO(1 , 3), SU(2), SU(3) and the sum of the two U(l )'s, all the spinors of 
the Standard model - that is the left handed SU(2) doublets and the right handed 
SU(2) singlets of (with the group SU(3) charged) quarks and (chargeless) lep- 
tons. There are the (two kinds of) spin connection fields and the vielbein fields in 
d = (1 + 13)— dimensional space, which might manifest - after some appropriate 
compactifications (or some other kind of making the rest of d — 4 space unobserv- 
able at low energies) - in the four dimensional space as all the gauge fields of the 
known charges, as well as the Yukawa couplings. 

The paper [11] analyzes, how do terms, which lead to masses of quarks and 
leptons, appear in the approach unifying spins and charges as a part of the spin 
connection and vielbein fields. No Higgs is needed in this approach to "dress" 
right handed spinors with the weak charge, since the terms of the starting La- 
grangean, which include y°y s , with s — 7,8, do the job of a Higgs field. 

Since we have done no analyses (yet) about the way of breaking symmetries 
of the starting group SO(1 , 13) to SO(1 , 7) x U(1 ) x SU(3) and further within 
our approach (except some very rough estimations in ref.[12]), we do not know 
how might symmetry breaking in the ordinary space influence the fields (spin 
connections and vielbeins), which in the starting action determine the Yukawa 
couplings. We also do not know how do nonperturbative and other effects (like 
boundary conditions) influence after the break of symmetries the couplings of 
spinors to the part of the starting gauge fields which in d = 1 +3 manifest at low 
energy as the Yukawa couplings. It is namely a long way from a simple starting 
action for a spinor with only one parameter and only the gravity as a gauge field 
to the observable quarks and leptons interacting with all the known gauge fields. 

We can accordingly in this investigation, by connecting Yukawa couplings 
with the experimental data, only discuss about the appearance of the "vacuum 
expectation values" of the spin connection fields which enter into the Yukawa 
couplings, trying to guess how all the complicated breaks of the symmetries have 
lead from the starting action to the observable massess and mixing angles. We 
also have no explanation yet why the second kind of the Clifford algebra objects 
do not manifest in d = 1 + 3 any charges, which could appear in addition to the 
known ones. 

Since the generators of the Lorentz transformations and the generators of 
families commute, and since only the generators of families contribute to non- 
diagonal elements of mass matrices (which means, that the off diagonal matrix 
elements of quarks and leptons are strongly correlated), the question arises, what 
makes leptons so different from quarks in the proposed approach. Can it be that 
at some energy level they are very alike and that there are some kinds of bound- 
ary conditions together with the nonperturbative effects which lead to observ- 
able properties, or might it be that Majorana like objects, not taken into account 
in these investigations up to now, are responsible for the observed differences? 
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To evaluate whether this way of going beyond the Standard model of the 
electroweak and colour interaction is in agreement with the nature, we must first 
make a rough estimate of what the approach predicts, before going to more so- 
phisticated and therefore also more trustable predictions. Starting with only the 
gauge gravity in d = 1 + 13 (or in any d) and then coming down to the "the 
physical world", is a huge project, which needs to be made in several successive 
steps. 

What turns out in our approach to be exciting is that one Weyl spinor in 
d=l+13 offers all the quarks and the leptons postulated by the Standard model 
and with just right quantum numbers answering the open question of the Stan- 
dard model, how it is at all possible that the weak charge "knows" for the hand- 
edness (left handed weak charged quarks and leptons and right handed weak 
chargeless quarks and leptons), since in the Standard model the handedness and 
the weak charge belong to totally separated degrees of freedom. 

Next exciting thing in the approach is, that it is a part of the starting La- 
grangean which does, what the Standard model requires from the Higgs: makes 
the nonzero matrix elements between the weak chargeless right handed quarks 
and leptons and weak charged quarks and leptons, correspondingly, which then 
manifest as the mass matrices. 

The third exciting thing of our approach is that it offers a mechanism for the 
appearance of families and accordingly the possibility to calculate the Yukawa 
couplings "from the first principle". 

There are also severe problems, which the approach is confronting in the way 
down to the d = 1 + 3 world. 

First severe problem is, how can at all exist at low energies any "non Planck 
scale" mass, if one starts from a very high scale? We are working hard[13J to over- 
come this "Witten's no go" theorem[14J, which concerns all the Kaluza-Klein-like 
theories. We found for a toy model (d = 1 + 5) a particular boundary condi- 
tion, which by requiring that spinors of only one handedness exist on the bound- 
ary makes that no Dirac mass can occur and therefore also the corresponding 
Yukawas do not appear after a particular break of a symmetry. We also found 
[15 1 that a Weyl spinor with no charge has no Majorana mass only in some di- 
mensions, and d=l+13 is one of those. In this paper we assume that the extension 
of the toy model can work also in our more generalized case. The justification is 
under consideration. 

One further severe problem is how to treat the "history" of spinors after the 
breaks of symmetries with all the perturbative but mostly nonperturbative ef- 
fects, which are "dressing" spinors (quarks and leptons), since even in the hadron 
physics a similar problem is not yet solved. In this paper we assume that the 
"dressing" manifests, after going beyond "the tree level" in different "vacuum 
expectation values" of the omega fields. 

Even proceeding in this way, this paper (making a first step towards more 
justified results by allowing several approximations and assumptions in order to 
find out, whether something very essential and unexpected can go wrong with 
our approach) was quite a work. We indeed come in this paper and in the previ- 
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ous one, which this one is following to promising results, that the approach might 
have a chance to go successfully beyond the Standard model. 

In this paper we try to understand properties of quarks and leptons within 
the approach unifying spins and charges treating quarks and leptons equiva- 
lently Within this approach we discuss also a possibility that the fourth family 
of quarks and leptons appears at low enough energies to be observable with new 
accelerators. 

In Sect l3.2l of this paper we present the action for a Weyl spinor in (1 + 13)- 
dimensional space and the part of the Lagrangean, which manifests at "physical 
energies" as an effective Lagrangean, with the Yukawa mass term included. This 
section is a brief repetition of the derivations presented in the ref. fTTl . 

Also Sect l3.3l is a short summary of the paper [11] in the part in which the 
explicit expression for the four mass matrices of the four families of quarks and 
leptons is presented, derived under several assumptions and simplifications from 
the starting action of the approach unifying spins and charges. In Subsect f3.3.1l 
we study properties of the mass matrices in the approximation, that the "vacuum 
expectation values" of the gauge fields of the second kind of the Clifford alge- 
bra objects entering into mass matrices are the same for all the quarks and the 
leptons. In Subsect f3.3.2l we study properties of the mass matrices relaxing this 
requirement. 

In Subsect f3.3.3l we discuss the problem of the appearance of negative masses 
in connection with the internal parity, defined within the presented approach. 

In Sect l3.4l we present the numerical results after fitting the free parameters 
of the mass matrices with the experimental data, predicting masses and mixing 
matrices of the four families of quarks and leptons. 

In Sect l6.7l we comment on the success of our approximate prediction of our 
approach. 

3.2 Weyl spinors in d = (1 +13) manifesting at "physical 
energies" families of quarks and leptons 

We assume a left handed Weyl spinor in (1 + 13)-dimensional space. A spinor 

carries only the spin (no charges) and interacts accordingly with only the gauge 

gravitational fields - with spin connections and vielbeins. We assume two kinds of 

the Clifford algebra objects and allow accordingly two kinds of gauge fields lll2l3l4l5l6l7l8l9llOIH"l . 

One kind is the ordinary gauge field (gauging the Poincare symmetry in d = 

1 +13). The corresponding spin connection field appears for spinors as a gauge 

field of S ab = j(y a y b — y b y a ], where y a are the ordinary Dirac operators. The 

contribution of these fields to the mass matrices manifests in only the diagonal 

terms - connecting the right handed weak chargeless quarks or leptons to the left 

handed weak charged partners within one family of spinors. 

The second kind of gauge fields is in our approach responsible for families of 
spinors and couplings among families of spinors - contributing to diagonal ma- 
trix elements as well - and might explain the appearance of families of quarks and 
leptons and the Yukawa couplings of the Standard model of the electroweak and 
colour interactions. The corresponding spin connection fields appear for spinors 
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as gauge fields of S ab (S ab = j(y a y b — y b y a )) with y a , which are the Clifford 
algebra objects|2 16J, like y a , but anticommute with y a . 

Following the ref.[HJ we write the action for a Weyl (massless) spinor in 
d(= 1 + 13) - dimensional space as follow^ 



S = 



d d x£, 



POa = Pa — 2^ aba,Qbct _ 2^ ab< ^ aba ' 

Here f a Q are vielbeins (inverted to the gauge field of the generators of trans- 
lations e Q a , e Q a f a b = Sg, e a a fP Q = bj>), with E = det(e a a ), while cu aba 
and Co aba are the two kinds of the spin connection fields, the gauge fields of 
S ab and S ab , respectively, corresponding to the two kinds of the Clifford algebra 
objects lli7llOI , namely y a and y a , with the properties 

{Y Q ,Y b }+=2r| ab ={y Q ,y b } + , {y",f b } + =0, (3.2) 

leading to {S ab , S cd }_ = 0. We kindly ask the reader to learn about the properties 
of these two kinds of the Clifford algebra objects - y a and y a and of the corre- 
sponding S ab and S ab - and about our technique in the ref.[HJ or the refs. 11171161 . 

One Weyl spinor representation in d = (1+13) with the spin as the only inter- 
nal degree of freedom, manifests, if analyzed in terms of the subgroups SO ( 1 , 3) x 
U(l)xSU(2)xSU(3) in four-dimensional physical space as the ordinary ( S O ( 1 , 3 ) ) 
spinor with all the known charges of one family of the left handed weak charged 
and the right handed weak chargeless quarks and leptons of the Standard model. 
The reader can see this analyses in the paper[ll [ (as well as in several references, 
like the one in the ref. lflOl ). 

We may rewrite the Lagrangean of Eq. (|3.1b so that it manifests the usual 
(1 + 3)— dimensional spinor Lagrangean part and the term manifesting as a mass 
termilll 

£ = ^y m (Pm-X g A T Ai A^)^ + 

A,i 

Y_ $Y s Vos A> + the rest. (3.3) 

s=7,8 

Index A determines the charge groups (SU(3), SU(2) and the two U(l )'s), index i 
determines the generators within one charge group. t Ai denote the generators of 



1 Latin indices a, b, .., m, n, .., s, t, .. denote a tangent space (a flat index), while Greek 
indices a, (3, .., \i,y, ..cr,T.. denote an Einstein index (a curved index). Letters from the 
beginning of both the alphabets indicate a general index (a, b, c, .. and <x, |3,y, .. ), from 
the middle of both the alphabets the observed dimensions 0, 1, 2, 3 (m, n, .. and ..), 
indices from the bottom of the alphabets indicate the compactified dimensions (s,t, .. 
and <t,t, ..). We assume the signature n ab = diag{1 , — 1 , — 1 , • • • ,— 1}. 
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the charge groups 

_Ai V - „Ax est 

s,t 

{ T Ai ,T Bi }_ =i5 AB f Aiik T Ak , (3.4) 

with s,t e 5, 6, .., 14, while A Al , m = 0,1,2,3, denote the corresponding gauge 
fields (expressible in terms of tu stm ). 

We have: Y = t 41 + t 21 , Y' = t 41 - t 2 \ with t 1 1 := }(S 58 - S 67 ), t 12 := 
1(5 57 +5 68 ),t 13 :=1(5 56 -5 78 ),t 21 :=1(5 56 +5 78 ),t 31 := 1(5 9 12 -S 10 11 ), 

T 32 . = 1(^9 11 + 5 10 T 33 . = ^9 10 _ 5 11 U ]f T 34 . = ^9 14 _ ^10 13^ 

T 35 . = f (5 9 13 + 5 10 14 }/ T 36 . = 1 14 _ ^12 13^ T 37 . = ^^11 13 + £12 14^ 

x 38 . = 1 ( 5 9 10 +5 11 12_ 2i ol3 14 ]/T 41 . = _ 1 (£9 1 + £1 1 12+^13 14) 

2v3 3 

The subgroups are chosen so that the gauge fields in the physical region 
agree with the known gauge fields. If the break of symmetries in the S ab sec- 
tor demonstrates the same symmetry after the break as in the S ab sector, then 
also the corresponding operators with f Al should be defined. 

Making several assumptions, explained in details in the ref.[ll| - we shall 
repeat them bellow - needed to manifest the observable phenomena (and can not 
yet be derived, since we do not yet know how the break of symmetries influ- 
ences the starting Lagrangean), we are able to rewrite the mass term of spinors 
(fermions) from Eq. d3.3t (Y.s=7 8 1 l'Y s Pos iK neglecting the rest) by assuming 
that they are small in comparison with what we keep at "physical energies") as 
Ly, demonstrating the Yukawa couplings of the Standard model 

£ Y =W{ (+) ( Y. v A ++ Y vH) + 

V=y,y y=N + ,NJ,-f 13 ,Y,Y' 

78 «— 

(-) ( Y y A -+ Y yA y ) + 

y=Y,Y' y = N + ,N- ,-t 13 ,Y,Y' 

7g oc bd 

(+) Y_ (k)(l) A kl ((ac),(bd)) + 

{(ac)(bd)},k,l 
7g ac bd 

(-) Y. (k)(l) A kl ((ac),(bd))}tJ;, (3.5) 

{(ac)(bd)},k,l 

with 



, aa 



ab 1 ^aa ab i c 

(k):= 2 (yQ + V y8) ' M=2 [r + ^k fb) (3 - 6) 

ab 

and with k = ±1 , if r] aa T] bb = 1 and ±i, if ri aa ri bb = -1 . While (k) are expressible 

ab 

in terms of ordinary y a and y b , (k) are expressible in terms of the second kind of 
the Clifford algebra objects, namely in terms of y a and y h . 

The Yukawa part of the starting Lagrangean (Eq. l|3.5b ) has the diagonal terms, 
that is the terms manifesting the Yukawa couplings within each family, and the 
off diagonal terms, determining the Yukawa couplings among families. 
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The operators, which contribute to the non diagonal terms in mass matri- 
ces, are superposition of S ab (times the corresponding fields tD abc ) and can be 
represented as factors of nilpotents 

ab cd 

(k)(t), (3-7) 

with indices ( ab) and (cd) which belong to the Cartan subalgebra indices and the 
superposition of the fields d> Qbc . We may write accordingly 

1 78 78 QC bd 

Y_ ~2 (±)S ab cI> ab± =- Y. (±)(k)(l) A^((ac),(bd)), (3.8) 

(a,b) (ac),(bd), k,l 

where the pair (a, b) in the first sum runs over all the indices, which do not char- 
acterize the Cartan subalgebra, with a, b = 0, . . . , 8, while the two pairs (ac) and 
(bd) denote only the Cartan subalgebra pairs (for SO(l,7) we only have the pairs 
(03), (12); (03), (56) ;(03), (78); (12), (56); (12), (78); (56), (78)) ; k and I run over 
four possible values so that k = ±i, if (ac) = (03) and k = ±1 in all other cases, 
while I = ±1. The fields A^. l ((ac), (bd)) can then be expressed by cu ab ± as fol- 
lows 

AJ + ((ab), (cd)) = ~[w ac ± - -w bc± -ia> ad± - -cl> bd± ), 



i 



r r 
i 1 



A ± ((ab), (cd)) = --(u3 ac± + -cD bc± + icD Qd± - -cD bd ±), 
2 r r 

A± + ((ab), (cd)) = -^(d> ac± + -d) bc ± -id> ad± + -u) bd± ), 
2 r r 

- ii 1 

Aj ((ab), (cd)) = --(cI3ac± - -d>bc± +i*ad± + -u>bd±), (3.9) 
2 r r 

with r = i, if (ab) = (03) and r = 1 otherwise. We simplify the index kl in the 
exponent of the fields A kl ±((ac), (bd)) to ±, omitting i. 

We must point out that a way of breaking any of the two symmetries - the 
Poincare one and the symmetry determined by the generators S ab in d = 1 + 13 
- strongly influences the Yukawa couplings of Eq. l|3.5|l , relating the parameters 
o'abc- Not necessarily any break of the Poincare symmetry influences the break 
of the other symmetry and opposite. Although we expect that it does. Accord- 
ingly the coefficients c Al ab determining the operators t Ai in Eq. (|3.41 and the co- 
efficients c Al a b determining the operators f Al in the relations 



f A 



1 = Y_ eAi ab s Qb , 



a,b 



( f Ai )t Bj } _ =i6 AB f ~Aijk f Ak (31Q) 

might even not be correlated. If correlated (through boundary conditions, for ex- 
ample) the break of symmetries might cause that off diagonal matrix elements of 
Yukawa couplings distinguish between quarks and leptons. 

We made, when deriving the mass matrices of quarks and leptons from the 
approach unifying spins and charges, several assumptions, approximations and 
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simplifications in order to be able to make at the end some rough predictions 
about the properties of the families of quarks and leptons: 

i. The break of symmetries of the group S0(1, 13) (the Poincare group in 
d = 1 + 13) into S0(l,7) x SU(3) x U(l) occurs in a way that only massless 
spinors in d = 1 + 7 with the charge SU(3) x U(l ) survive, and yet the two U(1) 
charges, following from SO (6) and SO(1 , 7), respectively, are related. (Our work 
on the compactification of a massless spinor in d = 1 + 5 into d = 1 + 3 and a 
finite disk gives us some hope that such an assumption might be justified! 13 1.) 
The requirement that the terms with S 5a u)5 a b and S 6a a>6 a b do not contribute to 
the mass term, assures that the charge Q = t 41 + S 56 is conserved at low energies. 

ii. The break of symmetries influences both, the (Poincare) symmetry de- 
scribed by S ab and the symmetries described by S>, and in a way that there 

56 56 

are no terms, which would transform (+) into [+]. This assumption was made 
that at "low energies" only four families have to be treated and can be explained 
by a break of the symmetry SO(1 , 7) into SO(1 , 5) x 11(1 ) in the S ab sector so that 
all the contributions of the type S 5a d>5 a \, and S 6a tU6ab are equal to zero. We also 
assume that the terms which include components p s , s = 5, .., 14, of the momen- 
tum p a do not contribute to the mass matrices. We keep in mind that any further 
break of symmetries strongly influences the relations among <v a bc, appearing in 
the paper [11 J as "vacuum expectation values" in mass matrices, so that predic- 
tions in Sect l3. 41 strongly depend on the way of breaking. 

iii. We make estimations on a "tree level". 

iv. We assume the mass matrices to be real and symmetric (expecting that 
complexity and nonsymmetric properties will not influence considerably masses 
and mixing matrices of quarks and leptons). 

3.3 Four families of quarks and leptons 

Taking into account the assumptions, presented in Sect l3.21 we end up with four 
families of quarks and leptons 

03 12 56 78 

I. I (+)(+) II- 

03 12 56 78 

II. [+i][+] I (+)(+) ||... 

03 12 56 78 

III. [+i](+) I ( + )[+] ||... 

03 12 56 78 

IV. | (+)[+] ||.... (3.11) 

The Yukawa couplings for these four families are for u-quarks and neutrinos pre- 
sented on Table [37T1 where a stays for u-quarks and neutrinos. 

The corresponding mass matrix for the d-quarks and the electrons is pre- 
sented on Table [3721 where (3 stays for d-quarks and electrons. 

The explicit form of the diagonal matrix elements for the above choice of 
assumptions in terms of cu a b C s's, 6 = a, (3 and A^.,u = Y and Y', cu Q b C 6 and A^ 1 
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a 


Td 


Hr 


TTTo 


TVd 

IVR 


II 




A++((03),(12)) = 

t -*- J 327ix _ '~ t *-'01 8 oc 


A++((03)J78)) = 

* 3 8 7 a + "J 7 S a 


-A++((12),(78)) = 

^ 2 7 7 a + ^ 1 87a 


2 


2 


2 


Hl 


A-;((03),_(12)) = 

327<x~' _ 01 & ot 


A„ = A^ + 

(,U>127a — U>038cl) 


A-+((12},(78)) = 


-A-+((Q3),(78)) = 

"J 3 8 7 a - * 7 8 a 


2 


2 


2 


nii. 


A--((03),(78)) = 

" 3 3S7ct + < I, 07Sa 


-A+-((12),(78)) = 

<l>277ct— * 1 87<* 


A-a 11 = A^ + 

(tU787a — tU038oc) 


A-+((03),(12)) = 

"J327a-' I, 01 Sa 


2 


2 


2 


IV L 


A--((12),(78)) = 

"i277« + "il 87ct 


-A+-((03),(78)) = 

& 3 8 7 a - <I> 7 S a 


A+-((03),(12)) 

<t>327« — "-'Ol 8a 


Aa V = A„ + 
(d>127a + 0>787a) 


2 


2 


2 



Table 3.1. The mass matrix of four families of u-quarks and neutrinos, obtained within the 
approach unifying spins and charges under the assumptions i.-iv. (see also the ref.|ll[). 
According to Eq. J3.121 and Table |3~T1 and 13.21 there are 13 free parameters, expressed in 
terms of the fields A„ and cl> aa bc, if tlWbc are the same for different a. They then accord- 
ingly determine (due to assumptions i.-iv.) all the properties of the four families of the two 
types of quarks and the two types of leptons. 



p 


Ir 


IIr 


IIIr 


IVr 


II 


A 1 , 


A++((03),J12)) = 

" 5 327|3-" J 018|3 


-Aj+((03),(78)) = 

UJ 387p _UJ 078p 


A++((12),(78)) = 

UJ 277p + di 187p 


2 


2 


2 


II l 


Ap-((03),(12)) = 

" J 327|3-" i 01 SB 


Ap 1 = Ap + 
(CUT27P + CU038p) 


-Ae + ((12),(78)) = 

a, 277p-' i 187|3 


A p + ((03),(78)) = 

* 3 8 7 p + * 7 8 p 


2 


2 


2 


III l 


-A p -((03),(78)) = 

" J 387|3-" J 078p 


A+-((12),(78)) = 

tD 277|3- UJ 187|3 


A p "=A p + 

(U>787(3 + U>038(3 ) 


A p + ((03),(12)) = 

" J 018p + " J 327p 


2 


2 


2 


IV l 


-Ap-((12),(78)) = 

t6 277|3 + "jl87p 


A+-((03),(78)) = 

"j387 p + tI >078 p 


A+-((03),(12)) 

"j£> 1 8 p +"^327 p 


Ap v Ap + 

(d>1270 + *787|3) 


2 


2 


2 



Table 3.2. The mass matrix of four families of the d-quarks and electrons, |3 stays for the 
d-quarks and the electrons. Comments are the same as on Table l3.ll 
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is as follows 

= -Al u - -A^ + d.L u , A4=-Al^+(BLv, 

A|j = — "j A +d + + <*>+d> = — A +e + &+e> 

A « = A a + U* 3-<x + 0>12-a), A" = Ap + (iu> 3 + |3 + Cl> 12+ p), 
A L H = A L + + tt>78-a), A" 1 = Ap + (icl>03+|3 + *78+p), 

A L V = A L + (&12-a + UJ78-«), A P V = A P + + ">78+p), 

(3.12) 

with a = u,"v, (3 = d, e and — d) I ± = j(iu3o3± + cut 2± + &56± + d>78± + ^ A ± )■ 
We put the index oc = u,~v and (3 = d, e to manifest that breaking of symme- 
tries, boundary conditions and nonperturbative and other effects influence the 
"vacuum expectation values" gauge fields. If we assume that uj a bc& are equal if 
they belong to different 5, the assumption that all the matrix elements are real 
relates w 1 + — jWo^g + d>, Lb 1 - — — -ju3o38 + <&> where cu is (in case that breaking 
of symmetries does not influence quarks and leptons differently) one common 
parameter. 

If the break of symmetries does not influence the quarks and the leptons in a 
different way then under the assumptions i.-iv. the off diagonal matrix elements 
of mass matrices for quarks are the same as for the corresponding leptons (the off 
diagonal matrix elements of the u-quarks and the neutrinos are the same, and the 
off diagonal matrix elements for the d-quarks and the electrons are the same) and 
since the diagonal matrix elements differ only in a constant times a unit matrix, 
the predicted mixing matrices of the quarks and the leptons would be the same. 

We must ask ourselves at this point: Can we find a way of breaking sym- 
metries - allowing some special boundary conditions and taking into account all 
the perturbative and nonperturbative effects - which would lead to so different 
properties of quarks and leptons as experimentally observed or must we take the 
Majorana like degrees of freedom into account additionally? 

In this paper, we are not yet able to answer this question. We can only make 
some estimates trying to learn from the approach unifying spins and charges 
about possible explanations for the properties of quarks and leptons. 

We proceed by relating the experimental data and the mass matrices from 
the approach. Knowing from the experimental data that the first two families of 
quarks and leptons are much lighter than the third one, while in the refs. 1181191201 
the authors, analyzing the experimental data, conclude that the experimental data 
do not forbid masses of the fourth family of quarks to be between 200 GeV and 
300 GeV, of the fourth electron to be around 1 00 GeV and of the fourth neutrino to 
be at around 50 GeV we make one more assumption, which seems quite reason- 
able also from the point of view of the measured matrix elements of the mixing 
matrix for quarks. Namely, we assume that the mass matrices of the four families of 
quarks and leptons are diagonalizable in two steps, so that the first diagonalization trans- 
forms them into block-diagonal matrices with two 2x2 sub-matrices. This assumption, 
which means, that a real and symmetric 4x4 matrix is diagonalizable by only 
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three rather than six angles, simplifies considerably further studies, making con- 
clusions very transparent. Let us comment that such a property of mass matrices 
could be a consequence of an approximate break of symmetry in the S ab sec- 
tor from SO (1,5) to SU(2) x SU(2) x U(1), which makes, for example, all the 
terms S 7a d>7 a b6 and S 8a c&8abs contributing small terms to mass matrices. The 
exact break of this type makes that the lower two families completely decouple 
from the higher two. (Similarly we have required, in order to end up with only 
four rather than eight families, that SO ( 1 , 7) breaks to SO ( 1 , 5) x U( 1 ) so that all 
S 5a d>5a±6 and S 6a Wea±& contribute nothing to mass matrices.) 

It is easy to prove that a 4 x 4 matrix is diagonalizable in two steps only if it 
has a structure[21| 

(A B \ 
\B C — A + kB J 

Since A and C are, as assumed on Table l3.ll and Table l3.2[ symmetric 2x2 matrices, 
so must then be also B. The parameter k is assumed to be an unknown number. 

The assumption (|3.13b requires: i. CU277S = 0, 0)3275 = — ftl>i876, &7&7b = 
f&3876, ^0386 = ^ 7&078S/ 6= "u, "V, d, e, and, in the case that cu a b C 5 do not de- 
pend on 5 but only on abc. ii. k u = — and k^ = — k e , where k u and kv are two 
independent parameters. (If k = in Eq. <|3.13b , the angle of rotation is 45° - then, 
if also all the 2x2 matrices would have the same structure (namely equal diago- 
nal and equal nondiagonal elements[22 23[), the corresponding mixing matrices 
for quarks and leptons would be the identity.) 



3.3.1 Evaluation of mass matrices under assumption that cD abC 6 do not 
distinguish among quarks and leptons 

Let us first assume that neither boundary conditions nor nonperturbative or other 
effects influence the fields cu Q b C & in a way that they would differ for different 6. 

We shall present in what follows some simple relations which demonstrate 
transparently properties of mass matrices. After the one step diagonalization de- 
termined by the angle of rotation 



tancp a = ±(yi +(j) 2 ±j), tancpp = ±(y 1 + (^) 2 T ^), 

k 2 

with tan(cp a - (pp) =±-, (or ± -) (3.13) 

we end up with two by diagonal matrices, with k = k u for quarks and k = ky for 
leptons, while a concerns the u-quarks and y, and (3 the d-quarks and electrons. 

The first by diagonal mass matrix of the u-quarks (a = u) and neutrinos 
(a = y) is as follows 



l(d>018 - A/ 1 + (¥) 2 tf'w), Q« + CU127 + J1 + (T^ 2 ^078 



with a u = §A Y - lA Y ' + &- -jtDo38 + j(*f - y 1 + {^f) 1 ){&07& + &> iS7 ) and 
Oy = -A Y ' + w - jWqss + j{^r - V 1 + ( J t) 2 )( cI '078 + c&387)- The mass matrix 
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for the second two families of u-quarks (a = u) and neutrinos (a = ~v) is equal to 



A* 



'1 + (^f ) 2 (tU 78 + CU387), ^((1)018 + \J 1 + i^fV t&ISj) 

T(u>o™ + Jl + { h f) 2 tuisz), a« + cD 12 7 + Jl + { ]S fY 0>387 



Accordingly we find for the first two families of d-quarks (|3 = d) and elec- 
trons ((3 = e) 



A P 



a P) -i(d>oi8 - y 1 + (t) 11,187 

i(c&018 - v'^+iW ^187). Q P + ^127 + a/ 1 + (¥) 2 ^078 



with a d = -±A Y + §A Y ' + cu + ltu 38 - i(¥ + V 1 + (^Mi&ozs - d3 387 ) and 



a e =-A Y + tu + lcu 38 -i( J x + yl + ( ) 2 ) ( cuo78 - C&387 ) ■ k« in H314) is k u for 
d-quarks and ky for electrons. For the second two families of d-quarks (|3 = d) 
and electrons (|3 = e) it follows 



A b 



a P + V 1 + ( T) 2 (^078 -C&387), -ji&OlS + \J 1 + f^) 2 a>187 
-i(tU 18 + a/ 1 + i^f) 2 ^187), a p +CU 127 - v/l +( J x) 2 CU387 



Again, k a in l|3.14t is k u for d-quarks and ley for electrons. 

There are three angles, which in the two step orthogonal transformations 
rotate each mass matrix into a diagonal one. The angles of rotations for u— quarks 
and d— quarks, and accordingly for neutrinos and electrons, are related as seen 
from Eq. ll3.13t for the first step rotation. It follows namely that tan cp a = tan -1 cpp 
and accordingly 

71 7T (p 71 (p 

with cp = cp« — (pp. (3-14) 

Similarly also the two angles of rotations of the two by two diagonal matrices 
are related. Reminding the reader that in the unitary transformations (S^S = I) 
the trace and the determinant are among the invariants, while the angle of rota- 
tion, which diagonalizes 2 by 2 matrices (of the type l|3.13b ), and the values of the 
diagonal matrices are related as follows 



,C A , ~. C A . 
tanO = U / 1+ (^) 2 T ^, 



Ai ,2 = + A) ± yf (C - A) 2 + (2B) 2 ), (3.15) 

where for A, B, C the corresponding matrix elements from Eqs. (|3. 1413. 1413. 1413. 141 
must be taken, one easily finds that 

a r[a = - a rip, h r[ a = - b r)|3, a = u,-v, (3 = d,e, (3.16) 
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where index a denotes the first two families and b the second two families of 
either quarks (a = u, |3 = d) and leptons (a = v, (3 = e) and n = C 2B A . One then 
finds the relations, equivalent to those of Eq. ll3.14l l 

ah a , b a , b _7T Q ' b Cp Qb 71 a ' b (p 

Va — -J - <PP> <P« = -7 (Pp — ^ H 2~ ' 

with Q ' b (p = tt ' b (p a - tt > b (pp. (3.17) 

We accordingly find for the case, that cu a b C 6 do not depend on 6 the following 
relations among the masses of the quarks and the leptons |m U2 — m u , I = |ra d2 — 
m dl |, |m U4 - m-ujl = |m d4 - ra dj |, |ttu„ 2 - rrw, | = |m e2 - m ei |, |rrtv 4 - rrtyj = 
|m e4 - m e J, |(m U4 + m U3 ) - (m U2 + m u , )| = |(m d4 + m d3 ) - (nx d2 + m d , )|, 
|(rrw 4 + ra-yj - (ra-v 2 + ra-v, )| = |(m e4 + vn ei ) - (nx e2 + nx ei )|, |m U4 + m U3 | « 

1 + (%»-) 2 UJ 3 87 ~ |m d4 + TTldJ, |TOv 4 + TTVy, I ~ 2 J] + (%^) 2 U! 



-j - ./ ~ i M >-a 4 i '"-c^b i M ^v 4 i '"-V3I ~ y • t I ^387 ~ 

|m e4 + TTVe 3 1. We take the absolute values of the sums and the differences, since 
whenever an eigenvalue A12 (Eq l3.15l l appears to be negative, an appropriate 
change of a phase of the corresponding state transforms the negative value into 
the positive one by changing simultaneously the internal parity of the particular 
state, as it will be discussed in Sect l3.3.3l 

The above relations among masses of quarks and leptons do not agree with 
the experimental data, as expected. We can take them as a very rough estimation 
in the limit when masses of the fourth family are much higher than the mass rat, 
knowing that rat is more than 30 times larger than the mass rab • 



3.3.2 Predictions with relaxed assumptions 

We shall make in this subsection the evaluation of the mass matrices and their 
properties by allowing that cD a b C 6 depend on the type of quarks and leptons, 
assuming that boundary conditions connected with breaking of symmetries, per- 
turbative, nonperturbative and other effects, appearing after each break of sym- 
metries influence the "vacuum expectation values" of the gauge fields, entering 
into mass matrices so that cl) a bc5 are not the same for all the quarks and the lep- 
tons. To find out, how do they differ, one should make very sophisticated calcula- 
tions, which even under the assumption that one can treat gravity in the region far 
away from the Planck scale as all the other gauge fields, is a very tough work, not 
only because we do not yet know how to treat all the breaks of symmetries to end 
up with massless spinors before the Yukawa couplings take care of their mass, 
but also because the nonperturbative effects are not solved yet even in hadron 
physics. 

Learning from Subsect f3.3.1l that diagonalization of 4 x 4 matrix in two steps 
enables a transparent view on properties of the mass matrices and recognizing 
that such an assumption is at least not in disagreement with the known experi- 
mental data, we shall make calculations under the assumption that the diagonal- 
ization in two steps is a meaningful simplification. We shall also keep the relation 
from Subsect |3.3.11 which requires that the angles of rotations for the u-quarks and 
the d-quarks mass matrices, as well as for the neutrinos and the electrons mass matrices, 
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which determine the first and the second step of diagonalization, are simply related, just 
as presented in Eqs. ll3.14l3.17b . These assumptions reduce considerably the num- 
ber of free parameters (and might also help us to recognize the way of breaking 
symmetries in more sophisticated calculations). 

It follows then that in Eqs. (|3. 1413. 1413. 1413.141 all the fields cD abc6 carry an ad- 
ditional index 6 = a, (3, while we keep the relations k a = — kp, a = u,~v and (3 = 
d, e, where k a ^ define the first step orthogonal transformations leading to 2 by 2 
by diagonal mass matrices and the relations among the angles of rotations in the 
second step of orthogonal transformations determined by a,b r| a = a ' b ( r^rr )a, 
requiring that (Eq. J3.16l )') a,b r) a = — a,b r|p (where A, B, Care to be replaced by the 
corresponding matrix elements for the first two families determined by the ma- 
trix A£ p (Eqs.< 13. 141 13.141 1) and the second two families determined by the matrix 

A b ,p (Eqs.iHllII3)). 
Then it must be 



a £a (U>018|3 - \/l + (y 1 ) 2 U>187p) = (d> 18oc ~ \/ 1 1 ! ~f 1 ' C ° 1 sr ^ ! 



b £« ("Joisp + \/l + (~^) 2 U>187|3) = (d> 18a + \/ " i ! — !' ■ a'is; A i 



l £a (d>127[3 + \/l + (y 1 ) 2 U>078|3) = (a>127<x + \/ 1 + (^) 2 U) 78a) 



b £« (d>l 2 7|3 - Jl + (y 1 ) 2 U>078|3) = (*127oc- y 1 + (y 1 ) 2 Uj 07 8a), (3.18) 

where index a and b distinguish between the two by two matrices for the first 
two and the second two families, correspondingly while a = u,~v and |3 = d, e. 
The two mixing matrices for the quarks and the leptons have the form 



1 + 




2 1 


1 + 


I—) 2 




1 + 


I— ) 2 


, 2 j 



/ c((p)c( Q cp) -c((p)s(» -s(cp)c( Q cp b ) s(cp)s( a (p b ) \ 
c((p)s( Q cp) c(cp)c( a cp) -s(cp)s( a cp b ) -s(cp)c( a cp b ) 
s(cp)c( Q cp b ) -s((p)s( Q cp b ) c(cp)c( b cp) -c(cp)s( b cp) 

Vs((p)s(> b ) s((p)c(> b ) c((p)s( b (p) c((p)c(» / 



(3.19) 



with the angles (Eq |3.14|3.17| | described by the three parameters k a , a r| a , b r\ a as 
follows 



b Q cp + b cp 



cp = cp a -(pp, a cp- a cp a - Q cpp, a cp b = ^j—^. (3.20) 

We recognize that with the mixing matrix for either quarks or leptons describ- 
able by the three parameters each k a , a r\ a , h r\ oc ; a = u,"v, we have just enough 
free parameters to make any choice for the masses of the fourth family of quarks 
and leptons. To see this we just express A, B, C in any of the two 2 by 2 matri- 
ces in terms of the corresponding diagonal values that is in terms of the masses 
Tn-ai> mpi;i = 1 ,2,3,4; a — u,~v, |3 = d, e, and the parameters k a = — kp, a r\ a = 
— a r|p, b T|ot = — b r|p; a = u,"V,. The matrix elements of A£ ( a a a , a b a , a c a ) for u- 
quarks and neutrinos are expressible with the masses m a i , m«2 of the first two 
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families of u-quarks or neutrinos and the corresponding angles of rotations as 
follows 



while the expressions for the matrix A b with matrix elements b Q a , b c«, b b a fol- 
low, if we replace m a i with m a 3 and vn a 2 with rria 4 . Equivalently we obtain the 
mass matrices for the d-quarks and the electrons by replacing a by |3 in all expres- 
sions. Eq. ip.19t below demonstrates that if once the three parameters k a , a n a , b n a 
are chosen to fit the experimental data, any four masses for the fourth family of 
quarks and leptons agree with the proposed requirements. 

The starting mass matrices M«,p - the Yukawa couplings - for quarks and 
leptons, which are 4 by 4 matrices of Table [3TT1 and Table [372l are expressible with 
the matrices A£ ,b j of Eq.( t3.19|) as follows 

/ 1 lY A a _|_ A b 1_ l A a.|3 ~ A «,|3 1 A »,p~ A a,P \ 



rr a a. i Ab ^ i ( A «.p A S.p) k «.s i 



(One easily sees that the matrix M a p is equal to a democratic matrix with all 
the elements equal to ra a4 /4, with a = u,~v, d, e, if all the angles of rotations are 
equal to n/4 (that is for k a = 0, a,b r\ a — 0), while m cti) p i = 0, i = 1,3, and that 
the mixing matrices are then the identity.) 

Once knowing the matrices M a ,p one easily finds for the parameters cu a bca,p , 
with (abc) equal to (018), (078), (127), (187), (387), entering in Table ED and Ta- 
ble H2] the expressions 

1 r m «2 ^ Trial , m a 4 - Tn a3 n 



WM7<x 



2 /TTT^F Vi + ( b Ti«) 2 ' 

1 a r)« (m a2 - m a i ) b Ti« (m a4 - m« 3 ) 



1 r Q ri a (m a2 - Tn a i ) , b T)a (ra a4 - m a 3) n 



2 \ATT^uF A+IVF ' 

1 r ra a2 - ra«i , m.a 4 - m a 3 



2.^1 +(^)2 V 71 +( Q Tla) 2 V 1 +( b Tla) 2 

uj387« = — , 1 [(m a4 + Tn a3 ) - (m a2 +m a i)], 
2^1 +(^) 2 

a « = t (Trial + Ta a2 =====—]. (3.21) 

2 V 71 +( Q T1«) 2 

3.3.3 Negative masses and parity of states 

We have mentioned in the previous section that after the diagonalization of mass 
matrices of quarks and leptons, masses of either positive or negative sign can ap- 
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pear and that by changing the phase of the basis and accordingly also the internal 
parity of states we change the sign of the mass. Let us prove this. 

First we recognize that while the starting Lagrange density for spinors (Eq l3.1|l 
commutes with the operator of handedness in d(= 1 + 13)-dimensional space 
r (i,i3) ( r n,i3) = i2 7 s 03 S 12 S 56 ---S 1314 ), it does not commute with the opera- 



tor of handedness in d(= 1 + 3) -dimensional space r (1 ' 3) (P 1 



1,3) 



-i2 2 S 03 S 12 ). 



Accordingly also the term, which manifests at "physical energies" as the mass 
term fft 



y° frt = y° 



78 



78 



v=y,y 
( X v A - + 

y =Y, Y' 



V=NJ,NJ,t' 3 ,Y,Y' 



78 
( + ) 



78 



Y_ A£((ab),(cd)) 
(ac)(bd)},lc,V 

+ (-) X (^(t) A^ l ((ab),(cd)) 

{(ac)(bd)},k,l 

(Eq |3.5l l), does not commute with the P (1 ,3) , they instead anticommute 

{r (1 ' 3) ,T°fa) + =0. 

But the rest of the "effective" Lagrangean (Eq J3.3b commutes with the operator of 
handedness in d = (1 + 3) -dimensional space: 

{Y°y m (Pm-X g^A^r' 1 ' 3 '}- =o. 

A,i 

It then follows that the Lagrange density 



Y°Y m (Pm-X g A T Ai A^)-r (1 ' 3) Y Ar (1 - 3) 



A,i 



(r (1 > 3 ^) 

(3.22) 



for the Dirac spinor r ( 1 > 3 'tb differs from the one from Eq. d3.3t in the sign of the 
mass term, while the function P' 1 ,3 'i|; differs from tJ> in the internal parity, if i|j 
is the solution for the Dirac equation. Since the internal parity is just the conven- 
tion, the negative mass changes sign if the internal parity of the spinor changes. 
The same argument was used in the ref.([24J), while the ref.([25|) uses the equiv- 
alent argument, namely, that the choice of the phase of either the right or the 
left handed spinors can always be changed and that accordingly also the signs of 
particular mass terms change. 

Let us demonstrate now on Table l3~3l how does the operator of parity V , if 
postulated as 



P=Y°Y 8 Ix,with I x x m (I x ) 



(3.23) 
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Table 3.3. The 8-plet of quarks - the members of SO ( 1 , 7) subgroup, belonging to one Weyl 
left handed (r n ' 13) = -1 = r (1,7) x r (6) ) spinor representation of SO(1 , 13). It contains 
the left handed weak charged quarks and the right handed weak chargeless quarks of a 
particular colour ((1/2, l/(2\/3))). Here r' 1,3 ' defines the handedness in (1 + 3) space, 
S ' 2 defines the ordinary spin (which can also be read directly from the basic vector, since 

S ab (k)= | ( U k)andS ab M= \ M,ifS ab belong to the Cartan subalgebra set), t 1 3 defines 
the weak charge, t 21 defines the U( 1 ) charge, t 33 and t 38 define the colour charge and t 41 
another U(1 ) charge, which together with the first one defines Y and Y'. 

transform a right handed u-quark into the left handed u-quark: "Pur = ocui_, 
where a is the proportionality factor. 

One notices that, if the operator V is applied on a state, which represents the 
right handed weak chargeless (t 1 3 = 0) u-quark of one of the three colours and 
is presented in terms of nilpotents in the first row of Table l33l it transforms this 
state into the state, which can be found in the seventh row of the same table 
and represents the left handed u-quark of the same colour and spin and it is 
weak charged. Taking into account Eq. J3.23[) and Eq.(12,16) from ref.[llj one finds 
Pur = iuL, while VV — I. By changing appropriately the phases of this two basic 
states (ur and Ui_) we can easily achieve that "Pur = ul, "Pul = ur. We should 
in addition keep in mind that V must take into account also the appearance of 
families. We shall study discrete symmetries of our approach in the "low energy 
region" in a separate paper. 

3.4 Numerical results 

In this section we connect parameters cu a bcs of the Yukawa couplings follow- 
ing from our approach unifying spins and charges (after several assumptions, 
approximations and simplifications) with the experimental data. We investigate, 
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how do the parameters of the approach reflect the known data. We also investi- 
gate a possibility of making some predictions. 



3.4.1 Experimental data for quarks and leptons 

We present in this subsection those experimental data, which are relevant for our 
study: that is the measured values for the masses of the three families of quarks 
and leptons and the measured mixing matrices. 

We take in our calculations the experimental masses for the known three 
families from the ref. 1 261271 . 

m u . /GeV = (0.0015 - 0.004, 1.15 - 1 .35, 174.3 — 178.1), 
ra di /GeV= (0.004-0.008,0.08-0.13,4.1 -4.9), 
rOv t /GeV=(1 10 12 , 1 10~ n ,5 10~ n ), 

m ei /GeV= (0.0005,0.105,1.8). (3.24) 

Predicting four families of quarks and leptons at "physical" energies, we require 
the unitarity condition for the mixing matrices for four rather than three mea- 
sured families of quarks [ 26 1 

' 0.9730 - 0.9746 0.21 74 - 0.2241 0.0030 - 0.0044 \ 
0.213-0.226 0.968-0.975 0.039-0.044 . (3.25) 
0.0- 0.08 0.0-0.11 0.07-0.9993 / 



The experimental data are for the mixing matrix for leptons known very weakly[27 
'0.79-0.88 0.47-0.61 < 0.20 

0.19 - 0.52 0.42 - 0.73 0.58 - 0.82 | . (3.26) 
^0.20 - 0.53 0.44 - 0.74 0.56 - 0.81 

We see that within the experimental accuracy both mixing matrices - for quarks 
and leptons - may be assumed to be symmetric up to a sign. We then fit with these 
two matrices the six parameters k a , a r| a , b r| a , a = u,"v. 



3.4.2 Results 

We started with the explicit expressions for the Yukawa couplings suggested by 
the approach unifying spins and charges after making several additional assump- 
tions to the starting assumptions of the approach, also some approximations and 
simplifications, in order to be able to make some approximate predictions. We 
ended up with the four families of quarks and leptons, with the mass matrices 
symmetric and real and diagonalizable in two steps with accordingly three angles 
of rotation for each mass matrix and yet the angles of rotations for the u-quarks 
are related to those of the d-quarks and so are related also the angles of the two 
kinds of leptons. Accordingly the two mixing matrices are diagonalizable with 
three angles of rotations each. 
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Since to do the rigorous calculations is a huge project and the evaluations 
in this paper are only the first rough step towards more sophisticated very de- 
manded calculations, we on this step only are parametrizing the influence of 
breaking symmetries of either the Poincare group or the group defining families 
and of nonperturbative and other effect by parametrizing the "vacuum expecta- 
tion values" of gauge fields entering into mass matrices. The assumptions, which 
we made, take care of simplifying the evaluation as much as possible while mak- 
ing the properties of mass matrices and mixing matrices as transparent as possi- 
ble. 

Let we repeat that according to Subsect f3.3.2l (in particular Eq. l|3.21l l) any 
choice for the masses of the fourth family fits the experimental data, once twice 
the three angles of the orthogonal transformations, determining the two mixing 
matrices are chosen. Assuming that some kind of symmetry (the charge Y and Y', 
for example) makes the "vacuum expectation values" of the gauge fields entering 
into the mass matrices further related, we are going to test how does the require- 
ment that the ratios of the parameters cu Q b C 5 are for a chosen set [a, b, c] as close 
to rational numbers as possible influence the properties of the fourth family of 
quarks and leptons. 

We also repeat the recognition made in Subsect l3.4,ll (Eqs.l l3.25l l3~!26ll ) that 
experimental data agree, within the experimental accuracy, for either quarks or 
leptons, that mixing matrices are symmetric and determined with only three an- 
gles. (We do not pay attention on CP non conservation.) 

We shall now connect the parameters of the approach, which are left free, 
with the experimental data and try to find out what can we learn from the corre- 
sponding results. 

We fit for quarks and leptons the three angles of Eqs. ll3.14l3.17l l with the Monte- 
Carlo method under the requirement that the ratios of the parameters cu a t, c s are for a 
chosen set [a, b, c] so close to a rational number as possible. 

We allow the masses of the fourth family as follows: The two quark masses 
must lie in the range from 200 GeV to 1 TeV, the fourth neutrino mass must be 
within the interval 50 — 1 00 GeV and of the fourth electron mass within 50 — 200 
GeV. 

Fig- EI shows the results of the Monte-Carlo simulation for the three angles 
determining the mixing matrix for quarks. There are the experimental inaccura- 
cies, which determine the allowed regions for the three angles. 

The results for the quarks are presented on Table [3~4l and 13.51 (together with 
the corresponding values for leptons). 

Fig. l3.2l shows the Monte-Carlo fit for the three angles determining the mix- 
ing matrix for leptons. There are the experimental inaccuracies which limit the 
values of the three angles. Since in the lepton case the mixing matrix for the three 
known families as well as the masses for the three neutrinos are weakly known, 
the calculations for four families bring much less information than in the quark 
case. 

The results for the leptons are presented together with the results for the 
quarks on Table 13 .41 and 13.51 
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83 . 65 

Fig. 3.1. Figure shows the Monte-Carlo fit|21 1 of the experimental mixing matrix for quarks 
(Eq l3.25t with the three angles of Eq. d3.19t . The three angles define the three parameters 
k U) a r|u and b T|u (Eqs l3.13l3.15t . We make a choice among those values for the best fit, 
which makes the ratios cDabcu/cD a bcd as close to rational numbers as possible while as- 
suring that the masses of the three known families stay within the acceptable values from 
Eq. J3.24t , with no constraints on a a and the two quark masses of the fourth family lie in 
the range 200 - 1 000 MeV. 





u 


d 




e 


k 


-0.085 


0.085 


-1.254 


1.254 




-0.229 


0.229 


1.584 


-1.584 




0.420 


-0.440 


-0.162 


0.162 



Table 3.4. The Monte-Carlo fit to the experimental data|26 27] for the three parameters k, 
a r\ and b n determining the mixing matrices for the four families of quarks and leptons are 
presented. 

In Eq. fl3.27b we present masses for the four families of quarks and leptons as 
obtained after the Monte-Carlo fit 

m Ui /GeV = (0.0034, 1.15, 176.5,285.2), 

m di /GeV = (0.0046, 0.1 1 ,4.4,224.0), 

m Vi /GeV= (1 10~ 12 ,1 lO^ 1 1 ,5 ,84.0), 

m £i /GeV= (0.0005,0.106,1.8,169.2). (3.27) 

The results of the Monte-Carlo fit show[21] that the requirement, that the 
ratios of the corresponding parameters of UJ a bcS/ for a chosen set [a, b,c], for 
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Fig. 3.2. Figure shows the Monte-Carlo fit of the experimental data for the mixing ma- 
trix for leptons(Eq J3.26P . The three angles define the three parameters k^,, and b r\-v 
(Eqs l3.13l3.15t . Again we make a choice among those values for the best fit, which make 
the ratios cD a bcu/cu a bcd as close to rational numbers as possible. 
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|o>018 


21205 


42547 


0.498 


10729 


21343 


0.503 


1*078 


49536 
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0.490 


31846 


63201 


0.504 


|CU127 


50700 


101239 


0.501 


37489 


74461 


0.503 


|o>187 


20930 


42485 


0.493 


9113 


18075 


0.505 


|o>387 


230055 


114042 


2.017 


33124 


67229 


0.493 


a a 


94174 


6237 




1149 


1142 





Table 3.5. Values for the parameters d> a bc8 (entering into the mass matrices for the 
u— quarks, the d— quarks, the neutrinos and the electrons, as suggested by the approach 
unifying spins and charges after making the additional assumptions and simplifications 
as decribed in this paper) as following after the Monte-Carlo fit, relating the parameters 
and the experimental data. 

the quarks and the leptons should be as close to the rational numbers as possi- 
ble, forces the masses of the fourth family to lie pretty low. We recognize that 
the values agree with the experimentally allowed values as evaluated by the 
refs.[18 19 20 1 . Eq. (|3.21l l, however, tells us, that it is mainly the top mass and the 
masses of the fourth family which mostly influence these ratios. But since only the 
integer 2 and the half integer (1 /2) are involved (the ratios go to one only when 
all the masses of the fourth family are equal and are high in comparison with 
the top mass), we could use this result as a guide when looking for the way of 
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breaking symmetries on the way down to d = 1 + 3 in (much) more sophisticated 
calculations. 

The Monte-Carlo fit leads to the following mixing matrix for the quarks 



The estimated mixing matrix for the four families of quarks predicts quite a strong 
couplings between the fourth and the other three families, limiting (due to the as- 
sumptions and approximations we made, which manifest in the symmetric mix- 
ing matrices) some of the matrix elements of the three families as well. 

The estimated mixing matrix for the four families of leptons predicts very 
probably far too strong couplings between the known three and the fourth family 
(although they are not in contradiction with the report in[26|). 

Let us now repeat the number of input data and the number of predictions: 

i. We take as the input data the experimental masses (the three u-quark masses, 
the three d-quark masses, the three electron masses and the three (very weakly 
known) neutrino masses), 

ii. The quark mixing matrix and the (weakly known) lepton mixing matrix. 

By taking into account relations among the mass matrix elements for the 
four types of spinors (u-quarks, d-quarks, neutrinos and electrons) as suggested 
by the approach unifying spins and charges (after some additional assumptions, 
approximations and simplifications, which all seam reasonable from the point of 
view of the experimental data and the fact that we want to obtain at least some 
rough estimations and come to at least some rough predictions for the approach 
to see whether it is not in a very severe contradiction with the experimental data) 
for the "low energy region", we were able to fit within the experimental accuracy 
and using the Monte-Carlo procedure all the known experimental data and predict 
four masses (the masses of the quarks and leptons of the fourth family) and the 
corresponding mixing matrices for quarks and for leptons. 

Let us end up this section by repeating that all the predictions must be taken 
as a rough estimate, since they follow from the approach unifying spins and 
charges after several approximations and assumptions, which we made to be able 
to come in quite a short way to simple and also (very) transparent predictions, 
used as a first step to much more sophisticated calculations. 

3.5 Discussions and conclusions 



/ 0.974 0.223 0.004 0.042 \ 
0.223 0.974 0.042 0.004 
0.004 0.042 0.921 0.387 

V 0.042 0.004 0.387 0.921 / 



(3.28) 



and for the leptons 




0.486 0.697 0.497 0.177 
0.177 0.497 0.817 0.234 
\ 0.497 0.177 0.234 0.817/ 



(3.29) 



In this paper and in the previous one[llJ we study a possibility that the approach 
of one of us [I 2 3 4 5 6 7 8 9 10 J, unifying spins and charges, might be a new right 
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way for answering those of the open questions of the Standard model of the elec- 
troweak and colour interaction, which are connected with the appearance of fam- 
ilies of fermions, of the Yukawa couplings and of the weak scale: Why do only 
the left handed spinors carry the weak charge, while the right handed are weak 
chargeless? Where do the families of the quarks and the leptons come from? What 
does determine the strenghts of the Yukawa couplings and the weak scale? 

We can conclude that within the approach unifying spins and charges the answer 
to the question, why do only the left handed spinors carry the weak charge, while the right 
handed ones are weak chargeless, does exist: The representation of one Weyl spinor 
of the group SO(l,13), analyzed with respect to the properties of the subgroups 
SO(1 , 7) x SU(3) x U(1 ) of this group and further with respect to SU(2) and the 
second U(l), manifests the left handed weak charged quarks and leptons and the 
right handed weak chargeless quarks and leptons. 

The approach answers as well the question about a possible origin of (by the Higgs 
weak charge) "dressing" of the right handed quarks and leptons in the Standard model: 
The approach proposes the Lagrange density for fermions in d(= 1 + 13)-dimen- 
sional space (a simple one) in which fermions interact with only the gravity (the 
gauge fields of the momentum - vielbeins - and the two kinds of the Clifford 
algebra objects - spin connections). It is a part of the spin connection field of the 
Poincare group, which connects the right handed weak chargeless spinors with 
the left handed weak charged ones, playing the role of the Higgs field (and the 
Yukawa couplings within a family) of the Standard model. 

The approach is answering also the question about the origin of the families of quarks 
and leptons: Two kinds of the Clifford algebra objects gauging two kinds of spin 
connection fields, are assumed. One kind takes care of the spin and the charges 
and of connecting right handed weak chargeless fermions with left handed weak 
charged fermions. The other kind takes care of the families of fermions and con- 
sequently of the Yukawa couplings among the families contributing also to the 
diagonal elements. 

In the previous paper lfTTl we derived from the approach the expressions for 
the Yukawa couplings for four families of quarks and leptons. 

It is a long way from the starting simple Lagrange density for spinors carry- 
ing only the spins and interacting with only the vielbeins and spin connections 
to the observable quarks and leptons. To treat breaking of symmetries properly, 
taking into account all perturbative and nonperturbative effects, boundary con- 
nditions and other effects (by treating gauge gravitational fields in the same way 
as ordinary gauge fields, since the scale of breaking SO(l , 13) is supposed to be 
far from the Planck scale) is a huge project. 

The purpose of this paper is to estimate whether has the approach unifying 
spins and charges at all a chance to be the right way beyond the Standard model. 
Accordingly we tried to estimate in a rough way what does the approach predict 
for a low energy physics in d = 1 + 3. To be able to make any predictions (in a 
simple enough way) we made several approximations, assumptions and simpli- 
fications, which look acceptble from the point of view of the known experimental 
data and the fact that only a very preliminary prediction is looked for. 
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Approximations, assumptions and simplifications we made enable simple 
and transparent view on the mass matrices for the four families of quarks and 
leptons in terms of the spin connection fields of the two kinds and allow to pre- 
dict masses of the fourth family of quarks and leptons and of the corresponding 
matrix elements of the mixing matrices. We treat quarks and leptons equivalently, 
no Majorana leptons are taken into account in this study. 

The assumptions, approximations and simplifications we made (which are 
not the starting assumtion of our approach unifying spins and charges) are pre- 
sented in the ref.[llj and in Sect l3.2l 

i. The break of symmetries of the group S0(1, 13) (the Poincare group in 
d = 1+13) into S0(l,7) x SU(3) x U(l) occurs in a way that only massless 
spinors in d = 1 + 7 with the charge SU(3) x U(l) survive. (Our work on the 
compactification of a massless spinor in d = 1+5 into d = 1 + 3 and a finite 
disk gives us some hope that such an assumption might be justified[13|.) The 
requirement that the terms with S 5a a>5 a b and S 6a cu6ab do not contribute to the 
mass term, assures that the charge Q = t 41 + S 56 is conserved at low energies. 

ii. The break of symmetries influences both, the (Poincare) symmetry de- 
scribed by S ab and the symmetries described by S ab , and in a way that there 
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are no terms, which would transform (+) into [+]. This assumption can be ex- 
plained by a break of the symmetry S0(1,7) into S0(1 + 5) x U(1] in the S ab 
sector. We also assume that the terms which include components p s ,s=5,..,14, 
of the momentum p a do not contribute to the mass matrices. 

iii. We make estimates on a "tree level", taking effects bellow the tree level 
into account by allowing the matrix elements to depend on the type of f ermions in 
a way, that the corresponding ratios are (very close to) rational numbers - which 
seems to be acceptable by the approach. 

iv. We assume the mass matrices to be real and symmetric and diagonaliz- 
able in two steps, first in two by two by diagonal matrices and then further to 
diagonal ones (which is suggested by the fact that this would happen if SO(l , 5) 
breaks (approximately) to SU(2) x SU(2) x 11(1 ) in the S sector, while the break of 
SO(1 , 5) to SU(3) xU(l] instead, would make the fourth family (approximately) 
decoupled from the first three and would accordingly strongly change our - very 
preliminary - results. (While such a break seems to be even acceptable when de- 
scribing properties of leptons, it would predict for quarks much too strong cou- 
plings between the third and the first two families with respect to the measured 
values.) 

Taking all the above assumptions into account we then relate the free pa- 
rameters of the mass matrices with the measured experimental data within the 
experimental accuracy treating quarks and leptons equivalently (and not taking 
into account a possible existence of Majorana neutrinos). Making a rough predic- 
tion of the properties of the fourth family of the two kinds of quarks and the two 
kinds of leptons: of their masses and the corresponding matrix elements of the 
mixing matrices, we found that our results are in agreement with the analyses of 
refs. 1181191201 . We predict the fourth family masses m U4 = 285 GeV, m d4 = 224 
GeV, rrvv 4 = 65 GeV, m e4 = 129 GeV. Predictions for the couplings between the 



3 On the Origin of Families of Fermions and Their Mass Matrices 49 

fourth and the other three families seem reasonable for quarks, while for lep- 
tons the corresponding mixing matrix elements might suggest that either differ- 
ent break of symmetries in the S ab sector from the assumed one, or the Majorana 
neutrinos, or both effects should at least be further studied. 

To try to answer within the approach unifying spins and charges the open 
question of the Standard model: Why the weak scale appears as it does? a more 
detailed study of the breaks of symmetries in both sectors (S ab and S ab ) is needed. 

What we can conclude, after making in this paper a first rough step towards 
more justified results (by allowing several approximations and assumptions) in 
order to find out, whether something very essential and unexpected can go wrong 
with our approach is, that the approach of one of us unifying spins and charges 
might have a real chance to go successfully beyond the Standard model. 

Acknowledgments 

We would like to express many thanks to ARRS for the grant. It is a pleasure 
to thank all the participants of the workshops entitled "What comes beyond the 
Standard model", taking place at Bled annually in July, starting at 1998, for fruit- 
ful discussions, in particular to H.B. Nielsen. 

References 

1. N. S. Mankoc Borstnik/'Spin connection as a superpartner of a vielbein" Phys. Lett. B 
292, 25-29 (1992). 

2. N. S. Mankoc Borstnik, "Spinor and vector representations in four dimensional 
Grassmann space", /. Math. Phys. 34, 3731-3745 (1993). 

3. N. Mankoc Borstnik, "Poincare algebra in ordinary and Grassmann space and super- 
symmetry",/. Math. Phys. 36, 1593-1601(1994). 

4. N. S. Mankoc Borstnik, "Unification of spins and charegs in Grassmann space?", Mod- 
ern Phys. Lett. A 10, 587-595 (1995). 

5. N. S. Mankoc Borstnik and S. Fajfer, "Spins and charges, the algebra and subalgebra 
of the group SO(l,14) and Grassmann space, N. Cimento 112B, 1637-1665(1997). 

6. A. Borstnik, N. S. Mankoc Borstnik, "Are Spins and Charges Unified? How Can One 
Otherwise Understand Connection Between Handedness (Spin) and Weak Charge?", 
Proceedings to the International Workshop on "What Comes Beyond the Standard Model, 
Bled, Slovenia, 29 June-9 July 1998, Ed. by N. Mankoc Borstnik, H. B. Nielsen, C. Frog- 
gatt, DMFA Zaloznistvo 1999, p.52-57, |hep-ph/9905357] 

7. N. S. Mankoc Borstnik, H. B. Nielsen, "Dirac-Kahler approach connected to quantum 
mechanics in Grassmann space", Phys. Rev. 62 (04010-14) (2000). 

8. N. S. Mankoc Borstnik, "Unification of spins and charges", Int. }. Theor. Phys. 40, 315- 
337 (2001) and references therein. 

9. A. Borstnik, N. S. Mankoc Borstnik, "Weyl spinor of SO(l,13), families of spinors of 
the Standard model and their masses", Proceedings to the International Workshop on 
"What Comes Beyond the Standard Model", Bled 2000,2001,2002 Volume 2, Ed. by N. 
Mankoc Borstnik, H. B. Nielsen, C. Froggatt, D. Lukman, DMFA Zaloznistvo 2002, 
p.27-51, hep-ph/0301029, and the paper in preparation. 



50 



M. Breskvar, D. Lukman and N.S. Mankoc Borstnik 



10. A. Borstnik, N. S. Mankoc Borstnik, "The approach unifying spins and charges 
in and its predictions", Proceedings to the Euroconference on Symmetries Beyond the 
Standard Model, Portoroz, July 12-17, 2003, Ed. by N. Mankoc Borstnik, H. B. 
Nielsen, C. Froggatt, D. Lukman, DMFA Zaloznistvo 2003, p.27-51, |hep-ph/0401043"} 
|hep-ph/040 1055 hep- ph/0301029"1 

11. A. Borstnik Bracic, N. S. Mankoc Borstnik, "On the origin of families of fermions and 
their mass matrices", hep-ph/ 0512062 to appear in Phys. Rev. D in September. 

12. N. S. Mankoc Borstnik, H.B. Nielsen, "Coupling constant unification in spin charge 
unifying model agreeing with proton decay measurements", Proceedings to the 6 th 
Workshop "What Comes Beyond the Standard Models", Bled, 2000,2001,2002, Ed. by N. 
Mankoc Borstnik, H. B. Nielsen, C. Froggatt, D. Lukman, DMFA Zaloznistvo 2002, 
p.94-105,[hep-ph/0301029| 



13. N. S. Mankoc Borstnik, H. B. Nielsen, Phys. Lett. B 633771 (2006), hep-th/0311037, 
|hep-th/0509i0"T] 

14. E. Witten, "Search for realistic Kaluza-Klein theory", Nucl. Phys. B 186 412 (1981); 
"Fermion quantum numbers in Kaluza-Klein theories", Princeton Technical Rep. 
PRINT -83-1056, October 1983. 

15. N. S. Mankoc Borstnik, H. B. Nielsen, "Fermions with no fundamental charges call 
for extra dimensions", sent to Phys. Lett. B., hep-th/0608006] 

16. N. S. Mankoc Borstnik, H. B. Nielsen, "How to generate families of spinors", /. of 
Math. Phys. 44 (2003) 4817-4827, |h~ep-th/0303224] 

17. N. S. Mankoc Borstnik, H. B. Nielsen, "How to generate spinor representations in 
any dimension in terms of projection operators, J. of Math. Phys. 43, 5782-5803, 
|hep-th/0111257| 

18. V.A. Novikov, L.B. Okun,A.N. Royanov, M.I. Vysotsky "Extra generations 
and discrepancies of electroweak precision data", Phys. Lett. B 529(2002) 111, 
|hep-ph/01lT028"j 

19. M.Maltoni, V.A. Novikov, L.B. Okun, A.N. Rozanov, M.I. Vysotsky, "Extra quark- 
lepton generations and precision measurements", Phys. Lett. B476 (2000) 107, 
|hep-ph/991 1535^ 1. 

20. S.S. Bulanov, V.A. Novikov, L.B. Okun, A.N. Rozanov, M.I. Vysotsky, "Z-lineshape 
versus 4-th generation mases", Phys. Atom.Nucl. 66 (2003) 2169; Yad.Fiz. 66 (2003) 
2219, |hep-ph/0301268 v2. 

21. M. Breskvar, Diploma work, University of Ljubljana, 2004. 

22. A. Kleppe, D. Lukman, N.S. Mankoc Borstnik, "About families of quarks and lep- 
tons", Proceedings to the 6 th Workshop "What Comes Beyond the Standard Models" , Bled, 
2000,2001,2002, Ed. by N. Mankoc Borstnik, H. B. Nielsen, C. Froggatt, D. Lukman, 
DMFA Zaloznistvo 2002, p.94-105, |hep-ph/0301029j 

23. M.Breskvar, J. Mravlje, N.S. Mankoc Borstnik, "Prediction for Four Generations of 
Quarks and Leptons Suggested by the Approach Unifying Spins and Charges", Pro- 
ceedings to the 7 th Workshop "What Comes Beyond the Standard Models" , Bled, July 19- 
31, 2004, Ed. by N. Mankoc Borstnik, H. B. Nielsen, C. Froggatt, D. Lukman, DMFA 
Zaloznistvo 2004, p.l-16, [hip : ph/0412208] 

24. H. Fritsch, Phys. Lett.73 B,317 (1978). 

25. Ta Pei Cheng, Ling Fong Li, "Gauge theory of elementary particle physics", Claren- 
don Press, Oxford 1984, p.358. 

26. S.Eidelman.et al., Phys. Lett. B 592(2004)1. 

27. M.C. Gonzalez Garcia, "Global analysis of neutrino data", Phys.Scripta T121 (2005) 
72,|hepjph/0410030vl. 



Bled Workshops 
in Physics 
VOL. 7, NO. 2 



Proceedings to the 9th Workshop 
What Comes Beyond ... (p. 51) 
Bled, Slovenia, September 16-26, 
2006 



4 Cosmoparticle Physics: Cross-disciplinary Study 
of Physics Beyond the Standard Model 

M.Yu. Khlopov 

Center for Cosmoparticle Physics "Cosmion" 
Miusskaya PI. 4 
125047 Moscow 
Russia 

e-mail: Maxim . KhlopovSromal . inf n . it 

Abstract. Cosmoparticle physics is the natural result of development of mutual relation- 
ship between cosmology and particle physics. Its prospects offer the way to study physics 
beyond the Standard Model and the true history of the Universe, based on it, in the proper 
combination of their indirect physical, astrophysical and cosmological signatures. We may 
be near the first positive results in this direction. The basic ideas of cosmoparticle physics 
are briefly reviewed. 

4.1 Cosmoparticle physics as the solution of Uhroboros puzzle 

Cosmoparticle physics originates from the well established relationship between 
microscopic and macroscopic descriptions in theoretical physics. Remind the links 
between statistical physics and thermodynamics, or between electrodynamics 
and theory of electron. To the end of the XX Century the new level of this re- 
lationship was realized. It followed both from the cosmological necessity to go 
beyond the world of known elementary particles in the physical grounds for in- 
flationary cosmology with baryosynthesis and dark matter as well as from the 
necessity for particle theory to use cosmological tests as the important and in 
many cases unique way to probe its predictions. 

The convergence of the frontiers of our knowledge in micro- and macro 
worlds leads to the wrong circle of problems, illustrated by the mystical Uhroboros 
(self-eating-snake). The Uhroboros puzzle may be formulated as follows: The the- 
ory of the Universe is based on the predictions of particle theory, that need cosmology for 
their test. Cosmoparticle physics (TJ, (2), [3| offers the way out of this wrong cir- 
cle. It studies the fundamental basis and mutual relationship between micro-and 
macro-worlds in the proper combination of physical, astrophysical and cosmo- 
logical signatures. 

4.2 Cosmological pattern of particle physics 

Let's specify in more details the set of links between fundamental particle prop- 
erties and their cosmological effects. 
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The role of particle content in the Einstein equations is reduced to its con- 
tribution into energy-momentum tensor. So, the set of relativistic species, domi- 
nating in the Universe, realizes the equation of state p = e/3 and the relativistic 
stage of expansion. The difference between relativistic bosons and fermions or 
various bosonic (or fermionic) species is accounted by the statistic weight of re- 
spective degree of freedom. The very treatment of different species of particles 
as equivalent degrees of freedom physically assumes strict symmetry between 
them. 

Such strict symmetry is not realized in Nature. There is no exact symmetry 
between bosons and fermions (e.g. supersymmetry). There is no exact symmetry 
between various quarks and leptons. The symmetry breaking implies the differ- 
ence in particle masses. The particle mass pattern reflects the hierarchy of sym- 
metry breaking. 

Noether's theorem relates the exact symmetry to conservation of respective 
charge. The lightest particle, bearing the strictly conserved charge, is absolutely 
stable. So, electron is absolutely stable, what reflects the conservation of electric 
charge. In the same manner the stability of proton is conditioned by the conser- 
vation of baryon charge. The stability of ordinary matter is thus protected by the 
conservation of electric and baryon charges, and its properties reflect the funda- 
mental physical scales of electroweak and strong interactions. Indeed, the mass 
of electron is related to the scale of the electroweak symmetry breaking, whereas 
the mass of proton reflects the scale of QCD confinement. 

Extensions of the standard model imply new symmetries and new parti- 
cle states. The respective symmetry breaking induces new fundamental physical 
scales in particle theory. If the symmetry is strict, its existence implies new con- 
served charge. The lightest particle, bearing this charge, is stable. The set of new 
fundamental particles, corresponding to the new strict symmetry, is then reflected 
in the existence of new stable particles, which should be present in the Universe 
and taken into account in the total energy-momentum tensor. 

Most of the known particles are unstable. For a particle with the mass m 
the particle physics time scale is t ~ 1 /m, so in particle world we refer to par- 
ticles with lifetime t > 1/m as to metastable. To be of cosmological signifi- 
cance metastable particle should survive after the temperature of the Universe 
T fell down below T ~ m, what means that the particle lifetime should exceed 
t ~ (mpi/m) • (1/m). Such a long lifetime should find reason in the existence of 
an (approximate) symmetry. From this viewpoint, cosmology is sensitive to the 
most fundamental properties of microworld, to the conservation laws reflecting 
strict or nearly strict symmetries of particle theory. 

However, the mechanism of particle symmetry breaking can also have the 
cosmological impact. Heating of condensed matter leads to restoration of its sym- 
metry. When the heated matter cools down, phase transition to the phase of bro- 
ken symmetry takes place. In the course of the phase transitions, corresponding 
to given type of symmetry breaking, topological defects can form. One can di- 
rectly observe formation of such defects in liquid crystals or in superfluid He. In 
the same manner the mechanism of spontaneous breaking of particle symmetry 
implies restoration of the underlying symmetry. When temperature decreases in 
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the course of cosmological expansion, transitions to the phase of broken symme- 
try can lead, depending on the symmetry breaking pattern, to formation of topo- 
logical defects in very early Universe. The defects can represent the new form of 
stable particles (as it is in the case of magnetic monopoles), or the form of ex- 
tended structures, such as cosmic strings or cosmic walls. 

In the old Big bang scenario the cosmological expansion and its initial condi- 
tions was given a priori. In the modern cosmology the expansion of the Universe 
and its initial conditions is related to the process of inflation. The global properties 
of the Universe as well as the origin of its large scale structure are the result of this 
process. The matter content of the modern Universe is also originated from the 
physical processes: the baryon density is the result of baryosynthesis and the non- 
baryonic dark matter represents the relic species of physics of the hidden sector 
of particle theory. Physics, underlying inflation, baryosynthesis and dark matter, 
is referred to the extensions of the standard model, and the variety of such exten- 
sions makes the whole picture in general ambiguous. However, in the framework 
of each particular physical realization of inflationary model with baryosynthesis 
and dark matter the corresponding model dependent cosmological scenario can 
be specified in all the details. In such scenario the main stages of cosmological 
evolution, the structure and the physical content of the Universe reflect the struc- 
ture of the underlying physical model. The latter should include with necessity 
the standard model, describing the properties of baryonic matter, and its exten- 
sions, responsible for inflation, baryosynthesis and dark matter. In no case the 
cosmological impact of such extensions is reduced to reproduction of these three 
phenomena only. The nontrivial path of cosmological evolution, specific for each 
particular realization of inflational model with baryosynthesis and nonbaryonic 
dark matter, always contains some additional model dependent cosmologically 
viable predictions, which can be confronted with astrophysical data. The part 
of cosmoparticle physics, called cosmoarcheology, offers the set of methods and 
tools probing such predictions. 

4.3 Cosmoarcheology of new physics 

Cosmoarcheology considers the results of observational cosmology as the sample 
of the experimental data on the possible existence and features of hypothetical 
phenomena predicted by particle theory. To undertake the Gedanken Experiment 
with these phenomena some theoretical framework to treat their origin and evo- 
lution in the Universe should be assumed. As it was pointed out in [4J the choice 
of such framework is a nontrivial problem in the modern cosmology. 

Indeed, in the old Big bang scenario any new phenomenon, predicted by 
particle theory was considered in the course of the thermal history of the Uni- 
verse, starting from Planck times. The problem is that the bedrock of the modern 
cosmology, namely, inflation, baryosynthesis and dark matter, is also based on 
experimentally unproven part of particle theory, so that the test for possible ef- 
fects of new physics is accomplished by the necessity to choose the physical basis 
for such test. There are two possible solutions for this problem: a) a crude model 
independent comparison of the predicted effect with the observational data and 
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b) the model dependent treatment of considered effect, provided that the model, 
predicting it, contains physical mechanism of inflation, baryosynthesis and dark 
matter. 

The basis for the approach (a) is that whatever happened in the early Uni- 
verse its results should not contradict the observed properties of the modern Uni- 
verse. The set of observational data (especially, the light element abundance and 
spectrum of black body radiation) put severe constraint on the deviation from 
thermal evolution after 1 s of expansion, what strengthens the model indepen- 
dent conjectures of approach (a). 

One can specify the new phenomena by their net contribution into the cos- 
mological density and by forms of their possible influence on parameters of mat- 
ter and radiation. In the first aspect we can consider strong and weak phenomena. 
Strong phenomena can put dominant contribution into the density of the Uni- 
verse, thus defining the dynamics of expansion in that period, whereas the con- 
tribution of weak phenomena into the total density is always subdominant. The 
phenomena are time dependent, being characterized by their time-scale, so that 
permanent (stable) and temporary (unstable) phenomena can take place. They 
can have homogeneous and inhomogeneous distribution in space. The amplitude 
of density fluctuations 6 = 6p/p measures the level of inhomogeneity relative to 
the total density, p. The partial amplitude Si = bpi/ pt measures the level of fluctu- 
ations within a particular component with density pt, contributing into the total 
density p = Y.i Pi- The case 8\ > 1 within the considered i-th component cor- 
responds to its strong inhomogeneity. Strong inhomogeneity is compatible with 
the smallness of total density fluctuations, if the contribution of inhomogeneous 
component into the total density is small: pi <C p, so that 6 <C 1 ■ 

The phenomena can influence the properties of matter and radiation either 
indirectly, say, changing of the cosmological equation of state, or via direct inter- 
action with matter and radiation. In the first case only strong phenomena are rel- 
evant, in the second case even weak phenomena are accessible to observational 
data. The detailed analysis of sensitivity of cosmological data to various phe- 
nomena of new physics are presented in [3J. This set of astrophysical constraints 
confronts phenomena, predicted by cosmophenomenology as cosmological con- 
sequences of particle theory. 

4.4 Cosmophenomenology of new physics 

To study the imprints of new physics in astrophysical data cosmoarcheology im- 
plies the forms and means in which new physics leaves such imprints. So, the 
important tool of cosmoarcheology in linking the cosmological predictions of 
particle theory to observational data is the Cosmophenomenology of new physics. It 
studies the possible hypothetical forms of new physics, which may appear as cos- 
mological consequences of particle theory, and their properties, which can result 
in observable effects. 
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4.4.1 Primordial particles 

The simplest primordial form of new physics is the gas of new stable massive 
particles, originated from early Universe. For particles with the mass m, at high 
temperature T > m the equilibrium condition, n • ov • t > 1 is valid, if their 
annihilation cross section cr > 1 /(mrapi) is sufficiently large to establish the equi- 
librium. At T < m such particles go out of equilibrium and their relative concen- 
tration freezes out. More weakly interacting species decouple from plasma and 
radiation at T > m, when n • crv • t ~ 1, i.e. at Tdec ~~ (cnupi) 1 • The maximal tem- 
perature, which is reached in inflationary Universe, is the reheating temperature, 
T r , after inflation. So, the very weakly interacting particles with the annihilation 
cross section cr < 1 /(T r rapi), as well as very heavy particles with the mass m 3> T r 
can not be in thermal equilibrium, and the detailed mechanism of their produc- 
tion should be considered to calculate their primordial abundance. 

Decaying particles with the lifetime T, exceeding the age of the Universe, tu, 
t > tu, can be treated as stable. By definition, primordial stable particles sur- 
vive to the present time and should be present in the modern Universe. The net 
effect of their existence is given by their contribution into the total cosmologi- 
cal density. They can dominate in the total density being the dominant form of 
cosmological dark matter, or they can represent its subdominant fraction. In the 
latter case more detailed analysis of their distribution in space, of their condensa- 
tion in galaxies, of their capture by stars, Sun and Earth, as well as of the effects 
of their interaction with matter and of their annihilation provides more sensitive 
probes for their existence. In particular, hypothetical stable neutrinos of the 4th 
generation with the mass about 50 GeV are predicted to form the subdominant 
form of the modern dark matter, contributing less than 0,1 % to the total density. 
However, direct experimental search for cosmic fluxes of weakly interacting mas- 
sive particles (WIMPs) may be sensitive to the existence of such component |5|. 
It was shown in | 6I7I8| that annihilation of 4th neutrinos and their antineutrinos 
in the Galaxy can explain the galactic gamma-background, measured by EGRET 
in the range above 1 GeV, and that it can give some clue to explanation of cosmic 
positron anomaly, claimed to be found by HEAT. 4th neutrino annihilation in- 
side the Earth should lead to the flux of underground monochromatic neutrinos 
of known types, which can be traced in the analysis of the already existing and 
future data of underground neutrino detectors 151 . 

New particles with electric charge and / or strong interaction can form anoma- 
lous atoms and contain in the ordinary matter as anomalous isotopes. For exam- 
ple, if the lightest quark of 4th generation is stable, it can form stable charged 
hadrons, serving as nuclei of anomalous atoms of e.g. anomalous helium [91101 . 
Massive negatively charged particles bind with 4 He in atom-like systems, as soon 
as helium is formed in Big Bang Nucleosynthesis 1111 . Since particles Q~ with 
charge -1 form positively charged "ion" [ 4 HeQ~] + , which behave as anomalous 
hydrogen, the grave problem of anomalous hydrogen overproduction is inevitable 
for particle models predicting such particles, as it is the case for the model of 
"Sinister" Universe |T2) . Particles Q~ 2 with charge -2 form neutral O-helium 
"atom" [ 4 HeQ~ 2 ], in which the charge of a particle is shielded I13I14I . It removes 
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Coulomb barrier between a particle and nuclei, giving rise to new paths of nu- 
clear transformations, catalyzed by O-helium. 

Primordial unstable particles with the lifetime, less than the age of the Uni- 
verse, t < tu, can not survive to the present time. But, if their lifetime is suffi- 
ciently large to satisfy the condition t ^ (mpi/m) • (1 /m), their existence in early 
Universe can lead to direct or indirect traces. Cosmological flux of decay products 
contributing into the cosmic and gamma ray backgrounds represents the direct 
trace of unstable particles. If the decay products do not survive to the present 
time their interaction with matter and radiation can cause indirect trace in the 
light element abundance or in the fluctuations of thermal radiation. If the particle 
lifetime is much less than 1 s the multi-step indirect traces are possible, provided 
that particles dominate in the Universe before their decay. On the dust-like stage 
of their dominance black hole formation takes place, and the spectrum of such 
primordial black holes traces the particle properties (mass, frozen concentration, 
lifetime) [22 \ . The particle decay in the end of dust like stage influences the baryon 
asymmetry of the Universe. In any way cosmophenomenoLOGICAL chains link 
the predicted properties of even unstable new particles to the effects accessible in 
astronomical observations. Such effects may be important in the analysis of the 
observational data. 

So, the only direct evidence for the accelerated expansion of the modern 
Universe comes from the distant SN I data. The data on the cosmic microwave 
background (CMB) radiation and large scale structure (LSS) evolution (see e.g. 
0) prove in fact the existence of homogeneously distributed dark energy and 
the slowing down of LSS evolution at z < 3. Homogeneous negative pressure 
medium (A-term or quintessence) leads to relative slowing down of LSS evolu- 
tion due to acceleration of cosmological expansion. However, both homogeneous 
component of dark matter and slowing down of LSS evolution naturally follow 
from the models of Unstable Dark Matter (UDM) (see |3] for review), in which the 
structure is formed by unstable weakly interacting particles. The weakly interact- 
ing decay products are distributed homogeneously. The loss of the most part of 
dark matter after decay slows down the LSS evolution. The dominantly invisible 
decay products can contain small ionizing component (23| . Thus, UDM effects 
will deserve attention, even if the accelerated expansion is proved. 

4.4.2 Archioles 

Parameters of new stable and metastable particles are also determined by a pat- 
tern of particle symmetry breaking. This pattern is reflected in the succession of 
phase transitions in the early Universe. Phase transitions of the first order proceed 
through bubble nucleation, which can result in black hole formation (see e.g. fl5l 
and [16 J for review and references). Phase transitions of the second order can lead 
to formation of topological defects, such as walls, string or monopoles. The ob- 
servational data put severe constraints on magnetic monopole and cosmic wall 
production, as well as on parameters of cosmic strings. A succession of phase 
transitions can change the structure of cosmological defects. More complicated 
forms, such as walls-surrounded-by-strings can appear. Such structures can be 
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unstable, but their existence can lead the trace in the nonhomogeneous distribu- 
tion of dark matter and in large scale correlations in the nonhomogeneous dark 
matter structures, such as archioles [17], which arise in the result of cosmological 
evolution of pseudo-Nambu-Goldstone field. 

A wide class of particle models possesses a symmetry breaking pattern, which 
can be effectively described by pseudo-Nambu-Goldstone (PNG) field and which 
corresponds to formation of unstable topological defect structure in the early Uni- 
verse (see lTl6l for review and references). The Nambu-Goldstone nature in such 
an effective description reflects the spontaneous breaking of global U(l) symme- 
try, resulting in continuous degeneracy of vacua. The explicit symmetry breaking 
at smaller energy scale changes this continuous degeneracy by discrete vacuum 
degeneracy. 

At high temperatures such a symmetry breaking pattern implies the succes- 
sion of second order phase transitions. In the first transition, continuous degen- 
eracy of vacua leads, at scales exceeding the correlation length, to the formation 
of topological defects in the form of a string network; in the second phase transi- 
tion, continuous transitions in space between degenerated vacua form surfaces: 
domain walls surrounded by strings. This last structure is unstable, but, as was 
shown in the example of the invisible axion 11171181191 , it is reflected in the large 
scale inhomogeneity of distribution of energy density of coherent PNG (axion) 
field oscillations. This energy density is proportional to the initial value of phase, 
which acquires dynamical meaning of amplitude of axion field, when axion mass 
is switched on in the result of the second phase transition. 

The value of phase changes by 2n around string. This strong nonhomogene- 
ity of phase leads to corresponding nonhomogeneity of energy density of coher- 
ent PNG (axion) field oscillations. Usual argument (see e.g. [20J and references 
therein) is that this nonhomogeneity is essential only on scales, corresponding to 
mean distance between strings. This distance is small, being of the order of the 
scale of cosmological horizon in the period, when PNG field oscillations start. 
However, since the nonhomogeneity of phase follows the pattern of axion string 
network this argument misses large scale correlations in the distribution of en- 
ergy density of field oscillations. 

Indeed, numerical analysis of string network (see review in [21]) indicates 
that large string loops are strongly suppressed and the fraction of about 80% of 
string length, corresponding to "infinite" strings, remains virtually the same in 
all large scales. This property is the other side of the well known scale invariant 
character of string network. Therefore the correlations of energy density should 
persist on large scales, as it was revealed in [17 18 19]. 

The large scale correlations in topological defects and their imprints in pri- 
mordial inhomogeneities is the indirect effect of inflation, if phase transitions take 
place after reheating of the Universe. Inflation provides in this case the equal con- 
ditions of phase transition, taking place in causally disconnected regions. 

4.4.3 Primordial clouds of massive PBH 

If the phase transitions take place on inflational stage new forms of primordial 
large scale correlations appear. The example of global U(l) symmetry, broken 
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spontaneously in the period of inflation and successively broken explicitly after 
reheating, was recently considered in ||24| . In this model, spontaneous U(l) sym- 
metry breaking at inflational stage is induced by the vacuum expectation value 
(t|>) = f of a complex scalar field M 7 = tyexp(iQ), having also explicit symme- 
try breaking term in its potential V e b = A 4 (l — cos 0). The latter is negligible in 
the period of inflation, if f ^> A, so that there appears a valley relative to values 
of phase in the field potential in this period. Fluctuations of the phase 9 along 
this valley, being of the order of A0 ~ H/ (27tf ) (here H is the Hubble parame- 
ter at inflational stage) change in the course of inflation its initial value within 
the regions of smaller size. Owing to such fluctuations, for the fixed value of 060 
in the period of inflation with e-folding N = 60 corresponding to the part of the 
Universe within the modern cosmological horizon, strong deviations from this 
value appear at smaller scales, corresponding to later periods of inflation with 
N < 60. If < % the fluctuations can move the value of 0n to 0n > n in 
some regions of the Universe. After reheating, when the Universe cools down to 
temperature T = A the phase transition to the true vacuum states, corresponding 
to the minima of V e b takes place. For 9n < 7t the minimum of V e b is reached 
at Q V ac = 0/ whereas in the regions with On > tc the true vacuum state corre- 
sponds to 9 VQC = 2tc. For Qso < 7t in the bulk of the volume within the modern 
cosmological horizon VQC = 0. However, within this volume there appear re- 
gions with vac = 2n. These regions are surrounded by massive domain walls, 
formed at the border between the two vacua. Since regions with Q va c — 2n are 
confined, the domain walls are closed. After their size equals the horizon, closed 
walls can collapse into black holes. The minimal mass of such black hole is deter- 
mined by the condition that it's Schwarzschild radius, r g = 2GM/c 2 exceeds the 
width of the wall, I - f/A 2 , and it is given by M m j. n - f (mpi/A) 2 . The maximal 
mass is determined by the mass of the wall, corresponding to the earliest region 
On > 7t, appeared at inflational stage. This mechanism can lead to formation of 
primordial black holes of a whatever large mass (up to the mass of AGNs |25|). 
Such black holes appear in a form of primordial black hole clusters, exhibiting 
fractal distribution in space |24|. It can shed new light on the problem of galaxy 
formation. 

4.4.4 Antimatter stars in Galaxy 

Primordial strong inhomogeneities can also appear in the baryon charge dis- 
tribution. The appearance of antibaryon domains in the baryon asymmetrical 
Universe, reflecting the inhomogeneity of baryosynthesis, is the profound sig- 
nature of such strong inhomogeneity |26|. On the example of the model of spon- 
taneous baryosynthesis (see [27] for review) the possibility for existence of an- 
timatter domains, surviving to the present time in inflationary Universe with 
inhomogeneous baryosynthesis was revealed in [28J. Evolution of sufficiently 
dense antimatter domains can lead to formation of antimatter globular clusters 
t29ll . The existence of such cluster in the halo of our Galaxy should lead to the 
pollution of the galactic halo by antiprotons. Their annihilation can reproduce 
l30lt the observed galactic gamma background in the range tens-hundreds MeV. 
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This observed background puts upper limit on the total mass of antimatter stars 
in Galaxy (M < 10 5 M Q ). The prediction of antihelium component of cosmic 
rays [31J, accessible to future searches for cosmic ray antinuclei in PAMELA and 
AMS02 experiments, as well as of antimatter meteorites [ 32 1 provides the direct 
experimental test for this hypothesis. In this test planned sensitivity of AMS02 ex- 
periment will reach the lower limit for the mass of antimatter stars in Galaxy. This 
limit (M > 10 3 M Q ) follows from the condition that antimatter domain should be 
sufficiently large to survive and sufficiently dense to provide star formation. 

So the primordial strong inhomogeneities in the distribution of total, dark 
matter and baryon density in the Universe is the new important phenomenon 
of cosmological models, based on particle models with hierarchy of symmetry 
breaking. 

4.5 The encircled pyramid 

New physics follows from the necessity to extend the Standard model. White 
spots in representations of symmetry groups, considered in extensions of Stan- 
dard model, correspond to new unknown particles. Extension of gauge symme- 
try puts into consideration new gauge fields, mediating new interactions. Global 
symmetry breaking results in the existence of Goldstone boson fields. 

For a long time necessity to extend the Standard model had purely theoreti- 
cal reasons. Aesthetically, because full unification is not achieved in the Standard 
model; practically, because it contains some internal inconsistencies. It does not 
seem complete for cosmology. One has to go beyond the Standard model to ex- 
plain inflation, baryosynthesis and nonbaryonic dark matter. Recently there has 
appeared a set of experimental evidences for the existence of neutrino oscilla- 
tions, for effects of new physics in anomalous magnetic moment of muon (g — 2), 
for cosmic WIMPs and for double neutrinoless beta decay (see for recent review 
0). Whatever is the accepted status of some of these evidences, they indicate that 
the experimental searches may have already crossed the border of new physics. 

In particle physics direct experimental probes for the predictions of parti- 
cle theory are most attractive. Predictions of new charged particles, such as su- 
persymmetric particles or quarks and leptons of new generation, are accessible 
to experimental search at accelerators of new generation, if their masses are in 
lOOGeV-lTeV range. However, predictions related to higher energy scale need 
non-accelerator or indirect means for their test. 

Search for rare processes, such as proton decay, neutrino oscillations, neu- 
trinoless beta decay, precise measurements of parameters of known particles, ex- 
perimental searches for dark matter represent the widely known forms of such 
means. Cosmoparticle physics offers the nontrivial extensions of indirect and 
non-accelerator searches for new physics and its possible properties. In exper- 
imental cosmoarcheology the data is to be obtained, necessary to link the cos- 
mophenomenology of new physics with astrophysical observations (See (U). In 
experimental cosmoparticle physics the parameters, fixed from the consitency of 
cosmological models and observations, define the level, at which the new types 
of particle processes should be searched for (see [33)). 
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4.5.1 New quarks and leptons 

The theories of everything should provide the complete physical basis for cosmol- 
ogy. The problem is that the string theory l34l is now in the form of "theoretical 
theory", for which the experimental probes are widely doubted to exist. The de- 
velopment of cosmoparticle physics can remove these doubts. In its framework 
there are two directions to approach the test of theories of everything. 

One of them is related with the search for the experimentally accessible ef- 
fects of heterotic string phenomenology. The mechanism of compactification and 
symmetry breaking leads to the prediction of homotopically stable objects |35| 
and shadow matter [36], accessible to cosmoarcheological means of cosmoparti- 
cle physics. The condition to reproduce the Standard model naturally leads in the 
heterotic string phenomenology to the prediction of fourth generation of quarks 
and leptons 11371 with a stable massive 4th neutrino [6 1, what can be the subject of 
complete experimental test in the near future. The comparison between the rank 
of the unifying group Eg (r = 6) and the rank of the Standard model (r = 4) 
implies the existence of new conserved charges and new (possibly strict) gauge 
symmetries. New strict gauge U(l) symmetry (similar to U(l) symmetry of elec- 
trodynamics) is possible, if it is ascribed to the fermions of 4th generation. This 
hypothesis explains the difference between the three known types of neutrinos 
and neutrino of 4th generation. The latter possesses new gauge charge and, be- 
ing Dirac particle, can not have small Majorana mass due to sea saw mechanism. 
If the 4th neutrino is the lightest particle of the 4th quark-lepton family, strict 
conservation of the new charge makes massive 4th neutrino to be absolutely sta- 
ble. Following this hypothesis [37] quarks and leptons of 4th generation are the 
source of new long range interaction (y-electromagnetism), similar to the electro- 
magnetic interaction of ordinary charged particles. New strictly conserved local 
U(l) gauge symmetries can also arise in the alternative approach to extension of 
standard model [14J based on almost commutative geometry [38 J. 

It is interesting, that heterotic string phenomenology embeds even in its sim- 
plest realisation both supersymmetric particles and the 4th family of quarks and 
leptons, in particular, the two types of WIMP candidates: neutralinos and massive 
stable 4th neutrinos. So in the framework of this phenomenology the multicom- 
ponent analysis of WIMP effects is favorable. 

The motivation for existence of new quarks and leptons also follows from 
geometrical approach to particle unification [39] and from models of extended 
technicolor I40I41I . 

New quarks and charged leptons can be metastable and have lifetime, ex- 
ceeding the age of the Universe. It gives rise to a new form of stable matter around 
us - to composite dark matter, whose massive charged constituents are hidden in 
atom-like systems. Inevitable by-product of creation of such matter in Universe 
is the existence of O-helium "atom", in which positive charge of a particle is 
shielded by negative charge of massive component. O-helium can be a fraction 
fl4l or even dominant form of dark matter [13J and search for its charged con- 
stituents at accelerators and cosmic rays [42] acquires the significance of direct 
experimental test for this form of dark matter. 
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In the above approach some particular phenomenological features of sim- 
plest variants of string theory are studied. The other direction is to elaborate the 
extensive phenomenology of theories of everything by adding to the symmetry 
of the Standard model the (broken) symmetries, which have serious reasons to 
exist. The existence of (broken) symmetry between quark-lepton families, the ne- 
cessity in the solution of strong CP-violation problem with the use of broken 
Peccei-Quinn symmetry, as well as the practical necessity in supersymmetry to 
eliminate the quadratic divergence of Higgs boson mass in electroweak theory is 
the example of appealing additions to the symmetry of the Standard model. The 
horizontal unification and its cosmology represent the first step on this way, il- 
lustrating the approach of cosmoparticle physics to the elaboration of the proper 
phenomenology for theories of everything [43 1 . 

We can conclude that from the very beginning to the modern stage, the evo- 
lution of Universe is governed by the forms of matter, different from those we are 
built of and observe around us. From the very beginning to the present time, the 
evolution of the Universe was governed by physical laws, which we still don't 
know. Observational cosmology offers strong evidences favoring the existence of 
processes, determined by new physics, and the experimental physics approaches 
to their investigation. 

Cosmoparticle physics [1] [2J, studying the physical, astrophysical and cos- 
mological impact of new laws of Nature, explores the new forms of matter and 
their physical properties, what opens the way to use the corresponding new 
sources of energy and new means of energy transfer. It offers the great challenge 
for the new Millennium. Its solution for the Uhroboros puzzle is as paradoxi- 
cal, as "encircled pyramid" - cosmoparticle physics implies complex cross-dis- 
ciplinary studies, offering the multi-dimensional exit from the plane with wrong 
circle of problems in the joint of cosmology and particle physics. 
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Abstract. I briefly stipulate here some ideas, which were considered in the Discussion 
Section. 

5.1 Some phenomenological aspects of 4th generation in 
geometrical approach 

In the presented realization of geometrical approach [1J 4th generation by con- 
struction is linked to 3d generation. In such realization mixing and transitions 
between these generations are inevitable. It makes all the particles of 4th genera- 
tion unstable. It would be interesting to estimate their lifetime. But, in any case, 
this feature is important for accelerator search for quarks and leptons of 4th gen- 
eration. Processes of creation and decay of these particles at accelerators have 
rather distinct experimental signatures and can be clearly discriminated. On the 
other hand, being unstable, hadrons and leptons of 4th generation should not be 
present in cosmic rays. 

The predicted values of 4th generation particles are in some cases accessible 
to test even with the use of existing experimental data [2]. For instant, unstable 
neutrino of 4th generation with the mass - 80GeV should have been seen in LEP. 

In the approach [1 1 another problem is also of interest: if it is possible to have 
the lightest particles of the 4th generation stable. In this case the 4th generation 
is decoupled from three known families and their possible mass pattern can not 
be directly deduced from mixing with known particles and from their known 
properties. 

5.2 Composite dark matter from Technicolor 

Cosmological aspects (3) of technicolor models [4] were concentrated on stud- 
ies of possible WIMP-like candidates for dark matter species. However, this ap- 
proach also provides the possibility of stable techniparticles T with charge -2. 

If the model provides the possibility to generate in early Universe excess of 
these particles, corresponding to cosmological dark matter density, atomic bound 



64 M.Yu. Khlopov 



states of these particles with primordial 4 He can play the role of composite dark 
matter in the form of techno-O-helium 4 HeT . Formation and evolution of this 
composite dark matter will follow the trend of O-helium dark matter, studied in 
13 for the case of stable quarks of 4th generation. Experimental search for stable 
techniparticles T in cosmic rays and at accelerators is possible similar to the 
case of new stable leptons II6I42I . 

5.3 Mass self-adjustment for CLEP neutrinos 

In the model of CLEP states, offered in [8J, the following mechanism of self- 
adjustment of neutrino mass can be realized. Near a galaxy as the neutrinos ex- 
pand to the outside their density decreases and neutrino mass increases. Then 
as the mass increases the neutrinos have now a tendency to clump; after they 
clump, they become dense and therefore the neutrinos are not anymore in the 
CLEP state, but rather in a higher density state with a corresponding lower mass, 
which is now preventing neutrinos to be clustered. Therefore neutrino gas ex- 
pands again to reach low density, to return in CLEP state. 

If in the course of cosmological expansion neutrino mass in CLEP state be- 
comes too big some other interesting possibilities arise. For instance a nontrivial 
realisation is possible for Unstable Neutrino Cosmology and other effects of non- 
equilibrium particles from CLEP state neutrino decays if this mass becomes too 
big some other interesting possibilities arise. For instance a nontrivial realisation 
is possible for Unstable Neutrino Cosmology and other effects of non-equilibrium 
particles from CLEP state neutrino decays[9J. 
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Abstract. The genuine Kaluza-Klein-like theories — with no fields in addition to gravity — 
have difficulties with the existence of massless spinors after the compactification of some 
space dimensions 1 1 1. We proposed in ref. |2| a boundary condition which allows massless 
spinors compactified on a flat disk to be of only one handedness. Massless spinors then 
chirally couple to the corresponding background gauge gravitational field (which solves 
equations of motion for a free field linear in the Riemann curvature). In this paper we study 
the same toy model: M ' 1 + 3 ' x M (21 , looking this time for an involution which transforms 
a space of solutions of Weyl equations in d = 1 + 5 from the outside of the flat disk in 
x 5 and x 6 into its inside (or conversely). The natural boundary condition that on the wall 
an outside solution must coincide with the corresponding inside one leads to massless 
spinors of only one handedness (and accordingly mass protected), chirally coupled to the 
corresponding background gauge gravitational field. We introduce the Hermitean oper- 
ators of momenta and discuss the orthogonality of solutions, ensuring that to each mass 
only one solution of equations of motion corresponds. 

6.1 Introduction 

The major problem of the compactification procedure in all Kaluza-Klein-like the- 
ories with only gravity and no additional gauge fields is how to ensure that mass- 
less spinors be mass protected after the compactification. Namely, even if we start 
with only one Weyl spinor in some even dimensional space of d = 2 modulo 4 
dimensions (i.e. in d = 2(2n+ 1 ) , n = 0, 1 , 2, • • • ) so that there appear no Majorana 
mass if no conserved charges exist and families are allowed, as we have proven 
in ref. [3|, and accordingly with the mass protection from the very beginning, a 
compactification of m dimensions gives rise to a spinor of one handedness in d 
with both handedness in d — m and is accordingly not mass protected any longer. 

And in addition, since a spin (or the total conserved angular momentum) 
in the compactified part of space will in d — m space appear as a charge and 
will manifest both values (positive and negative ones) and since in the second 
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quantization procedure anti particles of opposite charges appear anyhow, dou- 
bling the number of massless spinors of both — positive and negative — charges 
when coming from d(= 2(2n + 1 )) -dimensional space down to d = 4 and after a 
second quantized procedure is not in agreement with what we observe. Accord- 
ingly there must be some requirements, some boundary conditions, which ensure 
in a compactification procedure that only spinors of one handedness survive, if 
Kaluza-Klein-like theories have some meaning. However, the idea of Kaluza and 
Klein of having only gravity as a gauge field seems too beautiful not to have the 
realization in Nature. 

One of us [4 5 6 7 8] has for long tried to unify the spin and all the charges to 
only the spin, so that spinors would in d > 4 carry nothing but a spin and interact 
accordingly with only the gauge fields of the Poincare group, that is with viel- 
beins f a a U and spin connections cu Q b a , which are the gauge fields of the Poincare 
group. 

In this paper we take (as we did in the ref . [2]) the covariant momentum of a 
spinor, when applied on a spinor function tJ>, to be 

VOa = ^aVOot, P0a4> = Pa - ^S cd O) cda . (6.1) 

A kind of a total covariant derivative of e a a (a vector with both — Einstein and 
flat index) will be taken to be p 0a e a p = ie a (3 ;ct = i(e a p, a + cu a da e d p - r" |/ p a e a Y ), 
with the require that this derivative of a vielbein is zero: e a p ;(X = 0. 

The corresponding Lagrange density C for a Weyl spinor has the form C = 
El [(t|> WpoaM>) + Wpoa^) f ] and leads to 

£ = E^ Y V( Pa -ls cd cu cdQ )xK (6-2) 

withE = det(e a a )H 

The authors of this paper have tried to find a way out of this "Witten's no 
go theorem" for a toy model of M ,1+3) x a flat finite disk in (1 + 5) -dimensional 
space [2J by postulating a particular boundary condition, which allows a spinor 
to carry after the compactification only one handedness. Massless spinors then 
chirally couple to the corresponding background gauge gravitational field, which 

1 f a a are inverted vielbeins to e u K with the properties e a a f a b = b%, e a a f |3 a = &«■ 
Latin indices a, b, .., m,n, .., s, t, .. denote a tangent space (a flat index), while Greek 
indices a, 6, .., \i,y, ..u,t:.. denote an Einstein index (a curved index). Letters from the 
beginning of both the alphabets indicate a general index (a, b, c, .. and a, |3 , y, .. ), from 
the middle of both the alphabets the observed dimensions 0, 1,2,3 (m,n, .. and u.,~v, ..), 
indices from the bottom of the alphabets indicate the compactified dimensions (s,t, .. 
and cr,T, ..). We assume the signature r| ab = diag{1 ,— 1,— 1, • • • , — 1}. 

2 To generate more than one family, we actually observe up to now three families, a sec- 
ond kind of the Clifford algebra objects has also been introduced|5 6 9 10 [, which anti 
commute with the ordinary Dirac y a matrices ({y a , y h }+ = 0) and generate equivalent 
representations with respect to the generators S ab = ^-(y a y b — y b y a ) and are used 
accordingly to generate families|5 6 9 10[. In this work we shall not take families into 
account. 
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solves equations of motion for a free field, linear in the Riemann curvature, while 
the current through the wall is for a massless and massive solutions equal to zero. 
In the ref . [2 J the boundary condition was written in a covariant way as 

&l>Lan = 0, 

withn (p) = (0,0, 0,0, cos 4>, sin cf)), n ( *> = (0, 0, 0, 0, - sin 4>, cos 4>), which are 
the two unit vectors perpendicular and tangential to the boundary of the disk (at 
Po), respectively. The projector H can for the above choice of the two vectors n' p ' 
and n' * 1 be written as 

56 1 

£=[-] = ^(l-iyV). (6-3) 

ab ab 

The reader can find more about the Clifford algebra objects (±), [±] in the Ap- 
pendix (section [6l8|l . 

The boundary condition requires that only massless states (determined by 
Eq. d6.2b ') of one (let us say right) handedness with respect to the compactified 
disk degrees of freedom are allowed. Accordingly also massless states of only 
one handedness are allowed also in d = 1 + 3. 

In this paper we reformulate the above boundary condition as an involution, 
which transforms solutions of equations of motion from outside the boundary of 
the disk into its inside. We do this by the intention that the limitation of M2 on a 
finite disk would have a natural explanation, originated in a symmetry relation. 
We also define the Hermitean momentum p s and comment on the orthogonality 
relations of solutions of equations of motion, which fulfill the boundary condi- 
tions. 

We make use of the technique presented in ref. [9 10 1 when writing the equa- 
tions of motion and their solutions. It turns out that all the derivations and discus- 
sions appear to be very transparent when using this technique. We briefly repeat 
this technique in Appendix l6.8l 

6.2 Equations of motion and solutions 

We assume that a two dimensional space, spanned by x 5 and x 6 , is an Euclidean 
manifold M (2) (with no gravity) 

f . = 6°,, to 56 . = 0. (6.4) 

and accordingly with the rotational symmetry around an origin. 

Wave functions describing spinors in (1 + 5) -dimensional space demonstrat- 
ing M' 1 + 3 ' xM' 2 ' symmetry are required to obey the equations of motion 

Y°Y a PaV 6) =0, a = m, s; m = 0, 1 , 2, 3; s = 5, 6. (6.5) 

The most general solution for a free particle in d = 1 + 5 should be written as a 
superposition of all four (2 6/2_1 ) states of one Weyl representation. We ask the 
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reader to see Appendix 16.81 for the technical details how to write one Weyl rep- 
resentation in terms of the Clifford algebra objects after making a choice of the 
Cartan sub algebra set, for which we make a choice: S 03 , S 12 , S 56 . In our tech- 
nique [ 9 [ the four states, which all are the eigenstates of the Cartan sub algebra 

ab 

set, are expressed with the following four products of projections ( [k] ) and nilpo- 

ab 

tents ((k)): 

56 03 12 
Cpl = ( + )(+!)( + ) T|>0, 

56 03 12 
Cp] = ( + )[-!][-] t|> 0> 

56 03 12 

cp? = HHJX+Ho, 

56 03 12 

(pf = [-](+!)[-] ^o, (6.6) 

where 4>o is a vacuum state. If we write the operators of handedness in d = 1 + 5 
asr d+5) = y0 Y i y 2 Y 3 y 5 Y 6( =2 3iS 03 S 12 S 56 ),ind= 1+3asr (1+3) = -iyVyV 
(= 22iS 03 S 12 ) and in the two dimensional space as P' 2 ' = iy 5 y 6 (= 2S 56 ), we find 
that all four states are left handed with respect to r ,1+5 ', with the value —1, the 
first two are right handed and the second two left handed with respect to P' 2 ', 
with the values 1 and —1, respectively, while the first two are left handed and 
the second two right handed with respect to r' 1+3 ' with the values —1 and 1, 
respectively. 

Taking into account Eq. l|6.6l l we may write a wave function tJ;' 6 ' in d = 1 +5 

as 

56 56 

^) = { A (+)+B[-]W 4 \ (6.7) 

where A and B depend on x 5 and x 6 , while ij>' 4 ' determines the spin and the 
coordinate dependent part of the wave function ij>' 6 ' in d = 1 +3. 

Spinors, which manifest masslessness in d = 1 +3, must obey the equation 

y°y s PsV 6) =0, s = 5,6, (6.8) 
since what will demonstrate as an effective action in d = 1 + 3 is 

n^ m Troi23(fdx 5 dx 6 Tr 5 6(^ (6)t y°(y m Pm+y s PsH (6) )) = 

m 

ndx m Troi2 3 (V 4)t y°Y m PmV 4) ) - [n dxmTr °'23(V 4)t YW 4) ], (6.9) 

m m 

where integrals go over all the space on which the solutions are defined, ij;' 6 ' and 
i|;' 4 ' are the solutions in d = 1 + 5 and d = 1 + 3, respectively. Troi23 and Trs6 
mean the trace over the spin degrees of freedom in x0,x1,x2,x3 and in x 5 ,x 6 , 

56 

respectively. (One finds, for example, that Tr([±]) = 1.) For massless spinors it 
must be that J dx 5 dx 6 Tr 56 (V 6)t yVPsV 6) ) = if> ,4|t Y (-m.Ht> (4| = 0. 
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To find the effective action in 1 + 3 for massive spinors, we recognize that for 
the mass term we have 

56 56 56 56 

$Wiy°[-A* (+) f +B* HVp.M (+) +B [-W 4] = 

56 56 56 56 

V 4 V(--4* (+) f +B* [-?)(-m)(-A {+) +B [-W 4] , (6.10) 

56 56 56 56 

with s = 5,6, (±)t= — (=F) an d M^TtL while (*) means complex conjugation. 

56 56 56 56 

We took into account that y° (+)= - ( + ) y°, while y° [-]=[-] y°. We find 

56 56 56 56 56 56 

thatTr 56 ((+) t (+)) = Tr 56 (H) = 1 and Tr 56 ([-]*[-]) = Tr 56 (H) = 1- In order 

56 56 56 56 

that Jdx 5 dx 6 Tr 5 6(^ (4)t Y°(-^* (+) f +B* H t )Y*p.U (+) +B [-)W 4) will ap- 
pear in d = 1 + 3 as a mass term •4>' 4 'ty°(— mjtj)' 4 ', we must solve the equation 

56 56 56 56 

y s Vs (A (+) +B [-]) = (-vi)(-A {+) +B [-]). 

We can rewrite equations of motion in terms of the two complex superposi- 
tion of x 5 and x 6 : z := x 5 + ix 6 and z := x 5 — ix 6 and their derivatives, defined 
as Jj- := 7 (975- — i-afs")) ~§i '■— jijjz? + ^ef 5 ") an< ^ m terms of the two projectors 

56 

[±]:= 1(1 ± iy 5 y 6 ) as follows 

P) 56 3 56 56 56 56 56 

2iY 5 {^ H +— [+]}(A (+) +B [-]) = -m(-A (+) +B [-]). (6.11) 
oz oz 

Since in Eq.| |6.11| l t|>' 4 ' would be just a spectator, we skipped it. 

In the massless case the superposition of the first two states (i|} ( | 6 ' m ~° =(+) 

... . ... . 03 12 03 12 

^ , with ^ ' = (a +|3 K][-])^o) or the second two states 

u (6)m=0 ^)m= 0/ ^ ^(4)m=0 = (fl( +p (^[If])^) of the 

left handed Weyl representation presented in Eq. d6.66 must be taken, with the 
ratio of the two parameters a and |3 determined by the dynamics in x m space. 
In the massive case i|;' 6 > m is the superposition of all the states to which y 5 and 

56 56 

y° separately transform the starting state: \J)' 6 ' m = (A [+) +B [— ])4>± / with 

03 12 03 12 03 12 03 12 

x|4 4)m = fa[(+i)(+) ± [-!](+)] + p [t-i][-] ± (+i)[-]]}i|>o. The sign ± denotes the 
eigenvalue of y° on these states. 

56 

We shall therefore simply write (as suggested in Eq. l|6.7|l ) \j) (6) = (.4 (+) 

56 

+B [— ])i|)' 4 ' in the massless and the massive case, taking into account that in the 
massless case either A or B is nonzero, while in the massive case both are nonzero. 
Accordingly also 4>' 4 ' differs in the massless and the massive case. 

We want our states to be eigenstates of the total angular momentum operator 
M 56 around a chosen origin in the flat two dimensional manifold (M (2) ) 

M 56 =z|--z|-+S 56 . (6.12) 
oz oz 
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56 56 56 56 56 56 

Taking into account that y 5 (+)= - [-],y° (+)= - ( + ) y° andy 5 [-]=(+), 

56 56 

y° [—]=[—] y° (see Appendb j6.8t we end up with equations for A and B 

^ + > = 0- (6.13) 
For m = we get as solutions 

, (6)m=0 n f 5 , 6 -, , (4) 



n+1/2 

56 

-1/2) 



V2i = b ^ n H ^ , n > 0. (6.14) 



We required n > to ensure the integrability of solutions at the origin. The solu- 
tions have the eigenvalues of M 56 equal to (n+1/2) and —(n+1/2), respectively. 

Since in the massless case the contribution from (p 5 )2 compensates the one 
from (p 6 )2 for all the solutions from Eq. ||6.14| | with n > 1 and has therefore ob- 
viously one of the two contributions to the zero ml a negative real value unless 
n = 0, it seems natural to expect that the only massless solutions are the two 
solutions with the eigenvalues M.' 56 ' equal to 1/2 for the right handed spinor 

/2 — a o {+) A>+ ] ) and to 1/2 for the left handed spinor (4>L,/ 2 — 

bo [— ] 4'' 4 ')/ an d accordingly with the corresponding 4>' ) 4 ' m ~ and x| ) [ 4 ' m_0 G f 
the left and right handedness in d = 1 + 3, respectively. We shall reformulate 
the operator of momentum to be Hermitean on the vector space of solutions ful- 
filling the involution boundary condition in sect. 16.51 Having solutions of both 
handedness we must conclude that in such cases there is no mass protection. 
For m + 1 we get 

rr . 56 56 

*S?/2 = Qn(In (+) ~ iJ -+ ,ei * H^™*^™, forn > 0, (6.15) 

where J n is the Bessel's functions of the first order. The easiest way to see that J n 
and J n+ i determine the massive solution is to use Eq. d6.13l ), take into account that 
z = pe v *, define r = mp, p = \J (x 5 )2 + (x 6 )2, recognize that ^ = j£~ x ^[-§^ — 
i^r) and we find B = — J-M. Then for the choice A = I n e in * it follows that 

p ocp 1 xm oz ' IL 

B = -ie i ( n+1 '*(2.J n - which tells that B = -ij n +i e i(n+1 



6.3 Boundary conditions and involution 

In the ref . |2J we make a choice of particular solutions of the equations of motion 
by requiring that K4>| wa u = 0, where the wall were put on the circle of the radius 
Po of the finite disk (Eq. l|6.3| l). 

56 56 

This boundary condition requires that in the massless case (since [—](+)= 

56 56 56 ' 56 

while [-][-]=[-]) only the right handed solution (Eql6~Tl ^7 2 = Q o (+) 
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^ (4 ) m-o ^ a |. j g ^ e nanc j e( j w ith respect to S O ( 1 , 3 ) ) is allowed, while the left 
handed solution must be zero (b n = 0) making the mass protection mechanism 
work in d = 1 + 3. 

In the massive case the boundary condition determines masses of solutions, 
since only the solutions with T n+ i | p=Po = are allowed from the same reason as 
discussed for the massless case. This boundary condition determines masses of 
spinors through the relation m n+ i /2 po is equal to a zero of T n +i : 

Jrt+l (TTL n+1/2 Po) = 0. 

This time we look for the involution boundary conditions. 

First we recognize that for a flat M2 - {0} manifold, with the origin x 5 = = 
x 6 excluded, the Zz or involution symmetry can be recognized: The transformation 
P/Po — > (which can be written also as z/ po — > ^f) transforms the exterior of the 
disk into the interior of the disk and conversely. 

Then we extend the involution operator to operate also on the space of solu- 
tions 

d=(I-2£%/p ->P0/*. 

6 2 = I. (6.16) 

The involution condition & 2 — I requires, that H' is a projector 

(TZ>) 2 =1Z> (6.17) 

and can be written as W + H a dd, where 1t a dd must be a nilpotent operator 
fulfilling the conditions 

C^add) 2 = 0) "^add"^ = 0, HHadd = Hadd , (6.18) 
56 56 56 56 

We had H =[—], which is the projector. Since we find that [—](—)=(—] (see Ap- 

56 56 56 

pendix [oTBl , while (—)[—]= 0, we can choose 7t a dd = a (— ), where a is any 
function of z and Let us point out that H a dd is n °t a Hermitean operator, 

56 56 

since (— )t= — (+) and z^ = z, (-g^)* = Accordingly also neither 7c' nor & is a 
Hermitean operator. 

We now make a choice of a natural boundary conditions on the wall p = po 

{6l\> = iWlwall, (6.19) 

saying that what ever the involution operator is, the state 4> and its involution 
(fij; must be the same on the wall, that is at p = po- 

It is worthwhile to write the involution operator (9 and correspondingly the 
projector it' in a covariant way. Recognizing that n (p ' a y a n {p \, p b = [e 21 -* j(p 5 — 

56 56 

i-P 6 ) + jiv 5 +ip 6 )] H +ij(v 5 +T-P 6 ) + e~ 2itf 4(p 5 +ip 6 )] (+), we may write 

56 56 

1(1 -in(pLn(*W Q y b )(1 - |3n(P> Q Y a n(PU b ) =H (I + |3i[e 2i *£ + (-))• 
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This is just our generalized projector it', if we make a choice for a from Eq. d6.18b 

56 56 56 

as follows: a = (3i[e 2i *^ + ^ ( since HH=H)- We then have 

56 a a 56 

7t'= [-] (1+ [3i[e 2i *^- + — ] (-)), (6.20) 
oz oz 

where |3 is any complex number. 

6.4 Current through the wall 

The current perpendicular to the wall can be written as 

56 56 

56 56 

) ± = - 7 °{e-^ (+) (-)}. (6.21) 



We need to know the current through the wall, which for physically acceptable 
cases when spinors are localized inside the disk (involution transforms outside 
the disk into its inside, or equivalently, it transforms inside the disk into its out- 
side) must be zero. We find for the current 

{^UtWIwall = W+dltxdlWUl. ( 6 - 22 ) 

56 56 

Since 0t = I - 2(7$ + W add ) and 7t add = (a (-))t = -a* (+), it follows that 

56 56 

dt£id = -u -2a*T°e i4> [+] -2«r e- i * [+]. 
It must then be 

~ 56 

WH^Hwrfi = + 2y°(a*e i * + ae" 1 *) [+]}^)| wa n. (6.23) 

First we check the current on the wall for the "old" case, when a — and 3 — 

56 

I — lH, H —[—]. Not to be in contradiction with Eq. d6.23b the current on the wall 
must for either massless or massive case be zero. In the case of massless solutions 

(Eq. llOl ) only can fulfill this boundary condition (t|) (4 ' m =°tz n (-) (-) 

56 56 56 

{— Y°(e- 1 * (+) +e icf5 (-))}z n (+) il)V ] )lwail = 0, for each nonnegative n. The 
chosen boundary condition accordingly allows only the right handed solutions. 
We shall conclude when discussing Hermiticity of the operators that only n = 
is the physically acceptable solution. 

In the massive case the solutions of equations of motion (Eq. fl6.15l l) contribute 
no current through the wall, if T n+ i | wa n = 0, which is exactly what the boundary 
condition (Eq. fl6.19l l) (9i|)| wa u = ^Iwaii required. 

56 

Then we check the general case with & = 1—27?', where it' = H+it a dd =[—] 

56 

+ |3i[e 2l *-g^ + ■$=] (— ). For massless solutions it is not difficult to see that for any 
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nonzero choice of |3 only one handedness - the right handed one - survives and that only 
n = is allowed. In the massive case we find 

MP [e 2i * M + M] + B} | wall = o. (6.24) 
oz oz 

Since equations of motion require that B = — ^^jf and since -gt = je 1 *(^-^ + 
^■), we fulfill the involution condition on the wall for A — J n e in * only if B — 
-ijn+i e i(n+1 with the requirement that J n | waU = and (3 = ^. While J n | waH = 
can always be fulfilled, the second requirement (3 = means, since (3 can not be 
an arbitrary number, that our generalized condition is not written in an covariant 
form, and is accordingly not the acceptable boundary condition. 



6.5 Hermiticity of operators and the orthogonality of solutions 

In this section we comment on the Hermiticity properties of the operators, in par- 
ticular of p s and on the orthogonality properties of those solutions of the equa- 
tions of motion which fulfill the involution boundary conditions. We expect the 
solutions 

i) to be orthogonal (J d 2 xih f ((p 5 )2 + (p 6 )2)ij; j = J d 2 x^ i tTj) j m26 ij ) and that 

ii) on the space of these solutions the operators p s are Hermitean and have ac- 
cordingly expectation values of the operators (p s )2 positive contribution to ml 
for each s. 

Let us first check the orthogonality relations of the massive and massless so- 
lutions. We immediately see that the massive solutions ^^f+V/a belonging to dif- 
ferent n are all orthogonal due to the orhogonality of the functions e m *. We find 
Sd 2 xTr 56 (^l™} 2 ^™ 2 )=b n * <a n l^l 4 )^(4)m J Po (J;Tn + T;+iJn+i)pdp _ 

It also turns out that the massless solutions (^|f+i/2° (Eq. ^6.14t ) are orthog- 
onal to all the massive states (Eq. l|6.151 ) due to the properties of the J n Bessel's 
function. Namely, 



j2 t i .1 (6)m=0f | (6)m <- o* n ' 1 i.(4)m=0t 1 i,(4)m 

d xir 56 (\|; n+1/2 tp k+1/2 ) - 6nka n+1/2 ak^|^ <W 



Po 



p n J nP dp = 0, 



since Jq" p n jTipdp = p£ +1 J n +i (po), but J n+ i (po) must be zero in order that the 
massive state with n + 1 /2 obeys the involution boundary condition. Massless 
solutions are again due to the e ln * part orthogonal among themselves. 

So we conclude that all the states, which obey the equations of motion and the invo- 
lution boundary condition, are orthogonal. 

Are p s Hermitean on the space of these solutions? 

We know that p s = — i-gfr is Hermitean on the vectors space i|h if for any two 
functions i|h and from the vector space of solutions Tr 56 (J d 2 xp s (t|j{-4>j)) = 
(since then J dtxty^Vs^] + J d 2 x(-p s x|) i ) t ijjj = 0). 



We find that 



3 /cos cb it sincb-^-'. 
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for s = 5 (first row) and 6(second row). Since either messless W'n+V/a ' Eq. ^6.14) ) 
or massive W'jf+T/i' Eq. 16.15i ) states can be written as a product of e in * and the 
rest, say il) n; we see that J" d 2 xp s (i^ni^k) is nonzero only if |tl — k| = 1 . 

In this case we get that the integral Tr 56 (J d 2 xp s (i|) n 1 'tL) n±l )), s = x 5 ,x 6 , pro- 
portional to in{\ )|(n|4,Ata.+i | Po , with |pi|4*l>n.+i IS equal to 

i) a]^ =0 *aJ^p| (po) 2ln+1 ' in the case that two massless states are concerned, 

ii) a™ =0 *a™ +1 p^ +1 I n+ i (po) in the case that one massless and one massive state 
are concerned, 

iii) a^*a^ +1 po(JnJn+i + Jn+i Jn+2)lp m the case that two massive states are con- 
cerned. None of these integral is zero, since the two J n and J n+ i are not correlated 
( J n and J n+ i are correlated, if both belong to the solution with the same mass, de- 
termined by J n+ i (rap = mpo ) = 0). We conclude that for none of the solutions p s 
are Hermitean operators. 



One can check, however, that p" 



s 



3 1 X s 56 

0.=i{jr T -:T-5(p-Po)[+]}, (6-26) 
3x s 1 p 

are Hermitean operators on the space of massive and massless solutions, fulfilling the in- 
volution boundary conditions. It contains the part with the 5 function which corrects 
those parts of solutions, which are nonzero on the wall — the radial parts which 

56 

appear with ( + ) . It can be shown that the integral over the part with the 6 ( p — po ) 
function contributes just the terms which compensate the nonzero contribution 
in each of the three cases i)-iii). 

What we must check now is, what appears in this new definition of the 
operator of the momenta CEq. d6.25b '> as Y s PsY*Pt and whether now the integral 
Tr56(J d 2 xi|^ Y s p s y t Pt 4>), s = x 5 , x 6 , is still manifesting as just the mass term for 
those tJj which we accept as solutions of equations of motion (Eq. l|6. 1416.151 ). 

One finds 

Y^sY^t =p sP s 

13 1 9 i 3 56 

+ ^{[^6(p-Po) + -S(p-Po) + S(p-po)(^ rr]] W 

2 op p 3p p 3cp 

3 i 3 56 

+ 6(p-Po)(f- + -£-) H}. (6.27) 
3p p 3cp 

One notices that the first row of Eq.| |6.27| | represent the usual momentum squared. 
The last two terms are zero everywhere except on the wall. What we must check is 
the integral of the last two terms for all solutions fulfilling our involution bound- 
ary condition. 

56 56 

We find that the integral Tr56(J d 2 x (z (+))^ Y s PsY t Pt z n (+)] is for massless 
solutions (Eq. l|6.14b ) obeying the involution boundary condition proportional to 
n and it is zero only for n = 0. 

56 56 

The integral Trsg(J d 2 x (z (+))* Y s PsY t Pt z n (+)) demonstrates for massive 
solutions (Eq. (|6.15l l) the mass term squared originating in the first row of Eq. ll6.27l l, 
while the rest contributes zero. 
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The requirement that the integral Trs6(J d 2 x\|;^y s p s Y t pt 4>), s = x 5 ,x 6 , must 
be zero for massless solutions, makes a choice of only one among all possible massless 
solutions: the one with n = 0. 

Our the only possible solution is in the massless case "i^i/i ~~ ■ F° r the mas- 

56 _ 56 

sive solutions we have 4>i/2 = a i/2^j{Jo (+) — iTie 1 * [— ]), with m^Po as a 
zero of Ji,^™ 2 = a 3/ 2^(Ji -ifce 1 *)^*), with m. 3 / 2 Po asazeroof ]z,^_y 2 = 

56 56 

a_ 1/2 -^=(Ji (+) — ijo [— ] e _l *))e -1 -*), with ra_i/2Po equal to a zero of Jo and so 
on. 



6.6 Properties of spinors in d = 1 +3 

To study how do spinors couple to the Kaluza-Klein gauge fields in the case of 
M (1+5) , "broken" to M (1+3) x a flat disk with p and with the involution bound- 
ary condition, which allows only right handed spinors at po, we first look for 
(background) gauge gravitational fields, which preserve the rotational symmetry 
on the disk. Following ref . |2J we find for the background vielbein field 



6% e m a = 



f a - 

> 1 a — 



6V f'n 

= f^ s f» 



(6.28) 



with f a m — A^S^ m e. a r x r and the spin connection field 



'st\i 



(6.29) 



The U(1) gauge field depends only on x^. All the other components of the 
spin connection fields are zero, since for simplicity we allow no gravity in ( 1 +3) 
dimensional space. 

To determine the current, coupled to the Kaluza-Klein gauge fields A,^, we 
analyze the spinor action 



S 



d d xEil> (6) Y Q PoaV 6) = 
d d x4> (6 ' Y m (-)S sm cu sm ^ (6 



d d x1 |,(6) Y m 6 ^ p ^l 



6) 



,(6) 



d d xl J ) (6) y m 6 ^ mA ^ (£( T TX T pff + S 56^(6)_ 



(6.30) 



4>' 6 ' are solutions of the Weyl equation in d = 1 + 3 . E is for f a a from l|6.28b equal 
to 1. The first term on the right hand side of Eq. J6.30[) is the kinetic term (together 
with the last term defines the covariant derivative pojx in d = 1 + 3). The second 
term on the right hand side contributes nothing when integration over the disk is 



performed, since it is proportional to x CT (co sm|x 
We end up with 



-F e x a } 

2 1 irv u m^sff^ )• 



d 2 xii) (6) y m 5^mM 56 ^ (6) 



(6.31) 
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as the current in d = 1 +3. The charge in d = 1 + 3 is proportional to the total 
angular momentum M 56 = L 56 + S 56 on a disk, for either massless or massive 
spinors. 

6.7 Conclusions 

In this paper we were looking for what we call a "natural boundary condition" — a 
condition which would, due to some symmetry relations, make massless spinors 
which live in M 1+5 and carry nothing but the charge to live in M' 1+3 ' x a flat 
disk, manifesting in M' 1+3 ', if massless, as a left handed spinor (with no right 
handed partner) and would accordingly be mass protected. The spin in x 5 and 
x 6 of the left handed massless spinor should in M' 1+3 ' manifest as the charge 
and should chirally couple with the Kaluza-Klein charge of only one value to the 
corresponding gauge field, in order that after the second quantization procedure 
a particle and an antiparticle would not appear each of both charges. 
We found the involution boundary condition 



{dil>=iWUa, = (I-(I-m (p) a n<*> bY V')U^£o, o 2 =l, (6.32) 



which transforms solutions of the Weyl equations inside the flat disk into out- 
side of it (or conversely) and allows in the massless case only the right handed 
spinor to live on the disk and accordingly manifests left handedness in M' 1+3 '. 
The massless solution carries in the fifth and sixth dimension (only) the spin 1 /2, 
which then manifests as the charge in d = 1 + 3. 
We defined a generalized momentum p s 



which is the Hermitean operator in the case of our involution boundary condi- 
tion. 

The requirement that Y s f) s y t 'Pt manifests as the square of the mass leads 
in the massless case to to the solution with the total angular momentum 1/2 as 
the only solution, while the massive solutions carry all half integer angular mo- 
menta: ±1/2, ±3/2, • • • . The angular momenta in the fifth and sixth dimensions 
then manifests as the charge in the 1 + 3 dimension. The massless solution with 
the spin 1/2 is mass protected and chirally coupled to the corresponding Kaluza- 
Klein field. 

The negative charge of the massless 1 /2 charge state appears only after the 
second quantization procedure in agreement with what we observe. 

All the solutions fulfilling the involution boundary conditions are orthogonal 
and in this vector space and the generalized operators are Hermitean. 

The involution boundary condition of Eq. (|6.32l l are equivalent to the bound- 
ary condition, which we present in the ref . [ 2 1 . Both take care that massless solu- 
tions of one handedness appear in d = 1 + 3. 




(6.33) 
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We were looking for generalized boundary conditions with 

d = l-2lt>, 

56 r) r) 56 

£' = [-] U+[3i[e 2l *- + -] H), (6-34) 
oz oz 

where (3 is any complex number. This generalized boundary 7c' can be written in 
a covariant way as 

7c' = 1(1 -in (p Ln ( *W Q y b )(1 - |W p Q y Q n (p b p b ) 

56 a p\ 56 

= [-] (l + pi[e 2i *- + -] (-)). (6.35) 
oz oz 

But while in the massless case the generalized boundary condition {(fHj) = 4>}lwaii 
forbids all but s = 1 /2 solution, it fails in the massive case to demonstrate the 
covariance and is accordingly not an acceptable boundary condition. 

6.8 Appendix: Spinor representation technique in terms of 
Clifford algebra objects 

We define) 9 1 spinor representations as superposition of products of the Clifford 
algebra objects y a so that they are eigen states of the chosen Cartan sub algebra 
of the Lorentz algebra SO(d), determined by the generators S ab = i/4(y a y b — 
y b y a ). By introducing the notation 

(±i): = l(Y a TY b ), &=1(1 ±yV), for Ti a V b = -1 , 

ab 1 ab 1 

(±): = -(y Q ±i Y b ), [±]:=-(1 ±iY a Y b ), forn aa T 1 bb = 1, (6-36) 

it can be checked that the above binomials are really "eigenvectors" of the gener- 
ators S Qb 



ab V ab ab V ab 

S ab (k): = - (k), S ab [k]:= - [k] . (6.37) 



Accordingly we have 



03 1 03 1 

(±i):= 2(Y°TY 3 ), [±i]:= 5 (1±M 

(±):= ] -Iy' ±iy 2 ), [±]:=l(1±iyV) 

56 1 56 1 

(±): = -(y 5 ±iy 6 ), [±]:= -(1 ± iy 5 y 6 ) 



(6.38) 



03 03 12 12 

with eigenvalues of S 03 equal to ±j for (±i) and [±i], and to ±i for (±) and [±], 

56 56 

as well as for for (±) and [±]. 
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We further find 

ab ab ab ab 

y a (k) =r, QQ [-lc], y b (k) = -ik [-k], 

ab ab ab ab 

y a [k] = (-1c), y b [k]= -ikn aQ (-k) . (6.39) 

We also find 

ab ab ab ab ab ab ab ab ab ab 

(k)(k)=0, (k)(-k)=r, aa [k], [k][k]=[k], M[-k]=0, 

ab ab ab ab ab ab ab ab ab ab 

(k)[k]=0, [k](k)=(k), (k)[-k]=(k), [k](-k)=0. (6.40) 

To represent one Weyl spinor in d = 1 + 5, one must make a choice of the 
operators belonging to the Cartan sub algebra of 3 elements of the group SO(l , 5) 

S 03 ,S 12 ,S 56 . (6.41) 

Any eigenstate of the Cartan sub algebra (Eq. J6.41l l) must be a product of three 
binomials, each of which is an eigenstate of one of the three elements. A left 
handed spinor (r (1+5) = —1) representation with 2 6/2_1 basic states is presented 

03 12 56 

in Eq. d6.6|) . for example, the state (+i) (+)(+) iho , where \bo is a vacuum state (any, 
which is not annihilated by the operator in front of the state) has the eigenvalues 
of S 03 , S 12 and S 56 equal to 4, A and j, correspondingly. All the other states of 
one representation of SO(l , 5) follow from this one by just the application of all 
possible S ( ab), which do not belong to Cartan sub algebra. 
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7 How Can Group Theory be Generalized so 
Perhaps Providing Further Information About Our 
Universe? 

R. Mirman 
14U 

155 E 34 Street 
New York, NY 10016 
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Abstract. Group theory is very familiar, perhaps too much so. We are thus prejudiced 
about it, leading to views that are far too narrow. Yet it is significantly richer than usually 
realized (]13 [). Here we wish to understand the restrictions giving the familiar forms and 
how by changing these we can get added richness. Might these add to our knowledge of 
nature? A purpose of this note is to stimulate thinking about this. 

7.1 Geometry, through its transformations groups, is very 
information, but so far not enough 

It is clear that much (all?) of physics is determined by geometry, especially through 



its transformation groups (ED; ffl; [5|; [7]; |8|; \9\; [10); El; El; [3 1; |6|). Yet it is 



necessary to go much further. Can additional progress be made using group the- 
ory? This is a very open question but worth exploring. One aspect to be explored 
is whether group theory itself can be generalized. That could be of interest for 
purely mathematical reasons. And it has many applications. Generalizing it can 
thus be useful in various ways. This we wish to explore here. 

7.2 What is the best way to try to understand fundamental 
physics? 

What is the best approach to try to understand physics? The big fad nowadays is 
to come up with the wildest, most unlikely ideas, ones furthest from reality, ones 
totally unrelated to anything known, ones having no rationale whatever. History 
and common sense show that this approach is destined to lead nowhere except 
to even more wild ideas (as it has). 

Those who do that will find themselves badly cut by Occam's razor. Unfortu- 
nately a large part of the physics community is doing just that. Applying Occam's 
razor to the physics community will greatly help physics to advance. 

Another approach that is very likely to lead to failure, certainly if there is no 
other rationale for it, is to base laws on how we measure, on ourselves. We do 
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not determine the laws of nature (something many scientists, especially physi- 
cists, do not believe). How we measure is limited by physical laws, but does not 
limit them. Studying measurement can help us understand physics, but cannot 
determine it. 

What then shall we do, what approach shall we use? The best approach is 
the most conservative using requirements that are certain, or at least likely, to be 
correct, or ones that deviate the least from these. 

7.3 Reasonable requirements for developing theories 

What requirements can we impose? 

First is consistency. Fundamental physical theories must be consistent. (Phe- 
nomenological theories, classical physics is an example, can be inconsistent mish- 
mashes.) This is more difficult than it might seem, so can be quite powerful. 

Geometry imposes requirements, restrictions. Physics takes place in geom- 
etry so must be in accord with the rules it leads to. This also is powerful as we 
have seen (particularly) in the references. 

What can we say about geometry? We always assume that it is a manifold 
(locally flat) and that its coordinates are real numbers (rather say than complex 
numbers or quaternions). It is very unlikely that physics would be possible oth- 
erwise, but this can be investigated. A fundamental property of geometry is its 
dimension. However it has long been known that physics would be impossible 
unless the dimension is 3+1 (|4J; [13]). This is required by consistency, illustrating 
its importance, for only with this dimension is a consistent physics possible. 

7.4 What can we say about geometry? 

Is space curved? The curvature of space is given by a function over it, the con- 
nection ([5]). Can every space that is a manifold be regarded as flat but with a 
function, the connection, so that all curved spaces can be reduced to flat ones 
with different such functions? This is an interesting question that we raise but do 
not try to answer here. Also (many) curved spaces can be mapped (in reasonable 
ways) into flat ones ([8]). Thus we consider only flat spaces, but these questions 
should be looked into. 

What properties do flat spaces have? Beyond the reality of coordinates and 
the dimension is a most fundamental property: the transformation groups of the 
spaces (which are not symmetry groups (|13|), although it is quite interesting 
that they are that also). For our space, apparently the only one in which physics 
is possible, the largest symmetry group is the conformal group ([8]), which has 
subgroups the Poincare group, its subgroup the Lorentz group and the subgroup 
of that, the rotation group (SO(3)). The last gives that angular momentum must be 
integral or half-odd-integral (|3|), illustrating how transformations limit physics. 
The massless representations of the Poincare group determine electromagnetism 
and gravitation (0). Clearly these are quite informative, but clearly insufficient. 
It is possible that the conformal group can also be quite informative but how is 
less clear. 



82 R. Mirman 



7.5 How might group theory be generalized? 

Can we go further? What we wish to do here is study whether what is known 
about group theory can be generalized. We are all too familiar with semisimple 
groups, like the rotation and Lorentz groups. But group theory is far richer, even 
for these groups ([5J; |8j). Perhaps it is richer than we realize. That is what we 
consider here. 

This may not have anything to do with fundamental laws. But it helps to 
understand group theory, decreasing prejudice and broadening our views, and 
produces interesting mathematical results. And they are likely to lead to useful, 
even important, applications. 

7.6 Indexed groups 

Start by considering the curve, x = rcos0, y = rsinO, which describes a circle. We 
move around a circle using the two-dimensional rotation group 0(2). By putting 
a constant into the representation matrix we can generate an ellipse. But what 
about, say, the curve x = r(cos0)3, y = r(sin0)3. What set of transformations 
moves along this curve and why don't they form a group? Clearly there is an 
identity, we do not have to move, and for every transformation there is an in- 
verse. Moreover the product of two transformations is a transformation; if we 
move from A to B and the from B to C, we can find a transformation from A to 
C. However the transformations are not associative. It is for this reason that they 
do not form a group. The transformation from B to C depends on where B is (it is 
in a sense history dependent, depending on the previous transformations). That 
is the operator going from B to C has a form that depends on B, unlike rotations. 
This causes associativity to fail. 

Thus for a circle the product of the transformation matrix for 0i and for 02 is 
that for 81 + 02, which is not true for this curve. 

While there are transformations along any (reasonable) n-dimensional sur- 
face it is only in special cases that they form a group (of the form usually consid- 
ered). This emphasizes the relevance of associativity and the restrictions it places. 
Many of the properties of groups and their representations come from associa- 
tivity. While restricting, it also allows us to obtain properties that are so useful in 
applications of groups. 

How do we deal these more general transformations, say ones along arbi- 
trary surfaces? We introduce the concept of indexed groups. To do this we assume 
that the transformations can be mapped (at least) one-to-one onto a group of trans- 
formations. For each general transformation we have a corresponding group ele- 
ment and this element is the index of the general transformation. Thus for three 
dimensions the index can be an element of the rotation group. For a group the 
matrix representing the transformation that is the product of two is the matrix 
product of the matrices of the two transformations. It is here that indexed repre- 
sentations differ. For these the matrix labeling the product is the matrix product 
of the two labeling matrices. But the product matrix of the indexed transforma- 
tion is not the product of the matrices of the two transformations it is a product 
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of. It is the index that is given by the product, not the transformation. With this 
definition of group product, differing from the normal definition, these transfor- 
mations form a group. Associativity holds, since it follows from the associativity 
of the indexing group, but only because of the revised definition of a product. 

To give a group we list its members and their products (thus the spaces on 
which they act). But now with each set of members we have an infinite number 
of product rules (determined by the mapping of the transformations into the in- 
dexing groups, of which there may be several) thus an infinite set of groups. Each 
n-dimensional surface has its own group. 

7.7 Product rules determine groups 

This shows how properties of a group are dependent on the definition of its prod- 
uct and how by revising this definition we can generalize the type of structures 
that form groups. This adds to the richness of group theory. 

Groups which can be realized as matrices, thus whose products are matrix 
products, we call standard groups. Ones whose elements are indexed so whose 
products are given by the matrix products of their indices we call indexed groups. 

Indexed transformation groups exist for any surface that can be mapped 
(properly, in a way that must be investigated) to the defining space of a Lie group. 
For the three-dimensional rotation group SO(3) that is a sphere, and there is a 
third parameter which can be considered as giving the direction of a vector at 
each point of the sphere. Then each point is mapped to a point on one of the 
generalized (overlying) space, and each direction to one on that (which perhaps 
might be considered as an internal symmetry). Likewise we can use SO(2,l) to get 
another set of such surfaces. So we have associated with each group an infinite 
set of groups, each given by different product rules (from different mappings), or 
another way of saying this, an infinite set of realizations or representations. 

We assign to each group element a matrix, that of the regular (adjoint) rep- 
resentation. That this is possible follows from the group axioms. Then the group 
product is a matrix product, and this is the usual group product. We call this the 
standard product. 

The rotation group, besides its defining representation of 3 X 3 matrices has 
an infinite set of others. Consider the 5X5 one, say. This is a subrepresentation of 
SO(5). We can map a surface to the defining surface of SO(5) and the group de- 
fined over it (one for each of the infinite number of such surfaces, ignoring aspects 
like inversions), form representations of SO(3), but with additional transforma- 
tions which might be taken as internal ones. Since SO(3) has an infinite number 
of representations (and it is simple so other groups have infinite sets each of in- 
finite numbers of representations) this admits a huge number of transformation 
sets to be defined over it. 

But we can do more. The conformal group algebra is (isomorphic to) that 
of SO(4,2) and SU(3,1). The group algebras are the same but realized in terms of 
different variables. Instead of 4+2 real ones, or 3+1 complex ones, the conformal 
algebra is realized over 3+1 real ones ([8]). We might also realize it over more, 
rather than fewer, variables. These again can be taken as internal ones. 
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We now map surfaces (choosing from an infinite number) to the 3+1-dimen- 
sional real space on which the conformal group acts. The product of the trans- 
formations on these is given by the product of conformal transformations taking 
points and directions of the 3+1-space to others. These conformal transformations 
are the indices of the transformations on the preimage space. Note that we have 
three groups, SO(4,2), SU(3,1) and the conformal group, plus all their representa- 
tions, which can act as indices. This shows the great richness introduced. 

7.8 Groups of 3+1 space as illustrations 

Thus there are two groups defined by a 3+1-dimensional real space, one, ISO(3,l) 
is the Poincare group, the other is the conformal group. This itself shows rich- 
ness known but not understood in group theory. Take a surface (one of an infinite 
number) that is mapped to real 3+1-dimensional space. The transformations on it 
can be indexed by the transformations of the base space. However there are two 
sets, the Poincare group and the conformal group. So from ordinary space we 
can find two infinite sets of groups (realizations, representations). And the trans- 
formations of our space can be taken as subsets of larger groups, giving further 
labels and products, of infinite number. Some of this additional freedom might 
be relevant to internal transformations. 

Why should we consider these, aside from their showing the assumptions? 
Groups are useful in many ways and these can extend their usefulness. For exam- 
ple special functions are group representation basis states and many properties 
can be derived from this. Generalizing the concepts of representation can lead to 
other special functions, perhaps with useful properties. However it is not clear 
that properties can be found for these as they can for standard representations, or 
indeed that they have simple properties. This must be investigated. Associativ- 
ity is important in determining these properties, and allowing simple properties 
(that we can get rules for). This procedure allows such great generalization that it 
is likely that only a few cases, at most, can give simple rules. But there might be 
some and these could be useful. 

Also physical objects are statefunctions (L13J) that are group representation 
basis states. By expanding the set of these we may be able to expand the set of 
objects that are such states. There are clear limitations as the known states are 
those of standard representations. It may be that the requirement that objects be 
observers, and conversely ([7\), provides strict limits. Yet this is not known and 
these new representations allow study of this. And some may have physical ap- 
plications, perhaps to these fields. 

7.9 Why standard products are matrix product and why are 
these usually relevant? 

While these indexed groups (groups with indexed products) may seem unusual 
they raise the question why standard products are the relevant ones, for those 
cases in which they are? This has to be considered for each application. A gen- 
eral class is applications to geometry. The transformations, but only for certain 
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geometries, have standard products. This is true for lines, circles, planes, spheres 
and generalizations. For these, why are the products of transformations the stan- 
dard products? 

Here symmetry enters. The action of a group transformation on the base 
space is independent of the point in that space — since all these points are iden- 
tical. Thus a group transformation taking a point to another, acting on a second 
such transformation, gives a transformation with identical action (as can be seen 
with a circle), i From this associativity follows. These transformations thus form 
the regular representation ([3|). But this representation can be given by matrices, 
and the group product is the same as a matrix product, the standard product. 

Other spaces do not have symmetry so their transformations cannot be rep- 
resented by matrices. Here we see how symmetry gives group operations, and 
limits these. For spaces without symmetry we must use other products. 

7.10 Conclusion 

Groups are determined very much by their products. Usually these give matrix 
products. Here we have considered a set of different product rules. These illus- 
trate how such rules determine the properties of groups, and the role of symme- 
try in the standard product rules. Whether some of these generalized rules are 
useful has to be studied. Lie groups have an extensive structure including their 
algebras. Do these generalizations allow corresponding structures, including al- 
gebras? That is another field of study. There is much that can be done, some at 
least profitably. 
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Abstract. We want to unify usual equation of motion laws of nature with "laws" about 
initial conditions, second law of thermodynamics, cosmology. By introducing an imagi- 
nary part - of a similar form but different parameters as the usual real part - for the action 
to be used in the Feynman path way integral we obtain a model determining (not only 
equations of motion but) also the initial conditions, for say a quantum field theory. We set 
up the formalism for e.g. expectation values, classical approximation in such a model and 
show that provided the imaginary part gets unimportant except in the Big Bang era the 
model can match the usual theory. Speculatively requiring that there be place for Dirac 
strings and thus in principle monopoles in the model we can push away the effects of the 
imaginary part to be involved only with particles not yet found. Most promising for seeing 
the initial condition determining effects from the imaginary part is thus the Higgs parti- 
cle. We predict that the width of the Higgs particle shall likely turn out to be (appreciably 
perhaps) broader than calculated by summing usual decay rates. Higgs machines will be 
hit by bad luck. 

8.1 Introduction 

Usually when we talk about "theory of everything" as superstring theory is hoped 
to be, it is not really meant that the initial state of the universe is included in the 
model immediately. Rather one needs to make extra assumptions - cosmology, 
second law of thermodynamics [1 2 3 4J, etc. - about the initial conditions or one 
simply leaves it for the applicator of the theory to somehow himself manage to 
find out what the initial conditions are for the experiment he wants to describe 
with the theory. It is, however, the intention of the series of articles B5I6I7I8I9I to 
which this article belongs to set up assumptions telling the initial conditions in 
a way that can be called that these initial condition assumptions are unified! 10 1 
with the part of the theory describing the equations of motion and the particle 
content (the usually T.O.E.). Our unification may though be mainly a bit formal 
in as far as our main point is to use in the Feynman path integral an action which 
has both a real Sr and an imaginary part Si. Usually of course the action is real 
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and the imaginary part Si = (Total S = Sr + iSi). We may quickly see that the 
imaginary part gives a typically hugely different extra factor in the probability 
for different paths obeying equations of motion. Thus such an imaginary part es- 
sentially fix the path obeying equations of motion which should almost certainly 
be the realized one. In this way we can claim that to a good approximation an 
imaginary part of the action will choose/ settle the initial conditions. 

In the present article it is not the point to settle on any choice of the in usual 
sense "theory of everything". Rather we shall present our idea of introducing an 
imaginary part in the Lagrangian and thereby also in the action as a modification 
that can be made on any theory as represented by the real action S. 

We have already published a few articles on essentially a classical formula- 
tion of the present model. We sought in these articles to be a little more general by 
simply defining a probability weight called P(path) defined for all possible paths. 
In classical theory it is really only the paths which obey the classical equations of 
motion for which we need to define P. We already in the earlier articles suggested 
that this probability P(path) for a certain track, path, to be the one realized in na- 
ture should be given as the exponential of an expression depending on the track, 
path, of the form of a space-time integral over a locally defined quantity C\ de- 
pending on the fields in the development, path. Really this quantity L\ (really 
— C\) comes into determining the probability as if it were the imaginary part of 
the Lagrangian density. 

A major point of the present article is to set up the quantum formulation of 
our already published model, now really settling on taking the suggestive idea 
of just making the action complex, but with a priori a different set of coupling 
constants and m 2 for real and imaginary part separately. 

A genuine problem with our kind of model is that very likely it predicts that 
special simple configurations leading to big probability may be arranged at a pri- 
ori any time. That is to say, with our type of model it needs an explanation that 
one in particle almost never see any great arrangements being organized to occur 
later on. Really such arrangements might seem to us to be something like a hand 
of God, but they seem very seldom. Thus at first it looks that our type of model is 
already falsified by the non-appearance of arrangements. Really such a problem 
is almost obviously expected to occur in a model that like ours does not a priori 
make any time reversal asymmetric assumption at the fundamental level. Unless in the 
Hartle-Hawking no boundary postulate lHUl we add some timereversal asym- 
metry spontaneously other otherwise that theory will be up to similar problems 
1TTTT21 . 

A model-language describing how final states can be imposed by a density 
matrix Pf is put forward by Hartle and Gell-Mann |[14|. 

In the present work we hope for that a certain moment in the 'middle of 
times' will turn out to become dominant w.r.t. fixing the special solution selected 
as the realized one, and that this time can then be interpreted as a close to Big 
Bang time ( there may not really be a true big bang but just an inflation era com- 
ing out of a deflation era continuously). Then since we live in the time after this 
decisive Big Bang simulating era there is for us a time reversal asymmetry, nev- 
ertheless it is a problem that like ours is even timetranslational invariant w.r.t. 
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the law that finally settle the 'initial conditions' to explain that there are not more 
prearranged events than one seemingly see. 

However, we believe to have found some explanations able to suppress so 
many of these prearrangements that our model can be made compatible with 
present experience of essentially no prearrangements. 

For really avoiding it we shall assume consistency of Dirac strings, but let us 
postpone that discussion to section l8A3l below. 

Our model is really inspired from the considerations of time machines[13| 
and the troubles of needs for prearrangements in order to avoid the so called 
grand mother paradoxes, meaning the inconsistencies occurring when one seeks 
to go back in time and changes the events there. 

We shall present the work by making two attempts to assumption about how 
to interpret the Feynman path integrals with the imaginary part of the action 
non-zero. In the first part of the paper we start out from letting the average of a 
dynamical variable O be given by equation ||8.10|| below, but that this is a priori 
not so good is seen by it not being (safely) real even if the dynamical variable O 
is real. Therefore in section W7\ we restart the discussion so now from the side of 
the interpretation of the Feynman path way integrals in our model. 

First trial of interpretation 

In the next section l8!2l we shall put forward the basic formula for expectation 
values with our complex action model and the philosophy that this model even 
deliver the initial condition, or better the solution of equation of motion to be the 
one realized. 

In section [8~3l we review our earlier reasons for that future should have only 
little influence on what happens. 

In section 18.41 we then shall argue for some approximate treatment of the 
functional integral in the late times t, the future. 

In section l831 we shall make use of the approximation of the future to obtain 
the usual quantum mechanics expressions at least in the case where our imagi- 
nary part Si of the action can be ignored. (It should be stressed that we actually 
have used already a philosophy based on this S i being non zero, so it is not fully 
zero.) 

In section l8l>l it turns out that we - perhaps not completely convincing though 
- can make the effect be that we return to probability in practical scattering ex- 
periments say get conserved. 

Second trial of interpretation: 

In section [8~71 we restart the discussion of making the interpretation formula 
for the Feynman path way integral, which after some talk takes the way of using 
the classical approximation weighted with the exponential of minus 2 times the 
imaginary part of the action. In a subsection l8.7.5l we formally connect our model 
to our earlier one based on the probability weight P(path). 

In section [8781 we develop a rather general formula for the correlated proba- 
bility for a series of dynamical quantities or operators Oi at different moments of 
time take values inside small ranges specified. 

In section 18.91 we go a bit further in making the expressions like the ones 
one uses in practice in usual theories. Most importantly we again consider how 
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to approximate the future when the effects of the imaginary part of the action is 
very small. 

In section 18.101 we put the simplest example of the more general formula, 
namely a formula for the probability of just one operator at one time being in 
a given range. ( This question would be impossible to predict even in principle 
in other theories, but we in principle can, but in practice not usually). But the 
resulting formula has what we call "squared form" in the sense that the projector 
comes in twice as a factor in it. The finding of a reduction to an unsquared form 
is left to section l8.121 while we in section 15 . 1 1 1 then give an example of application 
of very interesting physical significance. In fact section [8 . 1 1 I predicts a broadening 
of the width of the Higgs particle due to the imaginary action. 

In section [8.12l we then bring about a connection between the postulated in- 
terpretation formulas for probabilities put forward in part I and part II. In fact we 
find that they coincide under rather suggestive assumptions. 

In section 18.131 we bring the promised argument for removing the effects 
of the imaginary action Si from the domain of older accelerators, since other- 
wise our model would have been falsified. The argument is based on assuming 
monopoles. 

In section l8.14l we conclude and give a bit of outlook. 

Part I, First Trial of Interpretation 

8.2 Philosophy and formula 

Our basic modification of introducing an imaginary part in the actions leads 
to that integrand e lS or e^ s of the Feynman path way integrand Je^f^cf) or 
Je7r s <^cf) (jf the Planck constant is written explicitly) varies a lot in magnitude, 
and not only in phase as usual. This effect is likely to make some regions in the 
space of paths - or we could restrict to the space paths with 6S = 0, i.e. the space 
classical solutions - get a very much bigger weight in the integral than others. 
Actually it can likely happen that only a very narrow range of paths or better 
solutions (= paths obeying 6S = 0) will quite dominate the integral 

e iS ^4>. (8.1) 

That should naturally be taken to mean that the presumed narrow range of 
dominating paths represent the paths being actually realized in nature. It is in this 
way that we hope our model to essentially predict the initial state for the realized 
solutions. It is important to have in mind that such an effect of the imaginary 
action Si of selecting narrow bunches of solutions can make the boundary con- 
ditions at an initial and a final time for a period to be studied say, superfluous. 
A bit optimistically we might imagine that the imaginary part of action makes 
the functional integral converge even without boundary condition specifications. 
Note that being allowed to throw boundary conditions away - having them re- 
placed by effects of Si - is a great/nice simplification. We consider this achieve- 
ment as an aesthetically very nice feature of our model! Supposing that this works 
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to deliver a meaningful Feynman-path integral ll8.ll even without boundary con- 
ditions this way we must now decide how one is supposed to extract information 
now in principle for the true expectation value as it should occur even without 
further input. Note here that we are - but only in principle - proposing an ex- 
ceedingly ambitious model compared to usual quantum field theories: 

We want to predict expectation values without any further input than the 
mere complex action! This of course corresponds to that our level of ambition 
is to in addition to the usual time-development laws of nature also predict the 
initial conditions, i.e. what really happens! 

To write down the formula for some physical quantity let us first exercise by 
a quantity C(cp|t) which is a function of the fields cp|t restricted to some time t, 
where cp is a general symbol for all the fields in the model. 

If we for instance use the Standard Model as the starting model, then pro- 
viding it with an imaginary part of the Lagrangian density, then the symbol 
cp(x) (x 6 R 4 ) is really a set 

cp = (A^ b ,H) (8.2) 

where the indices on the Fermion fields runs through the combination of fla- 
vor and color and/or W-spin components, while the index on the gauge fields 
run through the 12 gauge fields - 8 gluon color combination plus (B H ) the U(l )- 
component and 3 W's -. Finally H is the two complex component Higgs field. 

The quantity 0(cp| t ) can of course be considered a functional of the whole 
field development 0(cp) also, i.e. it could be consider a functional of the path of 
one wants. 

The simplest proposal for what the average quantity 0(tp) would be 

This would mean that we have a "sort of probability" given by 

"Probability of ObeingOo" = J j^^cp ~ (8A) 

Now, however, we must admit that conceiving of this expression as a probability 
is upset by the severe problem that it will typically be a complex number. There 
is no guarantee that it is positive or zero. 

Thus a priori one would say that this simple expression for the probability 
density is quite untenable. 

Nevertheless it is our intention to claim that we should - and that is then part 
of our model - use the simple expression | |8.4| | and the corresponding i8.3\ and the 
expression to be given below for more general operators O corresponding also to 
JO and fp. 

First let us again stress that it is our a priori philosophy that somehow the 
imaginary part Si managed to fix both a state in future and in past. Thereby ask- 
ing the average of some quantity O becomes much like in an already finished 
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double slit experiment (Bohr-Einstein) in which a particle already have been mea- 
sured on the photographic plate (presumable on an interference line) what were 
the average position of the particle when it past the double slit screen. Really ask- 
ing such a question concerning a quantity O that were not measured and could 
not have been measured without having disturbed the outcome of something 
later is one of the forbidden questions in quantum mechanics. Indeed it is by ask- 
ing this sort of questions which are not answerable by measurement that Einstein 
can find ammunition against quantum mechanics. In other words our proposal 
| |8.4| | for "probability distribution" is a priori - with our present philosophy of 
a future essentially determined by Si - an answer to a quantum mechanically 
forbidden question. Niels Bohr would say we should not ask it. 

In that light it may of course not be so serious that our formula gives a rather 
stupid or crazy answer, a complex probability! 

But now we have the problem of justifying that if we made a true measure- 
ment the answer would turn out to give positive (or zero) probability. 

Let us take as the important feature of a measurement of some quantity O 
that there is an apparatus which makes a lot of degrees of freedom, £, say (really 
macroscopic systems) develop in a way depending on value of O. Such an ampli- 
fication of the effect of the actual value of O is characteristic for a measurement. 
Unless somehow there are special reasons for that Si be insensitive to £, (as we 
shall actually later seek to show but do not assume to be the case) we expect that 
Si typically will depend on the macroscopically many d.o.f. £, being influenced by 
O-value measured. Now we argue like this: Since there is a huge (macroscopical) 
number of variables £, depending on the value of O "measured", the imaginary 
part Si of the action is likely to depend very strongly on this measured value - 
very rapidly varying. 

We here think of Si as the integral over the imaginary part of the Lagrangian 
Li over all times t e] — oo, oo[ 



Because of the great complications in an actual measuring apparatus, let 
alone the further developments depending the measured value, publications and 
so on, the imaginary action Si can easily be a very complicated function of the 
measured O value. Even if Si as function of the measured O value should in 
principle be continuous it may in practice vary so much up and down - caused 
by accidents influenced by the broadcasted measuring value - that very likely the 
smallest value of Si occurs for a seemingly accidental value of the measured O. If 
the Si-variation with the "measured O" is indeed very strong so that the Si vari- 
ations are big the exponential weight e~ Sl contained in d8-3|> and | |8.4|| will have a 
completely dominant value for only one measured O-value. 

In this way our model has the integral in the numerator of l|8.4t be much big- 
ger for one single value of Oo - If so, then the ratio l|8.4t is actually oc 1 for this Oq- 
value and negligible for all other Oo -values. This means that our model much like 
usual measurement theory (in Copenhagen interpretation) predicts that crudely 
only one value of a measured quantity is realized. In principle it is even so that, 




(8.5) 
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bearing a very special situation, the result of the measurement is calculable by 
essentially minimizing the imaginary action Si. In practice, however, such calcu- 
lation will only be doable in extremely rare cases. (If we impress a special result 
by threatening with a Higgs-producing machine). 

We postpone the argumentation for that the probability distribution to be 
obtained in practice shall be the one of usual quantum mechanical measurement 
theory partly to the later sections and partly to a subsequent paper. 

At the end of this section let us extend slightly our formula (|8.3b and thereby 
also ||8.4| |, to the case where the quantity O corresponds in usual quantum me- 
chanics to an operator that do not commute with the fields cp. 

An operator corresponding to a quantity measurable at a moment of time t 
will in general in the quantum field theory considered be given by a matrix with 
a columns and rows in correspondence with field functions cp[t restricted to the 
time t. I.e. O is given by a "matrix" 

(<p'lt|0|<plt) =0(«p'lt, (Pit). (8.6) 

What should be the formulas replacing 48.3b and J8.4I I in this more general 
case? 

Well, our main starting point were that we assumed our imaginary part Si 
to (essentially) fix both a further |B) and a past state |A). A natural notation to 
introduce is in fact - for the past - 



(cp| t |A> = A[cp| t ] 
and analogously 

(B|cpl t > = <cpltlB>*=B[cp| t r = 



e^-ootot^q, (8.7) 

ending at cp 1 1 



beginning at tp 



In this notation our previous formulas (|8.3b and J8.4II are for the in cp| t diag- 
onal operators 0(cp|t) become 

JV s O(<p| t )ff<p 
[ ' Je iS ^cp 

^ |t <B|<p| t )Q((p|t)(cp|t|A) 

^(BMOWtlA) 

(BIOIA) 



(B|A) 



(8.9) 



and d8.4b becomes 



"Probability of ObeingOo" = {m °~° o)lA) (8-10) 

(B|A) 

Now really we want to suggest that formula (|8.9t and J8- 10b can also be used 
for operators that are not simply functions of the fields <p|t at time t, used in the 
functional integral. 
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In order to justify that extension of our interpretation formulas we want to 
remark: 

1. Provided a Hermitean operator O has either |A) or |B) as eigenstate then the 
eigenvalue O ' in question can of course be extracted as 



_ (BjgjA) 
(B|A) 



(8.11) 



2. One can quite generally - by Fourier transformations at every step in a time 
lattice - rewrite a functional integral of the Feynman path way integral form 
from some set of variables cp to a conjugate set: 



e tS ^cp 



n 

te{t— lattice] 



^(3) (p| t e* ^- tSUatHce] 



where 



n ^(<pit+At,(Pi t )^ (3) cpk 

te{t— lattice} 



(cplt+At.cpIt) ^e^i'" 1 - ™ ), 



L dis „(c P i t , qjlt ^r' p ' t ) At 

(8.12) 



(8.13) 



can be rewritten into ^(TT|t+At,n| t ) matrices obtained from the ^(<p|t+At, <plt) 
by Fourier functional transformations 



^(n| t+ At,n| t ) d =i f 



(3) cplt+At e +i(plt+Atn|t+At ^((p|t+At,cplt)e" iipltn|t ^ (3) (plt 

(8.14) 



Now of course for long chains of W -matrices (ignoring end problems) you have 



except for end problems 



]1 ^(^lt+A t ,cp|t)^ (3) cp| t 

te{t— lattice} 

]1 ^(nit+At,n| t )^ ,3 »n| t 

te{t— lattice} 



(8.15) 



Supposedly you can put the right hand side into a form 

3 iS (in n) [n, An-defferences] 

Now you may argue with the same intuitive suggestion for getting 

° (n|t) - Je«"-"'^T 

as we did for l|8.3|l . By thinking of doing the just presented Fourier transformation 
partly we might argue for a similar average formula for any operator 

JV s 0(cp| t ,TT| t )^cp 



(8.16) 



(8.17) 



(0(<p| t , n_t)) = 



Je iS ^cp 
(B t |Q(<p|t,TT|t)|At) 
(Bt|A t ) 



(8.18) 
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Really this proposal looks very bad because of several lacks of good correspon- 
dence with usual quantum mechanics a priori: 

a) Obviously |A t ) is here (a sort of) wave function of the universe at time t, but 
our probability density ( 18.4b or 

"Probability for O being O " = ( B ^°-°°^) (8 . 19) 

(B t |A t ) 

is not quadratic in |A t ) as we expect from the usual corresponding formula 
"Probability for O being O usual" = ^ At|6( ^ . (8.20) 

b) As already stated the "probability density" | |8.19| | is even usual complex and 
needs the above measurement special case to become just positive. 

We shall below argue for an approximate treatment of the future part |B t ) of 
the integral thereby achieving indeed a rewriting into an expression which is of the 
form with |A t ) coming squared. Indeed we shall rewrite J8. 19b into l|8.20b below. 



8.2.1 Justification of philosophy from semiclassical approximation 

In semiclassical approximation one simply evaluates different contributions to 
the functional integral ( 1 ) by seeking the different extrema for e lS or equivalent 
S = Sr + iSj. Around such an externum it is extremely well known that one can 
approximate S by the Taylor expansion up to second order 



S = S(extremum) + - 



9 2 S 



3cpi (xi)9(p 2 (x2) 

c Pl (x 1 )-cpf r (x 1 ))(( P2 (x 2 )-cpf r (x 2 ))+-.- d 4 x,d 4 x 2 (8.21) 



where then the linear terms 

3S 



Qcpi (xi 



(qMxO-cpfVOKfxO (8.22) 



vanish because of the extremiticity condition. Here (pf tr (xi ) and (p| xtr (x 2 ) denote 
the fields at the extremum field configuration development. Such an extremum 
as is well known corresponds to a solution to 

6S = (8.23) 

i.e. solving the variational principle leading to classical equations of motion. 

The main term in the exponent iS (extremum) is in the usual real action case 
purely imaginary and thus only gives rise to a phase factor so that in this approx- 
imation the contribution has the same size for all the classical solutions, provided 
they can go on for real field configurations. With our Si included, however, we 
tend to get even to the approximation of the first term in the Taylor expansion 
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1 182]) a real term —Si into the exponent and thus the order of magnitude for one 
classical solution compared to another can easily become tremendous 

| giS (externum) g—S : (externum) (8 24) 

It is our philosophy that only relatively very few classical solution have terms 
e -Si (extremum) dominating violently the rest. In this sense we expect and assumed 
that such one or a very few classical solutions could be considered the only one 
realized. With very big size of Si - and that can easily come about for a couple 
of reasons - it gets relatively only exceedingly few classical solutions that are 
competitive in the sense that for most classical solutions (of l|8.23t ) you have ex- 
ceedingly small e~ Sl compared to the few dominant ones. As the reasons for Si 
being big when it is not forbidden by gauge invariance and the condition that 
Dirac strings shall be unobservable we can give: 

1. There is in analogy to the SR-term a |- -factor in front of Si. For practical pur- 
poses we know that we shall consider the Planck constant h to be very small. 

2. We could easily get Avogadros number come in as a factor in the Si because 
it would get such a factor a priori since there are typically in the world of 
macroscopic bodies of that order magnitude molecules. 



8.3 Approximate treatment of future part of functional integral 
(treatment of |B t )) 

In our earlier works|8J - in which we mainly worked in the classical approxi- 
mation - we presented some arguments that in the era which have been going 
on since short time of after some effective (or real) Big Bang the imaginary La- 
grangian or action l_i or Si effectively became very trivial. That should mean that 
under the times starting after some early Big Bang and extending into the future 
we could approximately take I_i and the part of Si coming from this era as inde- 
pendent of what are the practical possibilities for what can go on. Thus we should 
in this present era supposed to extend into even the infinite future be allowed to 
ignore in first approximation the imaginary parts Li or Si. 

The reasons, which we presented for that were that this present era including 
supposedly all future is dominated by two types of particles: 

1. Massless particles (really the entropy of the universe is today dominated by 
the massless microwave back ground radiation of photons). 

2. Non-relativistic particles carrying practically conserved quantum numbers 
(the nucleons and the electrons are characterized by their charges and baryon 
or lepton number so as to make their decays into lighter particles impossible). 

The argument then went that we could write the action - actually both real 
Sr and imaginary Si - for these particles, treated as particles, as a sum having 
each giving a contribution proportional to the eigentimes for them: 

Sr, Si= V Kp r r i • Tp . (8.25) 

particles P 
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That is to say that each of the particles contribute to Si say a contribution 
proportional to the eigentime 

Si from P oc Tp. (8.26) 

Now for massless particles any step in eigentime 

At p = (formassless) (8.27) 

and for nonrelativistic (~ slow) particles, such a step is 

At p = At (8.28) 

equal to the usual time. Since the number of the conserved quantum numbers 
protected particles are all the time the same the whole contribution to the Si from 
the present era becomes very trivial: 

Zero from the massless, and just a constant integrated over coordinate time 
for the conserved particles. 

In addition there are terms from interactions contributing a priori to say Si 
also. Since, however, in the era since a little after Big Bang the density of parti- 
cles were low in fundamental units presumably also the interaction contributions 
would be much suppressed in this after Big Bang era. 

So all together we estimate that it is only the very early Big Bang times that 
will dominate Si. Thus the solution to the equations of motion being in a model 
with an imaginary action Si selected to be the realized one will mainly depend 
on what happened in that solution in the early Big Bang era. This means that it 
will be in our era as if it were the initial state that were a rather special one de- 
termined by having an especially small contribution to Si from Big Bang times. 
This would mean a rather well determined starting state roughly which inter- 
preted as a macrostate would be one with low entropy. That is at least a good 
beginning for obtaining the second law of thermodynamics, since then there are 
supposedly no strong effects of Si any more to enforce the universe to go to any 
special macrostate. Rather it will go into bigger and bigger macrostates meaning 
that they have higher and higher entropy. 

Although we have now argued for approximately seeing no effects of Si in 
the era after Big Bang implying that our model should have no effects in this era, 
this is however, presumably not being quite sufficiently accurate. 

We shall, however, below in section 18.91 invent or find arguments that will 
allow us to get completely rid of the I_i or Si from the in the Standard Model 
already found particles. Only for the Higgs involving processes our arguments 
in section 9 based on gauge symmetry and the assumption of unobservability 
of Dirac strings associated with monopoles will not quite function. Thus we still 
expect that an Si -contribution pops up with Higgs-particles. But since Higgs- 
particles are so far not well studied such an effect of Si might well have been 
overlooked so far. 
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8.4 Treatment of |B t ) or Treatment of the future factor in the 
functional integral 

In equation d8.8|> above we defined what one could call "the future part" of the 
functional integral relative to the time t. It should however be kept in mind that 
it is a part in the sense that the full integral is a contraction (a sort of product) of 
the past part and this future part, 



Now we must remember that according to the second law of thermodynamics 
the state of the universe if at all obtainable (calculable) should be so by consider- 
ing the development in the past having lead to it. The future, however, should be 
rather shaped after what happened earlier. This suggests that we should mainly 
have the possibility to guess or know |A t ) but determined from the fundamen- 
tal Lagrangian as our model suggests. Really in order not to disagree drastically 
with the second law of thermodynamics the future should be shaped from the 
past and reflect the latter. However, there should not be - at least not much - ad- 
justment of the happenings at say time t in order to arrange something special 
simple happening in future. This means in or formalism that the by the Si future 
contributions determined |B t ) should according to second law better disappear 
quite from our formula for predicting probabilities for operator values, i.e. from 
a».4D or more generally ((87191 . 

Now, however, as we argued in foregoing section - section [8751 - reviewing 
previous articles working in the classical approximation it should be the state of a 
solution to the equations of motion in the early Big Bang time that dominates the 
selection of such a solution to be the realized one. The future on the other hand 
has only a small effect, if any, on choosing the true or realized solution. With the 
arguments to be given in section 15791 we argue for the effects of Si being even 
smaller in the future. Nevertheless we have if we talk exactly also effects of Si 
even in the future. Otherwise the hypothesis that the integral ( 18.8b defining |B t ) 
would be senseless since the e~ Sl -weighting is needed to suppress the integrand 
e~ Sl enough to make hope of a sensible practical convergence. 

However, we have in section 18.31 and will in section 18.91 argue for that S i 
varies much less in the future than in Big Bang era. 

It is now the purpose of the present section to use this only weak Si varia- 
tion with the fields in the future to argue for an approximation in density matrix 
terminology for the future part |B t ) of the functional integral. 

Let us indeed perform the following considerations for estimating the crude 
treatment of |B t ) which we shall use: 

a) Since S i has in practice only small non-trivial contributions in the future it is 
needed to involve contributions in the integral 



e iS ^cp = (B t |A t >. 



(8.29) 



dt dxL 



(8.30) 



from very large t > t'. 
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b) At these enormous t regions then at the end we get finally a rather restricted 
range of solutions. - we can think of classical solutions here, if we like - 

c) Now the solutions from the enormously late times under a) have to be devel- 
oped backward in time to the time t' say to deliver the state |Bf) (really we 
first get (B t < I4>) from equation l|8.8b ). 

d) Now we make the assumption that the system/world is sufficiently "ergodic" 
and the large times so large and so smeared out (also because of the smallness 
of the Li -effects) that we can take it that there is almost the same probability 
for finding the system in state |B t <) at any place in phase space allowed by 
the conserved quantum numbers of the theory practically valid in the future 
era. 

e) Ignoring for simplicity the conserved quantities we thus argued that with 
equal probability; equally distributed in phase space, we have that |B t <) will 
be any state. 

f) We can especially imagine that we have chosen a basis of wave packet states 
|w) in the field configuration space so that they fill smoothly the phase space 
- accessible without violating the conservation laws relevant -. Taking these 
to be - approximately - orthonormal (w|w') w 6ww' we clearly get for the 
average expectation of the projection operator 

P Bt , =|B t ')<B t ,| (8.31) 

the estimate 

av(|B t ,)(B t ,|) = 1 Y_ k)H * -1 (8.32) 

where N is the number of states in the basis 

|w),w = l,2,---,N. (8.33) 

That is to say we have argued for that our weak Si -influence in future com- 
bined with an assumed approximate ergodicity leads to that we can approxi- 
mate 

|Bt')<Bt'l«^l (8.34) 

in practice for all t' at least a bit later than the earliest Big Bang. 

The crude estimate that we could replace |B t ) (B t | by -^1 derived as formula 
(|B t )(B t | ~ ^qT] were based on that I_i were in practice small. 

8.5 Deriving a more usual probability formula 

We shall now make use of approximation d8.34b for the "future factor" in the 
functional integral in order to obtain an expression rewriting the formulas like 
< f873b , ((8~4l and 18.18} and 18.19} into expressions analogous to lOfl . 
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The calculation is in fact rather trivial, starting say from the most general of 
our postulated expressions 118.181 1: 



(0(<p| t , nit)) 



Je is 0(cp| t , TTjQgg 
Je iS ^cp 

(B t |Q((p|t, n| t )|At) 
(Bt|A t ) 

trivialstep (At|B t )(B t |0(<p| t , TT| t )|A t ) 



(A t |B t )(B t |A 



%§M (At|^lQ((p|t,n|t)|A t ) 
(Atl^l|A t ) 

(A t |Q((p|t,n|t)|At) 
(A t |A t ) 



(8.35) 



which is the completely usual quantum mechanical expression for the expecta- 
tion value of the operator 0((p| t , TT|t) in the wave functional state |A t ). 

With this expression we see that we should be allowed, as we anyway would 
expect, to use |A t ) as the quantum state of the universe. 

It should be noted though that our | At) is in principle calculable from the "the- 
ory" when as we shall of course, count also the Si -expression as part of the theory. 
In this way our model is widely more ambitious than usual quantum mechanics: 

We have - much like the Hartle-Hawking no boundary proposal - a func- 
tional integral l|8.7t delivering in principle the wave functional [At). In usual 
quantum mechanics the wave function is left for the experimental physicist to 
find out from his somewhat difficult job of preparing the state. In practice we 
would presumably have to let him be so helped by observation and arrangements 
under the preparation that we almost leave to him the usual job. We should, how- 
ever, have in mind that in preparing a state one will usually need to trust that 
some material is a rather pure chemical substance or that no disturbing cosmic 
radiation spoils the preparation. These kinds of trusts are usually based on some 
empirical experience which in turn makes use of that big assembles of pure sub- 
stances are easily /likely available and that generally cosmic ray has low intensity. 
Such trusts however, are at the very root connected with the starting state - the 
cosmology - of our world. But this starting state for practical purposes is in our 
model based on the activity of our l_i in early Big Bang times of the initial state of 
the universe. 

Thus it is even in the practical way of preparing a quantum state a lot of 
reference to our Si. 

If, however, somehow the universe develops into states where l_i is no longer 
negligible we should expect corrections to such an approximation (|B t )(B t | « 



-Li 

N — 



8 Future Dependent Initial Conditions from Imaginary Part in Lagrangian 101 

8.6 Time development and Si corrections to |B t ) 

From the definitions d8.7|l and d8.8|> of |A t ) and |B t ) it is trivial to derive the time 
development formulas for these Hilbert space vectors (say for t ' > t) 



I Af 



overtime— oo tot' 
iS 

over t to t ' 

W, t)|A t 



e' bttot 'A t [cp| t ]^(p 



(8.36) 



where ^(f, t) is the operator corresponding to the matrix (with columns and 
rows marked by cp| t configurations) 



Wit', 010 = 



over t to t ' with cp| t / — 0| t / and tp 1 1 — ^| 



(8.37) 



Similarly we have from | |8.8l l for t' > t again, first taking the complex conjugate 
of I 



(cp|t|B t ) 

and thus 

(cp|t|B t ) = 
which can be written 



Q ^Stto+OO 



beginning at cp 



e lb tt»+oo (cp| t ,|B t ')^cp 



|Bt)=^withL I ^-L I (t',t)+|B t ,). 



Here we used that e.g. 



Stt 



dt 



d 3 X(£ R +LC T ) 



(8.38) 



(8.39) 



(8.40) 



(8.41) 



where £ R and C\ are respectively the real and the imaginary parts of the La- 
grangian densities. So 



S* 
L t 



dt 



d 3 x(£ R -i£i), 



(8.42) 



and now restricting ourselves for {pedagogics/ simplicity} at first to boson fields 
we have (usually) that for them £ R and L\ are even order in the time derivatives 
which are under latticification 



3t<P{t, 



cp(t + At, x) — cp(t, x) 
At 



(8.43) 



Thus conceived as operators between the configuration at the two close by times 
t and t + At, i.e. with columns and rows marked by cp|t+At and cp| t we have e.g. 



(£ R + id) + =£ R -i£i 



(8.44) 
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because 

Cl = C R and Cj = Cj (8.45) 

and 

£r = C R and C\ = Cj. (8.46) 

In formula ([8.401 of course the meaning of the under symbol text in the expression 
^withLn-Li (f, t) + is that in addition to taking the Hermitian conjugation of 
^(t', t) as defined by the matrix representation l|8.37b one shall shift the sign for 
all occurrences of the Li-part of the Lagrangian or of the Li-part of the Lagrangian 
density. One should have in mind that it is easily seen that 

«r(t', tr 1 =% fiaiU ^-ult', (8-47) 

Especially the "usual" case of Li = means that &[t' , t) becomes unitary. This 
relation JOT) together with <f8~40l > and (|836) ensures that 



(B t |A t ) 



e tb_oot„+oo (8A8) 



can be true independent of the time t chosen on the left hand side. 

Since (6.1) represents a completely usual time development of the 'wave 
function' [At) we have of course analogously to the usual theory 

. d^O = H|At) (8 49) 

at 

where then H is the to the action 



S = S R + Si (8.50) 

corresponding Hamiltonian. As we saw under point a) in section 18.21 formula 
d8.20|) we can consider 

|A t ) (8.51) 

the wave function for the universe essentially. But really because of the normal- 
izing denominator in (18.20ft it is rather the normalized |A t ), namely 

|A t )norm = |A t )/ v / (At|A t ) (8.52) 

which is the true wave function. 

It is important to remeark that precisely because we now find that we shall 
use the normalized wave function rather than |A t ) itself we do not get as could 
be feared a lack of conservation of probability due to the non-unitarity of the time 
development. Have in mind that the to a non-real action corresponding Hamil- 
tonian H will not be Hermitean! But with the normalizaion comming from the 
(A L | A t ) in the denominator in l|8.20l l the total probability will anyway remain 
unity. This result matches nicely with the from the slightly different start evalu- 
ated (9.22) below. 
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Part II, Second Trial of Interpreta- 
tion 

8.7 Second Interpretation of the functional integral 

Usually one only uses the functional integral over a time interval to evaluate a 
transition matrix element from an initial time t| to a final time tf 



U(iM*lf),MH(ct>k)) 

where 



^fixedtime^ 



^fixedtime^i.^giSt^ot^*] (g 53) 



Sti totf — 



£(x)d 3 xdt (8.54) 



and the functional integral over 2>§ is restricted to cf-functions (field develop- 
ments, or paths) which at times t| and t f respectively coincides with 4>|i and cflf 
respectively. 

In the present article we, however, have the ambition of having the func- 
tional integral determine a priori not only the development with time but also 
say something about the initial conditions so that we a priori might ask for the 
probability of some dynamical variable O say having certain value O at a cer- 
tain time without imposing any initial conditions. In order to obtain a formula or 
proceedure or how to obtain such probabilities for what shall happen we have to 
assume such a formula. 

We therefore need some intuitive and phenomenological guess leading to 
such a formula/prescription. 

In order to propose such a formula in a sensible way we shall first consider a 
semiclassical approximation for our functional integral supposed to be connected 
with and describing the development of the Universe 

^4)e iS[ * ] (8.55) 

where we remember that in our model the S [<\>] is not as usual real but is allowed 
to be complex. 



8.7.1 Semiclassical approach 

For first orientation let us imagine that the imaginary part of the action S[§] is 
effectively small in the sense that we can obtain the most significant contributions 
to the functional integral by asking for saddle points for the real part Sr. That is 
we ask for field development solutions to the variational principle 

6S R = 0. (8.56) 

Without specifying the boundary conditions at t — > ±00 in our functional inte- 
gral there should be (essentially) one solution for any point in the (classical) phase 
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space of the field theory described. For the enumeration of the various develop- 
ment solutions ct) we could use the field and conjugate field configuration at any 
chosen moment of time, to say. However, now our hope and speculation is that 
the imaginary part should give a probability weight distribution over the set (~ 
phase space) of these classical solutions. 



8.7.2 A first but wrong thinking 

It is clear that we must make a definition of an expectation value for function(al) 
O say of the field development 4> so that if a single (semi) classical solution <J) so i 
comes to be highly weighted then this expectation value should be [cb sol ] . 
We might therefore at first think of 

{ >~ J^d>e*s[*] • ( ^ y) 

If really a single classical path contributed completely dominantly to both 
numerator and denominator, then indeed we would obtain that this proposal 
would obey 

(O) =0[* so ldom] (8.58) 

where cfWdom is this single dominant solution. 

It is however likely that if will be more realistic to imagine that there is a huge 
number of significant classical solutions 4> so i. But then appears the "problem" that 
in the expansion of the numerator functional integral into contributions from the 
various (semi) classical solutions <^ S oli- 



0<|>0[<|>]e lsl * ] = Y_ e iS[ * s ° lt] 0[c|5 so ,dv / d^ (8-59) 

4* sol i all the classical solutions 

the various contributions contribute with quite different signs or rather phases 
due to the appearance of the phase factor e lSR '* sou] . The proposal just put for- 
ward thus is not as it stands a usual average, it lacks the usual requirement of an 
average of being performed with a positive weight. Rather the summation over 
the contribution becomes a summation with random phases to a good approxi- 
mation. That means that if we classify in some ways the different solutions 4> so h 
into classes, then what would sum up when such classes are combined would 
be the squared contributions rather than the contributions themselves. In other 
words, if we define a contribution to 

i^e iSl * ] 0[4>] = Y Vde^~V s| * s ° u] en4>soK] (8.60) 



where 
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from a certain class of semi classical solution ^ then the quantities such as 

(8.62) 



Oe iS 9<$> 
obey approximately 



from class S 



Y_ v/dih V SI *^ ] 0[4> soH ] 



Oe iS ^4> 



from class 1 



Oe iS 3>(S> 



from class % 



from class Vi U % 



(8.63) 



However we do not have a similar addition formula for numerical values as 

2 



from class ¥1 

when they are not squared. However, of course, we do have 



(8.64) 



0e iS ^4> 



from class Vi 



Oe lS i?4) 



from class ^2 



(8.65) 



from class Ifi U ¥2 

but this relation has terms of typically rather random phases. 
8.7.3 Approaching a probability assumption 

If we take O to be a "projection operator" in the sense of being a functional of 4> 



should give the chance 



only taking the values and 1 then J" O e lS 2)§ 

from class 

for solutions in the class to pass through the configuration-class for which 
C?[cp] = 1. Because of the (random) phase and the lack of simple numerical ad- 
ditivity mentioned if the foregoing subsection we are driven to assume that the 
probability for cf) being in the [cj>] — 1 region must be given by the squared 
contributions 



0e iS 5>4> 



from class "ift 



(8.66) 



Calling the region in the space of 4>'s consisting of the 4>'s obeying 0[<\>] =1 
with our "project O", for region M, we get 



Prob(M) oc 



Me iS S>4> 



from class "ift 



(8.67) 



for restriction to the class 
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This means using the probability of the complementary set CM of M 

2 



Prob(CM) oc 
and the additivity (|8.63t 

e iS 2><k 

from class '< 

we derive 

Prob(M) 



CMe iS ^4> 



from class S 



(8.68) 



2 










e^cf) 






M 



from class 



2 






+ 




e iS ^c() 






CM 



from class S 



(8.69) 



e iS ^4> 



M 



from class e - 



e iS ^4> 



from class ^ 
-1 



4> so i i in M 



from class c < 



using random phases 



Y_ e^Vdet" 

<t>sol i 



-1 



from class S 



e-^'der 1 



<5> sol i in Mn% 



4>soi i in 



2S, 



det" 



(8.70) 



Here in principle of a classical approximation the e~ 2s 1 factor is much more im- 
portant than the "quantum correction" det -1 . Thus we would ignore the deter- 
minant det -1 factor in first approximation. 

Then we arrived to the picture here that the probability distribution over 
phase space - at some chosen time, that due to Liouville's theorem does not matter 
- is given by e~ 2Sl '* so1 ' where 4> so i is the classical field solution associated with the 
point in phase space for which e~ 2Sl '* so1 ^ shall be the probability density. 

8.7.4 About the effect of Si in the classical approximation 

To appreciate the just given probability density e~ 2Sl '* so1 ^ over phase space 

P(4>lt , nitoWlto^nito «e- 2s ' l *»' ] ^| t0 ,^n| t0 (8.7i) 

one should have in mind that in the classical approximation of the universe de- 
veloping along a solution 4> so i to the equations of motion 



6S R =0, 



(8.72) 
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the development is given quite uniquely by the equations of motion. The only 
place in which the imaginary part then comes in is in weighting with various 
probability densities the various "initial state data" (4>|t , n| to ) - i.e. the phase 
space point -. Once you know the initial state of the (sub)system considered the 
equation of motion determines everything in the classical approximation deter- 
mined by S R just described, the Si gets totally irrelevant. In other words it is only 
to know something about the "initial state" that Si has relevance. Here the usual 
terminology of "initial state" shall especially in our model not be taken too se- 
riously in as far as it with the word "initial" refers to a beginning moment, the 
Big Bang start say. No, as we just mentioned one can use any moment of time to 
for the description of the phase space describing the set of classical solutions 4>soh 
This t time does not have to be the first moment - even if such one should exist 
-. Rather we can use any moment of time as t . In the usual theory we would 
tend to use t being the initial moment and the state at this moment should then 
be one of very low entropy describing our start of universe state. However, in our 
model there is the rather unusual feature that the probability weight e~ 2Sl ^ so '^ is 
given via a functional Si [4> so i] depending on how the solution <J) so i behaves at all 
different times t and not only at to. Since we by the classical equations of motion 
can calculate the whole time development 4) so i from the fields and their conjugate 
(cf)|t , n| to ) at some chosen time t , we can of course also consider e~ 2Sl [ * sol] as a 
function of only the data at t , {4>\t , n| to ). 

So it is only by the fact that in our model e~ 2s R '* SDl1 is a rather simple function 
of cf) so i and thus because of the often chaotic development of the fields by the clas- 
sical equations of motion typically a complicated function(al) of the time t data. 
With a more usual model one might think of the initial state in a "first moment" 
to would be specified by some sort of cosmological model or no boundary con- 
dition. In this case the probability density should be rather simple in terms of the 
first moment data. The simplicity of e~ 2SR '* so1 ^ as functional of the cf) so i-behavior 
even at late times to some extend is extremely dangerous for our model showing 
observable effects not observed experimentally. Indeed an especially high prob- 
ability for initial states leading to a special sort of happenings today say would 
look as a hand of God effect seeking to arrange just this type of happenings to 
occur. In practise we never know the state of the universe totally at a moment 
of time. So there would usually be possibilities to adjust a bit the initial condi- 
tions. That ciould then in our model have happened in such a way that events or 
things to happen in the future gets arranged, if it can be done so as to organize 
especially big e ~ 2SRl * sol] i.e. an especially low Si. So a priori there would in our 
model be such "hand of God effects". In a later section we shall, however, invent 
or propose a possible explanation that could make the era of today be of very 
little significance for the value of Si so that in the first approximation it mainly 
the early time features of a solution that counts for its probability density 



e -2SiL*«*J ^ f(4> sol k early t i m es") = f(early time part of 4> sol ) (8.73) 
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8.7.5 Relation to earlier publications 

We have earlier published articles working in the classical approximation seeking 
to produce a model behind the second law of thermodynamics by assigning a 
probability P over the phase space of the Universe. It were also there the point 
that this probability density P in our model should depend in the same way on 
the state at all times. We already proposed that this P were obtained by imposing 
an imaginary part for the action Si. According to the above we clearly have 

Poce- 2Sl . (8.74) 



8.8 Suggestion of the quantum formula 



We already suggested above that if M denotes a subset of paths, e.g. those taking 
values in certain subset of 4>| t -configuration space in a moment of time t, then the 
probability for the true path being in M would be 



Prob(M) 



t M e iS @<t>\ 2 
Ue iS f?4>| 2 



(8.75) 



We imagine the paths to be described by the field <t) as function over M 4 , the 
Minkowski space. Thus we could use such an M to describe e.g. the project of the 
possible development $ to some subspace of configuration space Mi for a series 
of moments tj., i = 1 , 2, ■ ■ - , n. In fact then we would have 



M = | c() e {paths} 4>| ti e Mi for alii j . 
It would in this case be natural to think of the functional integral 



e iS $><t> 



(8.76) 



(8.77) 



M 



as a product of a series of functional integral associated with the various time 
intervals in the series of times — oo < ti < t2 < • • • < t n < oo. In fact let us define 



^t t tot t+1 (*ltt + ,,*lt t ) = 



BOUNDARY 
<t> ti and cf> | ti+ 1 kept 



e iS t t tot i+1 [*] 



Here 



■tt+i 



£(x)d 3 xdt 



(8.78) 



(8.79) 



(8.80) 



Stitot i+ i [<M = 

and remember that we here have the complex d(x), 

£(x)=£ R (x)+iCi(x). 
We then can write 

e iS 5>4> 

^ (3) 4>| ti ^t n tooo (4>loo , 4>lt J8m J4>k J^ tn _, tot n (*lt n , *lt n _, ) 

9 Mn _, • -e Ml ^ItJ^l^lt,,*!^ ) (8-81) 



M 
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where 0i[(t>|tj is the function 

We can also write this expression in language of a genuine operator product 

i to t n ^n— 1 ' ' ' 9 1 ^ co to 1 1 • 

(8.83) 

M 

where the 0t's are now conceived of as projection operators on the space of wave 
functionals characterized by being zero outside IVU, 

Here ip is a possible/general wave functional for the state of the universe, in the 
formula presented at the moment ti. 

In this operator formalism our probability formula takes the form 

Prob(M) = Prob(M, , M 2 , • • • , M n ) 

n tO OO 

= j (O.OD) 

n tO OO ^t n _, tot n • • • ^-oo tot, I 

Since we are anyway in the process of arguing along to just make an assumption 
about how to interpret in terms of probabilities for physical quantities of our com- 
plex action functional integral, we might immediately see that it would be sug- 
gestive to extend the validity of this formula for probabilities for field variables 
to also be valid for distributions in the conjugate fields n| ti or in combinations, 



Prob(C, e Mi , e> 2 e M 2 , • • • , O n e M„) 

^-L (8.86) 



l^tntOOO ' ' ' ^-OO tot! I 

Provided this proposal is not inconsistent to assume, we will assume it be- 
cause it would be quite reasonable to assume that the analogous formula to Il8.85b 
should be valid for any change for variables between 4> and IT in the formulation 
of our functional integral. 



8.8.1 Consistency and no need for boundary conditions 

It should be kept in mind that we expect that due to the presence of the imagi- 
nary part Si in the action S it is not needed to require any boundary conditions 
at t — > ±oo so that we basically can remove as not relevant the 4>|oo and 4>|-oo 
boundaries which one would at first have considered to be needed in the expres- 
sions ^_oo tot, (<f>lt i , cpl-oo ) and %, n to +oo (4>loo , 4>lt n )■ The imaginary part Si is in 
fact expected to weight various contributions so strongly different that whenever 
the by this weighting flavored component in cploo say is at all allowed by a poten- 
tial choice of boundary condition then that contribution will dominate so much 
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that all over contributions will be relatively negligible. So after taking the ratio 
for normalization such as d8.86|) the choice of the boundary conditions for cb|oo 
and c^l-oo becomes irrelevant. This irrelevance of the boundary conditions would 
indeed allow us to formally put in according to our convenience of calculation 
whatever boundaries we might like provided it does not precisely kill the bound- 
ary wave function component flavored by the Si. For instance we could put in 
at the infinity density matrices taken to be unity, since it does not matter anyway 
what we put and a unit matrix pi = 1 and = 1 would not suppress severely 
any state such as the flavored one(s). 

By this trick we could write our formula for probability 

Prob(C, 6M,,02 6 M2,-,0u€ M n ) 

= Tr(®t„ tooo Po n eM^ < ^t n _, tot n Po ll _,gM R _ 1 '"PdeM, ^-oo tot. 



^-00 tot, p o,eM, • ■ 'V-iGM^, ^t n _, tot n p o n eivw^ 



to oo 



/Tr(^ tn to oo^t n _, tow- ^-oo tot, tot, ■ • • io t n %[ tooo ) • (8-87) 

Here the reader should have in mind that because of the imaginary part in the 
action S = Sr + iSi the different development operators 



^ ti tot i+ ,(4>lt i+ ,,4>ltJ 



e iSI * ] ^4> (8.88) 

WITH BOUNDARIES 
^ and t(j | -[- ^ att^_|_j and t ^ respectively 



are in general not as usual unitary. Therefore it is quite important to distinguish, 

in general 



^t T itoti+1 + ^t7tot i+1 - (8-89) 



8.9 Practical Application Formulas 
8.9.1 Practical application philosophy 

Although in principle our theory is so much a theory of everything that it should 
even tell what really happens and not only what is allowed by the equations of 
motion, we must of course admit that even we knew the parameters of both Si 
and Sr it would be so exceedingly hard to calculate what really happens that 
can not do that. 

We are thus first of all interested in using some reasonable approximations 
to derive (in a spirit of a correspondence principle) some rules coinciding un- 
der practical conditions with the quantum mechanics (or quantum field theory 
rather) rules we usually use. 

Now as is to be explained in this section we can by means of requirements of 
monopoles and using the Standard Model gauge symmetries and homogeneity 
of the Lagrangian in the fermion fields argue that there will in the present era 
where Higgs particles are seldom and Standard model applicable be only very 
small effects of Si. We also seem to have justified to make the same assumption 
for very huge time spans in the future so that also until very far out in future the 
influence of Si is small. Even if we imagine that in the very long run Si selects 
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almost uniquely the state or rather development - as we used above to argue that 
the boundary conditions 4>|-oo and 4>| +00 were unimportant - then in the practical 
(i.e. rather near) future we would expect the far future determination to deliver 
under an ergodicity assumption an effective density matrix proportional to the 
unit matrix. 



8.9.2 Insertion of practical future treatment into interpretation formula 

The above suggestion for the treatment of the practical future to be equally likely 
in "all" (practical) states is implemented by replacing what is basically taking the 
place of a future density matrix Pf in our interpretation formula | |8.87t namely 

Pf ~^tttooo^t n tooo (8.90) 

by a normalized unit density matrix 

Pf « ~1, (8.91) 

where N is the dimension of the Hilbert space. (In practice N is infinite) So the 
interpretation formula becomes 

Prob(d e M,,e> 2 € M 2 ,--- ,O n e M n ) 

— TrfPc^gjvU^tn-, '"^V-iEMu-i ' ' ' ^Oi EM, ^-oo tot, 

®A p - . . . p - <9/t i 

*-oo tot, r O, £M| r On-l6M n _, ^tn-ltotn- 1 

/Tr(^ tn tooo ■ • • ^ to tr • • to t J • (8-92) 



where we used that 

O n eM n r O n £M 



p o„eM„ - p o. f M. • (8.93) 



8.9.3 Conditional Probability 

With formulas like ([8.87b or (|8.92l l we can easily form also conditional probabili- 
ties such as 

Prob(O p+1 e Mp+i,--- , O n G MnlOi 6 Mi , • • • ,O v G M p ) 
= Prob(0, eft,,-,0 p e^,-,O n 6 M n ) 
/Prob(0, e Mi,--- ,O p e M p ) . (8.94) 

In order to determine what happens if we know the wave function in some mo- 
ment. Let us as an example consider the idealized situation of a case in which we 
know - by preparation set up - the whole state of the universe of one moment 
of time. This we could imagine being described by taking a series of P . er ^. pro- 
jection of the same moment of time, the moment in which we suppose that we 
know the wave function. For consistency and for being able to take the limit of 
them being at same time - and therefore with an ill-determined algebraic order in 
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J8.87|l we must assume these same time projectors to commute. If we consider the 
situation that we already know that the system has all these Oi in the small re- 
gions. Mi because we know the wave function at their common time t com then we 
are after that discussing only the conditional probabilities with the set of relations 
Oi E Mi, i = 1 , • • ■ , p, taken as fixed. 

For simplicity let us consider the simple case that we just ask for if a variable 
O n at a later time being in M n or not then the conditional probability is in J8.87II 
form 

Prob(O n eMnl^I 

= Tr(^ t ^ to oo Po n eM n ^t com tot n Poi eM, Pe> 2 eM 2 ' ' ' ^a v eM. v ^-°° tot com 

^-00 tot com Po p eM p ' ' ' PozeMz^Oi EM, ^t com to t n P O n eM n %. n to oo ) (8.95) 
l Trf^tcomtotnPo, G M, ' ' ' PopGMp^-oo totcom^-oo tot com 
P • • • P - <9/^ 

r O p GM p r O|€M, ^tcomtOtn • 

Herein we can substitute 

W)kM = Po, GM, Po 2 6M 2 • • • PopGMp (8-96) 
and obtain using as usual for traces Tr(AB)=Tr(BA) 
Prob(O n e MM) 

— (ty I ^tl om to t n eMn^tL to oo ^tn tooo Po n £M n ^t„,m to t„ \A>) 

/ ( to oo ^t rom tooo KP) (W-oo to to t com W ) 

= {M ^tc„m tot n Pe> n 6M n ^t n to oo ^t n to oo e M n ^tcom to t n l"^) 

/W^tLtooo^toooW- 

(8.97) 

We may rewrite this expression in a suggestive way of the how it is modified 
relative to usual quantum mechanics by defining the final state density matrix 
from time t n 

Pffromtn =%{ tooo *tn tooo (8.98) 

in the following way 

(^l^t com totnPo n eMn Pffromtn P O n SMn ^tl m to t n l^) 



Prob(O n e MnW)) 



(^Wotn Pffromtn ^tLtOtn^) 



(8.99) 



Here ^t com tot n is really a non unitary S-matrix or development matrix for the 
time interval t com at which we have tj) to t n at which we look for O n . It is easy to 
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see that we could have replaced P© e ;& by a large set of commuting projections 
at a second common time t f destined to the single wave function |ib f ) and thus 
allowing to replace the series of projections put in place of ?o n GM n l^fX^f I- 
Then we get for the probability of the transition from |ib) to |ib f ) 

Prob (lab >|kb» = ^l^ mto tj^f)(^fip f f rom tj^f)(^fl^l mt0t jtb) 

' (m^tot.PffromW^totJ^ 

2 (^flPffromtJ^f) 



IM^WotJ%> 



(MW tcom to t n Pf from t n %\ om to t n M>) 

(8.100) 



Now we can compare this expression with the usual transition probability ex- 
pression when S is only real = Sr, 



Prob usual (|a|)f>|N)» = (ib|^t com totJtb f ) • (iW<, mto t» 

= to t n KW)| 2 (8.101) 

Denoting the transition operator 

S = (8-102) 

this means we have 



Prob(|^ f )|W) = |<tb f |S|tb)| 2 • ^jPffaomtJ^f) (g l03) 

compared to the usual expression 

Prob usual (N> f )|N>» = |(^ f |S|tb)| 2 . (8.104) 
The deviations are thus the following: 

1. With our imaginary part in S there is no longer unitality, i.e. 

S^S-'. (8.105) 
The transition S is calculated by the Feynman path integral with the full 

S=S R +iSi. (8.106) 

2. There is the extra wright factor (ib f |p f f romt ^ describing the effect of the hap- 
penings and the Si in the future of the "final measurement" |ib f ). 

3. There is the only on the initial state |ib) dependent "normalization factor" in 
the denominator 

WSWomwSlib) (8.107) 

This denominator is indeed a normalization factor normalizing the total prob- 
ability for reaching a complete set - an orthonormal basis - of final states 
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|i4>f ic), k = 1 ,2, ■ ■ ■ which we for simplicity choose as eigenstate of Pf f r0 mtn so 
that (|il>fkPf fromt n |iWk) g ets a diagonal matrix. Then namely 



k=1,2 



Y_ Prob(M> f ,ic>|N>» 

1 ,2,- 

1 ^- 2 

= T^t^ CUlA ' Z_ K^f,kl S l^)l (^f,klpffromt„l^f,k) 

\W r Pf fromtn 

= TTT^T- 1 ^TT\ 5~<^l St l^f,0<^f,lclPf fromt n l%,k)(^f,k|S|^) 

Pf fromtn b \W) ^ 

(8.108) 

which by using that the off diagonal elements of Pf f r0 mtn were chosen to be 
zero can be rewritten as a double sum - i.e. over both k and k' - 

^Prob (W>f,k)|W) W>|Stp ffiomtn SN)> 

k 

k,k' 

= (^|S t p ffromt nS|^) (8.109) 



Thus we see that this denominator just ensures that the total probability for 
all that can happen at time t n starting from |iL>) at t com becomes just one. 

8.9.4 Simplifying formula for conditional probability by approximating 
future 

We have already suggested that we should approximate 

Pf fromtn « ^1 (8-110) 

provided the future after t n is so that we can practically consider the Si -effects 
small or so much delayed into the extremely far future that our above ergodic- 
ity argument can be used. By such an approximation we remove the deviation 
number 2 above given by (iL> f |p f fromtn KW) because we approximate this matrix 
element (iWIPf fromtnKW) by a constant as a function of |il> f ). After this approxi- 
mation we get 

We should have in mind that S + S ^ 1 in general since with the imaginary part 
of the action the Hamiltonian will be non-hermitean and S non unitary. Thus the 
usual | (il) f |S|il)) | would by itself not deliver total probability for what comes out 
of |iL>) to be unity. Only after the division by the normalization (^jj> S+S |t4^) would 
it become normalized to unity. 
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8.10 Can we make an unsquared form? 

The formulas for the extraction of probabilities from our Feynman path integral 
with imaginary part of action Si also were derived by considerations of statistical 
addition with essentially random phases of various classical path. But our crucial 
formula, say (|8.87b , for probabilities is seemingly surprisingly complicated in as 
far as each projection operator occurs twice in the trace in the numerator. Even 
the simplest example of asking if some variable O at time t falls into the range M 
gets the expression 

Prob(0 € M) 

= Tr("^t tooo PoeM^-oo to t^-oo totPoeM'^ttooo )/Tr(^_oo tooo ^-00 tooo ) 

(8.112) 

containing the projection operator Po 6 m twice as factor in the expression. If we 
make the approximation of no Si -effects in the t to oo time range by taking 



Pf fromtn 



tjtto. (8.113) 



and approximating it by being proportional to the unit matrix then, however, 
the two projection operators come together and we could formally replace their 
product by just one of them. So in the case of the in this way approximated future 
we could write 

Prob(G e M) = Tr(P 0gtt ^ tattoo totW^-ootot^ioo tot)- ( 8 - 114 ) 

If O were a variable among the variables used as the path-description in the Feyn- 
man path integral the formula ^8.1121 would by functional integral be written 

ProbfO G M) = (8.H5) 
|JV S ^4>| 2 

Strictly speaking these Feynman integrals should be summed over all end of 
time configurations, but with significant Si presumably this summation would 
be dominated by a few "true" initial and states at ±oo and the summation would 
not be so important. 

So strictly speaking we have 

Hinit.final /initial ^OeM^^'k 

Prob(C e M) = — rr - T rj . (8.116) 



Linit,fin a llJe iS ^4)| 2 



8.11 The Higgs width broadening 

As an example of application of the Si -caused modification of the usual transition 
matrices we may consider the decay of a particle - which we for reasons to be 
explained below take to be the Weinberg-Salam Higgs particle - which has from 
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Si induced an imaginary term in the mass (or energy). Let us say take this term to 
have the effect of delivering a term being a positive constant number multiplied 
by — i in the Hamiltonian. In the Schrodinger equation 

= (8.117) 

such a term will cause the wave function xl> to decrease with time so that it will 
decay exponentially with time t. If the particle in addition decays "normally" into 
decay products, say bb as the Higgs particles do the exponential decay rate will 
be the sum r norma i + P s , of the S i -induced width 1$ : and the "normal" decay width 
^normal- Let us for simplicity take as an approximation that the real part of the 
mass is very large compared to both the "normal" and the Si -induced widths so 
that we can work effectively non relativistically with a resting Higgs particle. We 
can let it be produced in a short moment of time which is short compared to the 
inverse widths jX- and r 1 j while still allowing the particle may be considered 
at rest approximately. 

If we at first used the "usual" formula |(iL>f |S|-4>)| 2 for the decay process and 
calculate the total probability for the particle to decay into anything one will find 
that this probability is only r r,m ™ p s because the average lifetime has been re- 
duced by this factor, namely from — to T 1— = — . Since of course the usual 

J J 1 normal 1 normal T I S t 

particle with Is , = will decay into something with just probability unity, we 
thus need a normalization factor (il>|S^S|il>) to rescale the total probability to be 
(again) unity in our imaginary action theory. 

By Fourier transforming from time t to energy the Higgs decay time distri- 
bution we obtain in our model again a Breit-Wigner energy distribution 

p(E) = — j= r "i + r Sl (8 . 118) 



2n 



!E-m H i gg s) 2 + (^=^)' 



If indeed we effectively should have such an Si -induced imaginary part in 
the mass of the Higgs, then the Higgs-width could be made bigger than calculated 
in the usual width r norma i. This is an effect that might have been already seen in 
the LEP-collider provided one has indeed seen some Higgses in this accelerator. 
Indeed there has been found an excess of Higgs-like events with masses slightly 
below the established lower bound for the Higgs mass of 114 GeV/c. 



8.11.1 The effect of the (ipf pf from t„ |ipf ) suppression factor 

As an example of a (perhaps realistic) case of an effect of the factor (ibf | Pf irom tn |iL>f ) 
we could imagine that two particles are coming together organized to hit head on 
- say in a relative s-wave - able to potentially form two different resonances of 
which say one is a Higgs which as above is assumed to have an imaginary term 
in its mass. Now it is easy to see that the Pffr mt n (here t n is the moment the 
collision just formed one of the two resonances thought upon as physical objects) 
will have 

(iW HiggslPf fromtjlr'f Higgs) < (^falllPf from t n l^f all) (8.119) 
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where |ipf mggs) ar >d K|Jf ail) represent respectively the two possible resonances the 
Higgs and the alternative resonance. Compared to the usual calculation of the 
transition to one of the resonances - essentially the square of the coupling con- 
stant - the Higgs-resonance will occur with suppressed probability because of 
the (i|>| 

Pf fromtn |ij>) -factor in the formula ((8. 100b - Really if the collision were safely 
organized that the collision occurs because of s-wave impact preensured the total 
probability for one or the other of the two resonances to be formed would be with 
properly normalized probability 1 because of the (4>|S^S|4>) normalization factor. 
However, the effect of (*WIPf from v KM would be to increase the probability to 
form the alternative resonance while decreasing the formation of the Higgs. 

8.12 Approaching a more beautiful formulation 

Taking the regions in which may lie or not M as infinitesimally extended we 
would the formula for the probability density in the form 



Prob(C> = O ) 



^i,f |JBOUNDARY:i,f elS[ ' t ' 1P CigM^4 ) | 

£ tl f|J, 



piS[(t>]@4> I 
BOUND ARY:i',f ' e I 



,2 



K ri.fUBOUNDARY:!,^ 15 ^-^)^! (fi m)) 

V I f e iS[iM o,m\ 2 

^— i,f |J BOUND ARY:i',f "^1 

We may claim that this kind of formula the probability density for finding 
taking a value infinitesimally close to Oo is a bit unaestetic because of having the 
projection operator P^m or the equivalent 6(0 — Oo) occurring twice while one 
might have said it would be simpler to have just one 6(0 — Oo) or PoeM factor 
in the expression. 

We should now seek to reformulate our expression with these factors occur- 
ring twice into a simpler one with only single occurrence of PoeM or M0 — Oo)- 
To perform this hoped for derivation we first argue that for nonoverlapping 0- 
value regions Mi and M2 



P 0e M,e iSI * ] ^4> 

BOUND ARY:i,f 



P 0e M 2 e iS[ * ] ^4) 

BOUND ARY:i,f 



i,f 

» for Mt n M 2 = . (8.121) 

This is argued for by maintaining that giving a different value at time t very 
typically by "butterfly effect" - Lyapunov exponent - will cause very different 
states f and i at =Foo respectively. If the two factors in (18.1211 have very different 
final f and initial i states dominate at the boundaries and even random phases 
the total sum is indeed much smaller than what one would obtain if Mi and 
M.2 were taken to be the same region Mi = M2 = M. If we now use the zero 
expression (|8. 1211 1 by adding such terms into the numerator and analogously in 
the denominator of J8.120I I we can formulate this expression for the probability of 
being in M into an expression involving a summation over the value or region 
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for O in one of the occurrences 
Prob(C = C' 2) ) = 



L 

. i,f 



BOUND ARY:i,f 



BOUND ARY:i,f 

b(0 -O { 2) )e iSW ®§ 



(8.122) 



BOUNDARY: i ' , f ' 



e iS[ * ] ^4) 



BOUND ARY:i',f ' 



e iSl * ] ^4) 



But now obviously we have 



dO [ ]) 6(O-O ) = 1 



(8.123) 



and thus we get 

Prob(e>-e><- 



21- 



BOUNDARY: 



if e iS| *)^4>)* J 



BOUNDARY: 



i f &{O-O [ 2) )e iS[ * ] ^c() 



Li' f ' (J: 



BOUNDARY: i ' , f ' 



BOUNDARY 



., f ,e iS I*^4) 



(8.124) 



In this expression we have achieved to have 5(0 — O } only occurring once as 
factor. We could therefore trivially extract from it an expression for the average of 
the O-variable 



(O) 



{ 2] Fmb(0-O l 2) )dO { o 



(2) 



Lt, f (J 



BOUND ARY:i,f 



J. 



BOUNDARY: i ' , f ' 



3 iS[ct>] 



BOUNDARY 



. f Oe iSl * ] ^4) 



BOUNDARY: 



v f ,e ts l*]^4) 



(8.125) 



If we could somehow remove the after all identical complex conjugate functional 
integrals 



and 



BOUND ARY:i,f 



BOUNDARY:i',f ' 



3 IS[«] 



(8.126) 



(8.127) 



only deviating by the dummy initial and final state designations respectively (i, f ) 
and (i', f '), then we could achieve the simple expression l|8.57l l. But in order to ar- 
gue for such removal being possible we would have to speculate say that some 
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- we could say the true - boundary condition combination for the functional in- 
tegrals (|8. 1261 18.127b completely dominates. This is actually not at all unrealistic 
since indeed the Si will tend to very few paths dominate. In such a case of dom- 
inance we would have a set of dominant (f,i). Presumably to make the chance 
that there should be such dominance we should allow ourselves to be satisfied 
with a linear combinations of i-state and of f-states to dominate. But now if in- 
deed we could do that and call these linear combinations (fdom, i-dom)/ then we 
could approximate 

i n \ _ J BOUNDARY f dum ,i dom ^ , R . 

J BOUNDARY f dom ,i dom e V 

Now we would like not to have the occurrence in this expression of the rather 
special states (fdom, i-dom) - However, these dominant boundary conditions are pre- 
cisely the dominant boundary conditions for the denominator integral, because it 
were really just the complex conjugate for the latter for which we looked for the 
dominant boundaries. 

So if we let the boundaries free then at least the denominator should become 
dominantly just as if we had used the boundaries (fdom, i-dom) - It even seems that 
because of the smoothness and boundedness of the variable as functional of <j) 
the dominant boundaries (i, f ) should not be much changed by inserting an extra 
factor so that also by letting the boundaries free in the numerator functional 
integral J" 0[4)]e lS '*'£^(f) would not change much the dominant boundaries from 
those of the same integral with the 0-factor removed. But the removal of this 
leads to the denominator functional integral, for which we already saw that the 
dominating boundary behavior were (fdom> i-dom) - Thus we have argued that we 
can rewrite l!8.128l l into 

(0) = J f (8.129) 

where it is understood that the boundaries for t — > ±oo are "free". Then we sug- 
gested they would automatically go to be dominated by (fdom, i-dom) thus fitting 
on to the formulas with double occurrence of 5(0 — 0o)'s. 

The argumentation that the factor does not matter for the dominant be- 
havior at ±oo may sound almost contradictory to our assumption using the "but- 
terfly effect" to derive the rapid variation of ||8.4| | which meant that an insertion 
of 6 (0 — Co) would drastically change behavior, including that of the boundary. 

It is, however, not totally unreasonable that a sharp function 6 ( — Oo ) which 
is zero in most places could modify the boundary conditions, while a smooth one 
0, almost never zero would not modify them. Basically we hope indeed for that 
the Si -caused weighting is so severely restricting the set of significant paths, that 
it practically means that a single path, "the realized path" is selected. In this case 
the insertion of the factor into the functional integral would just multiply it 
by the value of on "realized path". If you however insert 6(0 — Oo) and it 
as most likely the case Oq is not the value of on the realized path then we 
kill by the zero-value of 6(0 — Oo) at the realized path would totally kill the 
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dominant contribution. Then of course the possibility for a completely different 
path is opened and the orthogonality used in J8.12H gets realistic. 

As conclusion of the just delivered derivation of l!8.129l l we see that the start- 
ing point in the beginning the articles is indeed consistent under the suggested 
approximations with the forms derived from the semiclassical start. 

8.13 The monopole argument for suppressing the Si in the 
Standard Model 

We have already above in section l83l argued that due to the material in the present 
era, and the future too, being either massless or protected from decay by in prac- 
tice conserved quantum numbers and due to weakness of the interactions the 
contribution to S i from these eras must be rather trivial. 

It were also for the above argument important that the non-zero mass parti- 
cles were non-relativistic in these eras. That above argument may, however, not 
be sufficient for explaining that no effect of our I_i or Si would have been seen 
so far. We have indeed had several high energy accelerators such as ISR (inter- 
secting storage ring at CERN) in which massive particles - such as protons - have 
been brought to run for days with relativistic speeds. That means that they would 
during this running in the storage rings say have had eigentime contributions sig- 
nificantly lower than the coordinate time or rather the time on earth. This would 
presumably easily have given significant contributions to Si which going to the 
exponent would suppress - or priori perhaps enhance - the probability of devel- 
opments, solutions, to equations of motion, leading to the running of such storage 
rings. Since the protons have not already been made to run around dominantly 
relativistically we should deduce that most likely the storage rings would lead to 
increasing Si and thus lowering of the probability weight. Thus one would expect 
that the initial conditions should have been so adjusted as to prevent funding for 
this kind of accelerators, at least for them running long time. Contrary to Higgs 
producing accelerators which have so far not been able to work on big scale (may 
be L.E.R produced a few Higgses for a short time) the accelerators with relativistic 
massive particles have seemingly worked without especially bad luck. In order 
to rescue our model it seems therefore needed to invent a crutch for it of the type 
that there are actually no I_i -contributions involving the particles for which the 
massive relativistically running accelerators have been realized. We have actu- 
ally two mechanisms to offer which at the end can argue away our I_i or Si effects 
for all the hitherto humanly produced or found particles, leaving the hopes for 
finding observable effects - bad luck for accelerators, mysterious broadening of 
resonance peaks - to experiments involving the Higgs particle or particles out- 
side the Standard Model. The point is indeed that we shall argue away the effects 
of Si for gauge particles and for Fermions (coupled to them). 

The suppression rules to be argue for are: 

1) Supposing the existence of monopoles we deduce that the corresponding full 
gauge coupling constants must be real, basically as a consequence of the Dirac 
relation. 
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2) For fields which like the Fermion fields in renormalizable theories occur ho- 
mogeneously in the Lagrangian density Cr + \C\ this Lagrangian density can 
be shown to be zero by inserting the equations of motion. 



8.13.1 Spelling out the suppression rules 

Spelling out a bit the suppression rules let us for the monopole based argu- 
ment remind the reader that although we consider a complex Lagrangian den- 
sity £ R + iCi for instance the electric and magnetic fields and the four potential 
Au.(x) for electrodynamics are still real as usual. Now if we have fundamental 
monopoles there must exist corresponding Dirac strings which, however, as is 
well known must be unphysical. The explicit flux in the Dirac string must to 
have the Dirac string unobservable - to be unphysical - be compensated by an 
at the string singular behavior of the four potential around the Dirac string. 
The singular flux to compensate the extra flux in the Dirac string can, however, 
only be real since the A^-field is real and it is given by a curve integral § A^dx^ 
around the Dirac string. Now as is well known the fluxes mentioned equal the 
monopole charges. Thus the monopole charge g must be real. But then the Dirac 
relation 



eg = 27m, n e Z (8.130) 

tells that also the electric charge e must be real. Now, however, in the formalism 
with the electric charge absorbed into the four potential A^ = eA,^ the coefficient 
on the F^-term in the Lagrangian density is — so that the pure electromag- 
netic, kinetic, Lagrangian density 



becomes totally real. I.e. 



1 F 2 ^ (8.131) 

puree.m. 



~4e 2 ^ 



£ljpuree.m. = 0. (8.132) 

We may skip or postpone a similar argument for non-abelian, Yang Mills, theories 
to another paper, but really you may just think of some abelian subgroup and 
make use of gauge invariance. 

Concerning the rule 2) for homogeneously occurring fields, such as the Fer- 
mion fields in renormalizable theories the trick is to use the equations of motion. 
For example the part of the Lagrangian density Cr + iCi involving a Fermion 
field -ty is of the form £p = Z • \f>(iJZ)— m)t|) where Z is a constant and D^ the co- 
variant derivative and of course D = y^D^. This Fermionic Lagrangian density 
is homogeneous of rank two in the Fermion field The Euler-Lagrange equa- 
tions, the equations of motion for the Fermion fields are derived from functional 
differentiation w.r.t. the field tJj> 



(8.133) 
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and end up giving equations of motion of the form 



9£f 
IhjT 







(8.134) 



or 



3£f 







(8.135) 



(really these forms are only trustable modulo total divergences but that is enough) 
leading as is well known to 



But now it is a general rule that a homogeneous expression, £f say, can be recov- 
ered from its partial derivatives 



where rank is in the present case rank= 2. Such a recovering for homogeneous La- 
grangian densities, however, means that the Lagrangian density - at least modulo 
total divergences - can be expressed by the equation of motions, which are zero, 
if obeyed. But then at least in the classical approximation the Lagrangian density 
is zero at least modulo total divergences. This means that the total Sr + iSi con- 
tribution from the just discussed homogeneous terms end up zero. Especially the 
imaginary part also ends up zero, although its form does not have to be zero. It is 
only insertion of equations of motion that makes it zero. 

8.14 Conclusion 

We have put up a formalism for a non-unitary model based on extending the La- 
grangian and thereby the action to be complex by allowing complex coefficients 
in the Lagrangian density £r + iCi. 

We used two starting points for how to extract probabilities and expectation 
values from the Feynman path way integral in our ambitious model that shall 
even be able to tell what really happens rather than just the equations of motion. 
The first were the interpretation that an operator 0(t) should have the expecta- 
tion value 



aI>(iJD-m) =0 



(8.136) 



or 



(iJD-m)oI> = 0. 



(8.137) 




(8.138) 



(0) = 



j-0(t)e iSl * ] ^4> 
Je iS [*1^4> 



(8.139) 



but this expression is a bit dangerous in as far as it is a priori not guaranteed to 
be real even though the quantity 0(t) is real. The second approach would rather 
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have a series of projections onto small regions M| for operator £\(ti) denoted 
Pej.gjft. inserted into the functional integral but then this integral is numerically 
squared for any combination (i, f ) of boundary behaviors at respectively — oo and 
+00 times. That is to say that the insertions are to be performed into the integral 
Je iS[ * ] ^4> so as to replace the latter by J* Yli ^o^m.^^^^ j ust as m tne nrst 
approach, but then one forms the numerical square summed over the initial i and 
final f behaviors 



The probability distribution is then obtained by inserting the projection op- 
erators into both factors in (|8. 1401 1 and then as normalization divide the (|8.140b 
having these insertion with 08.1401 1 not having the insertions. 

Under some suggestive assumptions we argued that the two approaches ap- 
proximately will agree with each other. The most important formula derived is 
presumably the formula to replace usual unitary S-matrix or 1^ -matrix transition 
between two moments in time in our model. This formula turns out for transition 
an initial state to a final |i|;f ) to be 



We used that to derive the broadening in our model of the Higgs-width. 

As an outlook we may mention some of the expectations of our model used 
in a more classical language in our earlier publications: If the Higgs - especially 
freely running Higgses — decrease significantly the probabilty (7.21) then the 
initial state should be organized so that Higgs production be largely avoided. 
This would actually make the prediction that some how or the other an accident 
will happen and the LHC-accelerator will be prevented from comming to full 
energy and luminosoty. 
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Abstract. We discuss possible origins of Multiple Point Principle. Inspired by results from 
finite temperature lattice gauge theory we conjecture possible physical reasons which 
could lead to the Principle 

9.1 Introduction 

One of many plagues of the Standard Model as we know it is the excessive num- 
ber of parameters. It is usually understood that it is a manifestation of Stan- 
dard Model being some effective theory for the "true" underlying model, such 
as GUTs. One very promising way of explaining how parameters get fine-tuned 
is so-called Mulptiple Point Principle! 6 [ introduced and advocated by Nielsen et 
al. As it is widely discussed in the very same proceeding we will not focus on its 
details and implications but only mention the essential for our conjecture prop- 
erty. The core statement of MPP is the assumption that many degenerate vacua 
exist and that they are separated by first order phase transition. Couplings of the 
model do self-tune themself to the critical lines, or, in existence of many vacua - to 
multi-critical points. Here we suggest how this situation may indeed be realized 
in nature and discuss the implications. 

9.2 Dimensional Reduction - a detour 

Before we turn to the main subject it is important to discuss certain results which 
come from finite temperature gauge theory. We need to underline that while these 
results were largely obtained using Lattice Gauge Theory the mechanism is rele- 
vant to continuum and lattice here is just a way to non-perturbatively study phase 
diagram. For simplicity we will consider 2+1-dimensional spacetime, though the 
results are similar in 3+1 dimensions (nobody ever analyzed 4+1 or more). Our 
model, which we will regard as "physics" would be pure Yang-Mills SU(3). It can 
be formulated in terms of closed loops on the (hyper)cubic lattice. We will de- 
note its extentions in space and in time directions by L s and I_o correspondingly. 
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The lattice is periodic in all directions, however differently. Spacial extent is finite 
only for practical purposes and generally is sent to infinity. Time extent behaves 
as usually in finite temperature field theory and is compact, with its length being 
equal the inverse temperature 

L a=1/T (9.1) 

On each link of it we put a matrix U(x; \±) which belongs to the SU(3) group and 
is linked to the gauge fields A by 

U(x;h) = expig A(x) (9.2) 

where go is the bare strong coupling constant. 

S^ = -4- Y- Tr(U(x;^)U(x+ ^^JUtx+^^-'Ufx;^/)- 1 +H.c.) (9.3) 

yields usual square of the field strength tensor in naive continuum limit. Phe- 
nomena, symmetries and phases which occur in this model we will consider to 
be "physical". 

It has a deconfinement phase transition separating confined and deconfined 
phase with Z(3) being relevant symmetry. If we go to high temperatures, about 
twice the critical one we can do a following trick, known as dimensional reduc- 
tion. As you can see from Eq l9.1l the time extent shrinks drastically. The theory 
becomes perturbative and time-like degrees of freedom may be integrated out. 
This has no relation to lattice and integrating is often done in continuum. First 
we set time-like gauge field to be independent of the imaginary time coordinate 

A (x 0) x) = A (x) (9.4) 

and then eliminate the remnants of the gauge freedom by adding the so-called 
Landau constraint 

2 

^^[A i (x)-A i (x-at)]=0 (9.5) 

Together with l|9.4|l this condition is called Static Time- Averaged Landau Gauge. 
It is essential that all gauge degrees of freedom are fixed and what is left is, or 
should be, physical. 

Few Feynman diagrams later (for details reader may consult [1]) we get an 
effective theory, formulated in two dimensions with an extra field which lives on 
sites. The resulting action depends on the parameters of the original theory Lo 
and go. We make a change of variables 

A (x) = <Mxh/{£ (9.6) 

V L o 

to make effective model look in a more familiar way. 

Now we can express S^ f f as the function of this field and the two-dimensional 
gauge field U. It consists of the following three parts: 
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• pure gauge part S w , which we had already in "naive" reduction 

S w = ^ L oZ Tr ( ] - l^U(x;1)U(x + r,l)U(x + l;t)-'U(x;l)-^ (9.7) 

• gauge covariant kinetic term, depending both on U and cfj, obtained by ex- 
panding the Sly to the second order in Ao 

Su,* = XI Tr ( Di ( u ^' x " < 9 - 8 ) 

x i=l,2 

with Di being covariant derivative 

D t (U) <\>{x) = U(x;t)<J>(x + cnMxjtr 1 - 4>(x) (9.9) 

• self-interaction term S^, often called "Higgs potential" 

=h 2 ^Trcf) 2 (x)+h4^(Tr^ 2 (x)) 2 (9.10) 

X X 

where couplings are fixed by the parameters of the original model. 

Now we have an effective action for our original theory which lives in lower- 
dimensional space and reproduced very accurately results of the higher-dimensional 
model. To understand it better one can make the couplings free and study the 
model per se. Then, if you do a non-perturbative analysis of the effective the- 
orytwo unexpected results appear. The theory has new phase transition, which 
was not there in the original model with relevant symmetry being time-reflection, 
here manifested as cf — > — cp. It can be spontaneusly broken and Higgs field gets 
vacuum expectation value. The transition is of a strong first order. Further impor- 
tant property of this transition is that original coupling s are pretty much on the 
transition line. Moreover, they turn out to be in a wrong phase. This suggests that 
up to numerical ambiguity and errors from the perturbation theory in the course 
of dimensional reduction - they may be exactly on the transition line. This transi- 
tion is of course not a physical one, as it is not present in the underlying theory 
which we define as "physics". Now lets us go to the conjecture itself. 

9.3 The Conjecture 

While the above described scenario has technically nothing to do with nothying 
but pure gluodynamics we can reformulate it. First we note that there is noth- 
ing special in the definition of temperature and instead of talking about high- 
temperatures effective theory we will view the model as zero-temperature one 
with one highly compact dimension. This way we can translate forementioned 
results into something more generic. 

• Gauge theory in D dimensions of which one is highly compactified may be 
rewritten as effective theory in D — 1 dimensions. 
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• Effective theory will be formulated via original variables and an additional 
field in the algebra of the gauge group with new couplings. 

• There will be a first order phase transition with one of the phases being un- 
physical. 

• For an observer which lives in D— 1 dimensional world and somehow guessed 
the effective theory it will be a "miracle" that couplings self-tune themself to 
the critical line. 

This scenario has been extensively tested in 3 dimensions with SU(3) being a 
gauge group and in 4 dimensions with SU(2)|1|. It is a bold move to generalize 
it to arbitrary number of dimensions as renormalizability, which differs from di- 
mension to dimension, may change the situation. However forementioned cases 
are also very different so one may hope it is a universal behaviour. 

If so, and the number of compact dimensions is more than one it may give 
rise to a number of phases with multi-critical points then working as attractors 
for effective couplings. 

Then it can be the physical reason behind Multiple Point Principle. The dif- 
ference between our conjecture and the MPP is of course that only one phase 
is actually physical and the remaining ones are the artefacts of the reduction of 
the compactified space. Once one has a candidate theory and mapped its phase 
diagram one can get extra information about how the couplings are fixed. How- 
ever one should keep in mind that the dynamics in these phases is completely 
unphysical and cannot be used for phenomenological purposes. 
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